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PREFACE 


The  work  reported  herein  was  performed  for  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  under  Purchase  Order  DACA39-82-M-0074  with 
the  University  of  Maryland  during  the  period  May  1982  to  June  1983.  The  work 
was  funded  by  the  Office,  Chief  of  Engineers,  U.  S.  Army,  under  the  FY  82 
RDTE  Program,  Project:  4A161102AT22,  Task  AO,  Work  Unit  009,  "Methodology  for 
Considering  Material  Variability  in  Pavement  Design."  OCE  Project  Monitor 
was  Mr.  S.  S.  Gillespie.  The  work  was  conducted  and  report  prepared  by  Drs. 

M.  W.  Witczak  and  J.  Uzan  and  Mr.  M.  Johnson. 

Dr.  W.  R.  Barker,  Pavement  Systems  Division  (PSD),  Geotechnical 
Laboratory  (GL) ,  WES,  was  the  Principal  Investigator.  The  study  was  conducted 
under  the  general  supervision  of  Dr.  T.  D.  White,  Chief,  PSD,  and  Dr.  W.  F. 
Marcuson  III,  Chief,  GL. 

Commander  and  Director  of  WES  during  the  period  of  this  study  was  COL 
Tilford  C.  Creel,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 


1 


CONTENTS 


PREFACE  . 

VOLUME  Is  STATE  OF  THE  ART  VARIABILITY  OF  AIRFIELD 
PAVEMENT  MATERIALS 

VOLUME  II:  MATHEMATICAL  FORMULATION  OF  RELIABILITY  MODELS 
UTILIZED  IN  RIGID  PAVEMENT  STUDIES 

VOLUME  III:  PROBABILISTIC  ANALYSIS  OF  RIGID  AIRFIELD  DESIGN  BY 
THE  WESTERGAARD  FREE  EDGE  THEORY 

VOLUME  IV:  PROBABILISTIC  ANALYSIS  OF  RIGID  AIRFIELD  DESIGN  BY 
ELASTIC  LAYERED  THEORY 


VOLUME  I 


STATE  OF  THE  ART  VARIABILITY  OF 
AIRFIELD  PAVEMENT  MATERIALS 


TABLE  OF  CONTENTS 


LIST  OF  TABLES 
Chapter  1: 
Chapter  2: 


Chapter  3: 
REFERENCES 


.  ii 

INTRODUCTION .  1 

MATERIAL  VARIABILITY  .  3 

Portland  Cement  Concrete  Materials  .  3 

Thickness . 3 

Modulus  of  Elasticity . 3 

Poisson's  Ratio  .  6 

Modulus  of  Rupture .  11 

Asphalt  Concrete  Materials  .  11 

Thickness .  11 

Dynamic  Modulus  .  11 

Flexural  Stiffness  .  15 

Poisson's  Ratio  .  16 

Cement-Treated  Materials  .  16 

Thickness .  16 


Modulus  of  Elasticity  .  16 

Subgrade  Materials  .  16 


Modulus  of  SubgTade  Reaction  .  16 

Resilient  Modulus  .  20 

California  Bearing  Ratio  .  .  .  20 

SUMMARY  AND  RECOMMENDATIONS  .  23 

.  25 


i 


LIST  OF  TABLES 


Table  No. 


Page 


1.  Thickness  of  Portland  Cement  Concrete  Pavements .  4 


2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 

11. 

12. 

13. 

14. 


Summary  of  Statistical  Results  on  Thickness  of 

Pavement . . . 

Summary  of  Test  Results  for  Cored  Specimens 

from  PCC  Airfield  Pavements  . 

Individual  Sonic  E  and  u  Values  -  50  series . 

Individual  Sonic  E  and  u  Values  -  60  series . 

Individual  Sonic  E  and  y  Values  -  70  series . 

Average  7-day  Flexural  Strength  and  Standard 

Deviations,  Portland  Cement  Concrete  Pavement.  .  . 

Average  28-day  Flexural  Strength  and  Standard 

Deviations,  Portland  Cement  Concrete  Pavement.  .  . 

Lab  Variability  -  Asphalt  Concrete . 

Summary  of  Test  Results  for  Cored  Specimens  from 
Blackbase  Highway  Projects  in  Texas.  ....... 

Flexural  Stiffness  Measurements  on  Field  Samples 
of  Asphaltic  Concrete  Using  Pulse  Loading  Method  . 

Summary  of  Properties  for  Cement -Treated  Bases  for 
Airfield  and  Highway  Pavements  . 

Summary  of  Resilient  Moduli  for  Subgrade  Soils  for 
Airports . . . 


5 

7 

8 
9 

10 

12 

13 

14 

17 

18 

19 

21 


Subgrade  CBR  Variability 


22 


15.  Variability  Recommendations 


24 


ii 


Chapter  1 
INTRODUCTION 

At  present,  pavement  design  methodologies  are  deterministic, 
that  is,  a  unique  pavement  system  is  designed  for  the  set  of  input  var¬ 
iables  having  unique  values.  The  variations  of  these  design  parameters 
have  been  taken  into  account  eitheT  by  setting  limiting  values  on  the 
parameters  or  by  using  a  factor  of  safety.  Such  an  approach,  however, 
leaves  the  engineer  without  any  true  quantitative  estimate  of  overall 
design  reliability.  A  probabilistic  design  method  may  therefore  be  devel¬ 
oped  and  used  in  order  to  take  material  variabilities  into  account  within 
a  reliability  framework.  The  use  of  a  probabilistic  methodology  neces¬ 
sitates  the  determination  of  not  only  the  mean  (average)  value  of  each 
parameter;  but  a  measure  of  the  variability  of  the  specific  parameter 
in  question  (i.e.,  standard  deviation,  variance,  coefficient  of  variation). 

The  purpose  of  this  report  is  to  give  a  general  summary  of  values 
or  ranges  of  values,  found  from  a  literature  search,  that  describe  mater¬ 
ial  variabilities.  This  report  summarizes  a  state-of-the-art  review  in 
variabilities  of  material  properties  which  are  used  in  pavement  design. 

This  study  specifically  concentrated  on  normal  design  input  variables 
coimonly  used  in  rigid  pavement  airfield  studies.  As  such,  the  wealth 
of  variability  information  concerning  such  non-design  parameters  (i.e., 
soil  classification,  gradation,  component  variability  analyses,  etc.) 
was  not  included. 

It  should  be  noted  that,  in  some  instances,  the  standard  deviation 
is  the  given  measure  of  variability.  Since  the  coefficient  of  variation 
may  not  always  be  constant,  the  standard  deviation  may  give  a  better 


indication  of  the  variability  of  the  naterial.  If  a  standard  deviation 

is  given  and  the  coefficient  of  variation  is  required,  it  may  be  obtained 

o 

by  dividing  the  standard  deviation  by  the  mean  (CV  *  -3-  x  100%). 

Chapter  2  contains  brief  discussions  of  material  properties  and 
their  variation  while  Chapter  3  contains  the  summary  and  recommendations 
of  values  to  be  used  within  a  probabilistic  design  method  foT  rigid 
airfield  pavements. 
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Chapter  2 


MATERIAL  VARIABILITY 

Well  defined  variability  information  is  more  prevalent  from  highway 
pavements  than  for  airfield  pavements.  However,  there  is  little  reason 
to  believe  that  there  are  major  differences  in  the  variability  of  the 
two  types  of  pavements.  Variabilities,  therefore,  have  been  collected 
for  both  highways  and  airfields  where  there  is  insufficient  airfield 
information. 

PORTLAND  CEMENT  CONCRETE  MATERIALS 
Thickness 

There  is  very  little  variation  in  the  thickness  of  Portland  Cement  concrete 
pavements.  Table  1  shows  information  taken  from  several  highway  and  air¬ 
field  pavements.  The  coefficients  of  variation  are  generally  below  3% 
but  go  as  high  as  8%.  Table  2  shows  information  from  the  Louisiana 
Department  of  Highways.  Here,  the  coefficients  of  variation  are  all  less 
than  5%.  Information  could  only  be  found  for  two  airfield  pavement 
sections  (Table  1).  All  of  the  information  gathered  was  for  pavement 
thicknesses  of  less  than  15  inches.  An  insignificant  amount  of  data 
was  found  for  pavements  of  larger  thicknesses  and  therefore  no  conclusions 
may  be  drawn. 

Modulus  of  Elasticity 

There  was  a  sufficient  amount  of  information  available  on  the  modulus 
of  elasticity  for  airfields.  A  summary  of  the  mean  values  and  the 
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Table  1.  Thickness  of  Portland  Cement  Concrete  Pavements 
(Kennedy,  Hudson  and  McCullough 


Pavement  Thickness* 


Project 

I dent i- 

CV 

Design 

fication 

(%) 

(inches) 

2 -A 

- 1 - 

ATP 

38 

8.3 

2  .6 

8.0 

2-E 

ATP 

50 

8.3 

8.0 

17-B 

ATP 

SO 

8.2 

2  .6 

8.0 

Cluster 

10 

7.8 

1  .2 

8.0 

ATP 

47 

8.2 

2.6 

8.0 

17-M 

Cluster  1 

7 

■ 

8.0 

Cluster  2 

8 

8.0 

18-N 

ATP 

9 

— 

8.8 

mm 

8.0 

18-0 

ATP 

24 

9.5 

4.7 

9.0 

19 -A 

ATP 

34 

8.2 

8.0 

ATP 

31 

8.2 

mm 

8.0 

19-B 

Cluster  1 

10 

7.6 

8.0 

Cluster  2 

J 

9 

7.6 

WM 

8.0 

Weighted  Average 
excluding  clusters 

- 

2.9 

- 

Variation 

Limits 

-  ' 

2. 5-4. 7 

- 

Range 

- 

2.2 

- 

O' Hare 

Rwy  9R-27L 

ATP 

12 

10.9 

8.3 

10.0 

O'Hare 

Rwy  4R-22L 

ATP 

10 

14.8 

3.3 

14.0 

- - 

Weighted  Average 

_ 

6.0 

- 

Variation 

Limits 

3. 3-8. 3 

- 

UM 

5.0 

- 

*  Thickness  determined  by  measuring  height  of  core  in  laboratory. 


**  Along-the-pavement  or  cluster  samples  from  thin  section  where  thickness 
is  less  than  design  value. 
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TABLE  2.  Summary  of  Statistical  Results  on  Thickness 
of  Concrete  Pavement 
(Kennedy,  Hudson  and  McCullough) 


Project 

Number 

Number  of 
Samples 

Mean 

(inches) 

Variance 

(inches) 

Standard 

Deviation 

(inches) 

C V 
(*) 

8-In.  Uniform  Thickness 

1 

34 

8.66 

0.192 

0.435 

5.0 

2 

39 

8.42 

0.171 

0.415 

4.9 

3 

48 

8.35 

0.040 

0.200 

2.4 

4 

S8 

8.36 

0.077 

0.276 

3.3 

5 

61 

8.05 

0.035 

0.185 

2.3 

6 

66 

8.11 

0.089 

0.300 

3.7 

7 

73 

8.06 

0.112 

0.333 

4.1 

Pooled  Values 

— 

8.29 

0.088 

0.300 

3.6 

9-In.  Uniform  Thickness 

1 

35 

9.25 

0.0466 

0.210 

2.3 

2 

SI 

9.19 

0.121 

0.350 

3.8 

3 

58 

9.28 

0.048 

0.220 

2.4 

4 

65 

9.18 

0.060 

0.240 

2.6 

S 

74 

9.20 

0.185 

0.430 

4.7 

6 

88 

9.11 

0.029 

0.170 

1.9 

Pooled  Values 

— 

9.20 

0.083 

0.290 

3.1 

10-In.  Uniform  Thickness 

1 

64 

10.38 

0.061 

0.240 

2.3 

2 

124 

10.34 

0.079 

0.280 

2.7 

3 

132 

10.35 

0.079 

0.230 

2.2 

4 

141 

10.28 

0.083 

0.270 

2.8 

Pooled  Values 

--- 

10.34 

0.069 

0.270 

2.6 

5 
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coefficients  of  variation  nay  be  found  in  Table  3.  Mean  values  ranged 
from  2.90  x  10^  psi  to  3.89  x  10^  psi  with  an  average  of  3.45  x  10^  psi. 
Coefficients  of  variation  ranged  from  21  to  49  with  an  average  of  34.4. 
Tables  4a,  4b,  and  4c  are  the  results  of  a  series  of  laboratory  tests. 

These  results  show  that  the  modulus  value  reaches  a  slight  peak  at  an 
age  of  60  days  and  then  decreases  slightly  with  age.  The  mean  modulus 
value  was  4.54  x  106  psi  with  a  range  of  4.05  x  106  psi  to  5.16  x  106  psi. 
The  standard  deviation  increased  with  time  until  the  60  day  mark  was 
reached  and  then  it  decreased  slightly.  The  standard  deviation  ranged 
from  0.06  x  10^  psi  to  0.22  x  10^  psi.  The  coefficient  of  variation 
ranged  from  1.2%  to  4.8%. 

Poisson's  Ratio 

Tables  4,  5,  and  6  are  a  summary  of  data  from  Gibeaut  (4),  who 
investigated  the  Poisson's  Ratio  of  PCC.  Values  tended  to  decrease 
with  time  whereas  the  standard  deviation  tended  to  increase.  There  was 
no  true  trend  for  the  coefficient  of  variation.  The  mean  value  of  Poisson's 
ratio  ranged  from  0.162  to  0.212  with  an  average  of  0.188.  The  standard 
deviation  ranged  from  0.036  to  0.019  and  the  coefficient  of  variation  ranged 
from  9.3%  to  20.2%. 
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TABLE  3.  Simmary  of  Test  Results  for  Cored  Specimens  from 
PCC  Airfield  Pavements  (Kennedy) 


Airport 

Project 

Identification 

Number 
of  Tests 

Indirect  Modulus 
of  Elasticity 

Mean  CV 

(106  psi)  (%) 

O'Hare 

Runway 

9R-27L 

48 

47  . 

3.39*  34 

15 

3.33*  28 

Runway 

4R-22L 

20 

2.90**  26 

39 

3.89*  49 

Palmdale 

Runway  7-25 
(Existing  Pavement) 

11 

Runway  7-25 
(Overlay) 

38 

39 

3.86*  32 

Taxiway  A 

Overlay 

7 

8 

3.39*  40 

Taxiway  B 

Overlay 

20 

3.09*  21 

Midway 

Runway  4R-22L 
13R-31L 

11 

3.21**  43 

Richmond 

Taxiway  5-4 

Runway  2 

16 

3.02**  31 

*  Assumed  Poisson's  ratio 
-  0.20 

**  Assumed  Poisson's  ratio 
•  0.1S 

Weig 

hted 

— 

Average 

3.45  34.4 

Variation 

Limits 

2.90  -  3.89  21  -  49 

Range 

0.99  28 
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Table  4 


Individual  Sonic  E  and  v  Values 
(Expanded  from  Gibeaut) 


50  Series 

Age 

Days 

E  x  106 
psi 

E  x  106 
psi 

xl06psi 

V 

p 

H 

CV 

\ 

4.71 

14 

4.88 

4.86 

4.88 

0.15 

■ 

0.202 

0.019 

9.41 

5.07 

- 

0.215 

4.72 

0.147 

28 

4.83 

4.95 

4.90 

0.16 

0.175 

0.201 

0.184 

0.029 

15.76 

5.10 

0.211 

4.74 

60 

4.84 

5.00 

4.94 

0.18 

0.173 

0.218 

0.193 

0.036 

18.65 

5.16 

0.229 

4.68 

R9 

90 

4.83 

4.89 

0.19 

0.178 

0.036 

20.22 

4.89 

WSSm 

5.14 

HI 

4.71 

R5I 

159 

4.82 

4.76 

4.84 

0.16 

M 

0.172 

0.027 

15.75 

5.08 

HI 
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Table  6 


Individual  Sonic  E  and  p  Values 
(Expanded  from  Bibeaut) 


70  Series 

Age 

Ex  106 

Ex  106 

oc  x  106 

cv 

Days 

psi 

psi 

psi 

% 

mrrm 

0.263 

14 

HQ 

4.14 

O.OS 

0.203 

0.187 

0.212 

0.035 

16.51 

m 

0.193 

Kl 

0.235 

28 

il| 

4.16 

0.07 

0.190 

0.187 

0.206 

0.022 

10.68 

mm 

0.212 

■a 

0.240 

60 

mm 

4.17 

0.08 

0.193 

0.187 

0.212 

0.026 

12.26 

4.20 

0.226 

mm 

0.215 

90 

wm 

4.15 

0.08 

0.191 

0.172 

0.198 

0.020 

10.10 

4.17 

.  _ 

0.212 

. 

1 

m 

130 

4.13 

0.06 

89 

0.195 

0.025 

12.82 

4. IS 

0.211 
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Modulus  of  Rupture 

Current  rigid  pavement  design  methods  utilize  the  coefficient  of 
variation  in  determining  the  design  modulus  of  rupture.  Values  of  C.V. 
used  are: 

<  10%  for  excellent  construction  control 

10-1S%  for  good  control 
15-20%  for  fair  control 
>  20%  for  poor  control 

Tables  7  and  8  show  that  in  most  cases  excellent  construction 
control  may  be  maintained. 

ASPHALT  CONCRETE  MATERIALS 
Thickness 

As  for  the  thickness  of  PCC  pavements,  there  is  little  information 
on  the  wide  range  of  pavement  thicknesses.  Yoder  and  Witczak 
suggest  a  typical  standard  deviation  range  for  new  highway  construction 

0.3  to  0.8  inches  for  asphalt  concrete  pavements. 

Dynamic  Modulus 

Table  9  is  a  summary  of  laboratory  tests  which  were  conducted  to 
determine  the  dynamic  modulus  of  asphalt  concrete.  It  may  be  observed 
that  the  coefficient  of  variation  within  a  project  is  a  function  of 
temperature  and  that  there  is  no  specific  trend  of  frequency.  It  should 
also  be  noted  that  the  in-situ  (field)  variability  should  be  greater 
than  the  laboratory  variability.  In  general,  ranges  of  the  CV  values  as 
a  function  of  temperature  are:  (9%-16%)  (11 19%)  and  (20%-23%)  for  tem¬ 
peratures  of  40°F,  70#F  and  100*F,  respectively. 
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Table  7 


Average  7-Day  Flexural  Strength  and  Standard  Deviations, 
Portland  Cement  Concrete  Pavement  (Brown) 


No.  of 
Samples,  n 

Strength,  psi 

Standard 

Deviation 

5,  psi 

Coefficient 

of 

Variation  (%) 

414 

688 

73 

10.6 

74 

552 

SO 

9.1 

324 

713 

77 

10.8 

76 

642 

75 

11.7 

8 

679 

30 

4.4 

42 

645 

18 

2.8 

44 

705 

64 

9.1 

170 

736 

55 

7.5 

38 

680 

120 

17.6 

16 

591 

52 

8.8 

41 

676 

92 

13.6 

8 

640 

55 

8.6 

1255 

Avg.  662 

63 

9.5 
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Table  8 


Average  28-Day  Flexural  Strength  and  Standard 
Deviations,  Portland  Cement  Concrete 
Pavement  ( Brown 1 


No.  of 

Samples,  n 

Avg  Flexural 
Strength,  psi 

Standard 

Deviation 

6,  psi 

of  Variation (%) 

582 

781 

61 

7.8 

146 

719 

57 

5.1 

312 

862 

83 

9.6 

101 

753 

70 

9.3 

82 

774 

34 

4.4 

26 

734 

26 

3.5 

735 

739 

66 

8.9 

67 

828 

122 

4.7 

16 

688 

40 

5.8 

82 

840 

68 

8.1 

8 

717 

60 

8.4 

2157 

Avg  766 

61 

8.0 
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Flexural  Stiffness 


Table  9  gives  a  summary  of  laboratory  variability  of  the  flexural 
stiffness  of  asphalt  concrete.  For  a  test  temperature  of  68#F,  the 
mean  stiffness  value  was  326.0  x  103  psi  with  a  range  of  190.0  x  103  psi 
to  6S6.0  x  103  psi.  The  standard  deviation  ranged  from  4.0  x  103  psi 
to  92.0  x  103  psi  with  an  average  of  49  x  103  psi  while  the  coefficient 
of  variation  ranged  from  2.2%  to  23.8%  with  an  average  of  13.1%. 

Modulus  values  of  asphalt  concrete  used  in  base  materials  seem  to 
have  more  variability  than  the  modulus  values  of  asphalt  concrete  used 
in  surface  courses.  Test  results  from  blackbase  highway  projects  (Table  10) 
show  that  the  mean  modulus  value  ranged  from  35  x  103  psi  to  91.5  x  103  psi 
with  a  weighted  average  of  58.8  x  10  psi.  The  coefficients  of  variation 
ranged  from  25%  to  52%  with  a  weighted  average  of  40%.  There  is  no  refer¬ 
ence  made  to  the  temperature  used  for  these  tests.  Therefore,  since  the 
modulus  of  elasticity  is  a  function  of  temperature,  further  research 
should  be  conducted  before  using  the  values  obtained  from  this  report. 

Table  11  is  a  summary  of  flexural  stiffness  measurements  on  field 
samples.  These  values  were  found  by  using  the  pulse  load  method.  For  a 
test  temperature  of  68°F,  mean  values  of  flexural  stiffness  ranged  from 
1.34  x  103  psi  to  1.79  x  103  psi  with  an  average  of  1.58  x  103  psi.  The 
coefficient  of  variation  ranged  from  23.5%  to  27.6%  with  an  average  of 
25.0%.  For  a  test  temperature  of  40®F,  the  flexural  stiffness  ranged 
from  5.90  x  103  psi  to  7.12  x  103  psi  with  an  average  of  6.71  x  103  psi. 

The  coefficient  of  variation  ranged  from  18.8%  to  27,2%  with  an  average 
of  22.4%. 
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Poisson’s  Ratio 


Kennedy,  Hudson,  and  McCullough  (£)  investigated  the  variability 
of  the  Poisson’s  ratio  of  asphalt  concrete.  Using  the  results  from  IS 
specimens,  they  found  that  the  mean  value  of  was  0.40  and  that  its 
coefficient  of  variation  was  27%.  In  Table  10,  the  mean  value  of 
Poisson's  ratio  varied  from  0.16  to  0.34  with  an  average  of  0.25.  The 
coefficient  of  variation  ranged  from  38%  to  75%  with  a  weighted  average 
of  52%.  There  is  no  reference  made  as  to  what  temperature  that  was 
measured  at.  Therefore,  since  is  a  function  of  temperature,  further 
research  should  be  conducted  before  using  these  values. 

CEMENT  TREATED  MATERIALS 
Thickness 

For  cement-treated  bases  with  a  thickness  of  4  to  8  inches,  the 
standard  deviation  ranged  from  0.60  to  0.72  inches.  The  coefficient  of 
variation  ranged  from  7.5%  to  18%  (9). 

Modulus  of  Elasticity 

There  are  vary  large  variations  in  the  modulus  of  elasticity  for 
cement-treated  bases.  Table  12  shows  that  mean  values  of  the  modulus 
ranged  from  0.6  x  106  psi  to  1.90  x  10^  psi.  The  standard  deviation 
ranged  from  0.36  to  106  psi  to  1.19  x  106  psi.  The  coefficient  of 
variation  ranged  from  53%  to  83%. 

SUBCRADE 

Modulus  of  Subgrade  Reaction 

Treybig,  Hudson,  and  McCullough  (8)  analyzed  data  from  the  AASHO 
Road  Test.  This  subgrade  material  had  an  average  modulus  of  subgrade 
reaction  of  100  pci  and  a  coefficient  of  variation  of  16%.  This 
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TABLE  10.  Summary  of  Test  Results  for  Cored  Specimens  from 
Blackbase  Highway  Projects  in  Texas  (Kennedy) 


Project 

Identification 


Number  of 
Specimens 


Distance 

Covered 

(miles) 


Indirect 
Modulus  of 
Elasticity 

Mean  CV 

(103  psi)  {%) 


Weighted  Average 


CV  of  means  (%) 


Variation 


Indirect 

Poisson's 

Ratio 

CV 

Mean  (%) 


38.6 

32 

0.34 

39 

91.5 

29 

0.28 

40 

44.9 

46 

0.16 

58 

87.3 

62 

0.16 

73 

35.0 

40 

0.26 

57 

86.1 

59 

0.23 

47 

55.2 

44 

0.24 

41 

42.2 

24 

0.20 

64 

55.2 

33 

38 

64.7 

34 

67 

58.8 

40 

0.25 

52 

36 

.28 

-- 

35-91.5 

24-62 

.16-. 34 

38-73 

TABLE  11.  Flexural  Stiffness  Measurements  on  Field  Samples 
of  Asphaltic  Concrete  Using  Pulse  Loading  Method 
(Kennedy) 


Measured 

Stiffness  (psi  x  10^) 

68°F 

40°F 

No.  of 

Sample 

Speci- 

Std. 

CV 

Std. 

CV 

Group 

mens 

Mean 

Dev. 

(%) 

Mean 

Dev. 

(%) 

Specimens 

from  Surface 

Course 

1 

19 

1.79 

0.42 

23.5 

6.80 

1.53 

22.5 

2 

20 

1.65 

0.39 

23.6 

7.03 

1.91 

27.2 

3 

20 

1.52 

0.41 

26.9 

7.12 

1.41 

19.8 

4 

19 

1.34 

0.37 

27.6 

5.90 

1.11 

18.8 

Lab 

25.0 

22.4 

Compacted 

26 

1.29 

0.22 

17.0 

5.76 

0.74 

12.8 

Specimens 

from  Base  Course 

1 

12 

1.57 

0.42 

26.6 

5.95 

1.54 

25.° 

2 

12 

1.39 

0.26 

18.7 

5.66 

3.14 

55.5 

3 

8 

1.47 

0.41 

27.9 

4.40 

0.90 

20.4 

4 

10 

1.42 

0.42 

29.6 

4.96 

1.22 

24.6 

25.3 

33.0 

Lab 

Compacted 

29 

1.19 

0.19 

16.0 

5.31 

1.50 

28.2 
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HIGHWAYS  AIRFIELDS 


TABLE  12.  Summary  of  Properties  for  Cement -Treated  Bases  for 
Airfield  and  Highway  Pavements  (Kennedy) 


Distance 


Indirect 
Modulus  of 
Elasticity* 


Identification 

Type  of 
Material 

Number  of 
Specimens 

Covered 

(miles) 

Mean 

(10°  psi) 

CV 

(%) 

» 

C  CM 

3  CM 

PC  1 

Top 

- 

IS 

1.9 

53 

0)  TT 

fr¬ 
ee  X 

Bottom 

15 

1.8 

66 

o  * 

Combined 

30 

1.8 

59 

DFW 

King burg 
Gravel 

4 

- 

- 

Weighted  Average 

12-A 

Sand 

Shell 

32 

1.76 

19-A 

Soil 

Cement 

20 

1.4 

1.05 

83 

19- B 

Soil 

Cement 

19 

1.2 

0.73 

57 

20- A 

Burned 

Clay 

29 

1.5 

0.60 

60 

Weighted  Average 

0.09 

68 

*Assumed  Poisson's  Ratio-0.22 

CV  of  Means (%) 

50 

- 

Variation 

Limits 

0.60-1.76 

Range 

1.16 
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coefficient  of  variation  should  only  be  applicable  to  a  k  of  100  pci. 
Further  information  should  be  gathered  before  drawing  any  conclusions 
about  the  coefficients  of  variation  for  other  values  of  k. 

Resilient  Modulus 

The  resilient  modulus  tended  to  increase  with  a  decrease  in  deviator 

3  3 

stress.  Table  13  shows  mean  modulus  values  of  12  x  10  psi,  16  x  10  psi, 

3 

and  19  x  10  psi  for  deviator  stresses  of  8  psi,  5  psi,  and  2  psi,  re¬ 
spectively.  The  coefficient  of  variation  ranged  from  5.3%  to  52%.  The 
average  standard  deviation  was  3.75  x  10^  psi. 

California  Bearing  Ratio  (CBR) 

Table  14  shows  variability  data  for  the  subgrade  CBR.  Mean  CBR 
values  ranged  from  4.2  to  26.3.  The  standard  deviation  ranged  from  0.9 
to  8.4  while  the  coefficient  of  variation  ranged  from  17.9%  to  36.9%. 
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TABLE  13.  Summary  of  Rfesilient  Moduli  for  Subgrade  Soils  for  Airports  (Kennedy) 


Table  14 


Subgrade  CBR  Variability 
(reduced  from  Yoder  and  Witczak) 


Mean  CBR  S  CV(%) 


7.1 

1.6 

22.3 

4.2 

0.9 

21.4 

26.3 

8.4 

31.9 

20.3 

7.6 

36.9 

18.2 

4.8 

26.2 

7.8 

1.4 

17.9 

Remarks 


(In-situ;  compacted  subgrade) 
(Estimated;  after  moisture  equil) 
(In-situ;  compacted  subgrade) 
(Estimated;  after  moisture  equil) 
(In-situ;  compacted  subgrade) 
(Undisturbed  samples;  soaked) 
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Chapter  3 


SUMMARY  AND  RECOMMENDATIONS 

The  purpose  of  this  report  was  to  give  a  general  summary  of  values, 
or  a  range  of  values,  that  describe  material  variabilities.  It  summarizes 
a  state-of-the-art  review  in  variabilities  of  material  properties  which 
are  used  in  pavement  design. 

If  a  probabilistic  design  methodology  has  been  developed,  a  measure 
of  the  variability  of  the  specific  parameter  in  question  must  be  determined. 
This  variability  information,  however,  is  not  always  readily  available. 

This  report,  then,  is  to  be  used  as  a  guide  in  determining  the  level  of 
variability  to  be  used  for  a  design  parameter.  Table  15  is  a  table  of 
variability  recommendations  which  may  be  used  for  different  levels  of 
inherent  variability.  As  variability  data  becomes  available,  it  should 
be  used  to  upgrade  or  to  check  Table  15.  If  variability  information 
pertaining  to  a  particular  project  is  available,  it  should  be  used  in  lieu 
of  the  information  from  Table  15. 
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Table  15.  Variability  Recommendations 


Coefficient  of  Variation  (%) 


Inherent 

Low 

Variability 

Average 

High 

PCC  MATERIALS 

Thickness 

1-3 

4-6 

7-9+ 

Modulus  of  Elasticity 

20-30 

30-40 

40-50+ 

Poisson's  Ratio 

8-12 

13-16 

17-20+ 

Modulus  of  Rupture 

10-13 

14-17 

18-20+ 

ASPHALT  CONCRETE  MATERIALS 

Thickness 

1-5 

5-10 

10-15 

Dynamic  Modulus 

Temp:  40°F 

8-10 

11-13 

14-16+ 

70°F 

10-12 

13-15 

16-19+ 

100°F 

18-20 

21-22 

23-24+ 

Base  Modulus  of  Elasticity 

25-35 

35-45 

45-55+ 

Flexural  Stiffness 

Temp:  40°F 

15-20 

20-25 

25-28+ 

68°F 

20-23 

24-26 

27-30+ 

Poisson's  Ratio* 

35-48 

49-62 

63-75+ 

CEMENT- TREATED  MATERIALS 

Thickness 

6-10 

11-15 

16-19+ 

Modulus  of  Elasticity 

53-63 

63-73 

73-83 

SUBGRADE  MATERIALS 

Modulus  of  Subgrade  Reaction* 

10-20 

20-35 

35-50+ 

Resilient  Modulus 

10-20 

20-35 

35-50+ 

CBR 

15-22 

23-31 

32-40+ 

* 

See  Text  for  additional  information 
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Chapter  1 
INTRODUCTION 

For  the  past  several  years,  probabilistic/reliability  analyses 
have  slowly  been  introduced  in  pavement  design  schemes  to  account  for 
the  known  variation  of  the  salient  design  parameters.  In  existing 
deterministic  design  methods,  material  variability  effects  on  pavement 
performance  may  be  only  indirectly  taken  into  account.  Frequently,  the 
design  parameters  (unique  values)  are  determined  from  either  laboratory 
or  field  results.  If  multiple  values  are  determined,  the  engineer  may 
select  (from  the  distribution  of  values)  a  design  percentile  level 
(e.g.  .  .  .  85th)  for  use  in  the  deterministic  solution.  The  confidence 
level  adopted  and  reflected  by  the  percentile  value  is  currently  based 
on  the  engineering  judgment  and  experience  of  the  designer.  This  approach, 
however,  cannot  be  easily  or  directly  used  to  quantitatively  evaluate 
the  effect  of  material  variability  in  design  analysis,  cost  studies  or 
reliability  solutions.  The  probabilistic  approach  which  directly  imple¬ 
ments  the  variability  distribution  of  all  significant  parameters  provides 
the  only  true,  accurate  measure  of  the  reliability  of  the  entire  pavement 
system. 

The  probabilistic  analysis  of  pavement  performance  is  based  on 
concepts  and  formulations  developed  in  the  statistical  field  dealing  with 
system  reliability.  This  volume  presents  a  summary  of 

(a)  the  definitions  and  mathematical  formulation  of  the  general 
problem  applicable  for  any  probability  density  distribution; 

(b)  the  existing  approaches  for  developing  pavement  design  solutions 
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in  probability  and  reliability  terms.  It  includes:  (i) 
the  interference  theory  or  stress -strength  approach,  (ii) 
the  simulation  technique  and  (iii)  the  approximate  closed-form 
probabilistic  formulation,  based  on  the  Taylor  series  expansion 
of  the  dependent  variable. 

The  various  theories  are  discussed  from  the  viewpoint  of  their 
possible  implementation  in  the  rigid  pavement  design  procedure  of  the 
USACE  for  airfield  systems. 
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Chapter  2 


DEFINITIONS  AND  MATHEMATICAL  FORMULATION 

General 

This  chapter  presents  a  brief  overview  of  basic  definitions  and 
theorems  related  to  the  probabilistic  formulation  of  design  variables 
applied  to  the  pavement  design  solution.  Further  background,  details 
proof  of  theorem,  etc.  .  .  .  can  be  found  in  most  standard  textbooks 
on  Probability  and  Reliability  (see,  for  example,  Hines  and  Montgomery  (1)) 


Probability  and  Reliability  Concepts 

Reliability  is  defined  as  the  probability  of  success  of  an  event 
and/or  a  statement  of  the  error  or  precision  of  an  estimate.  This  event, 
which  can  take  a  finite  or  infinite  number  of  values,  is  called  a  random 
variable.  The  distribution  of  these  values  is  related  to  the  probability 
distribution  of  the  random  variable  by  the  concept  of  probability,  i.e. 


P(X  =  xk)  =  f(xk) 


V  =  1,  2,  .  .  . 


(2.1) 


with  f^)  dxk  =  1 

where  P(X  *  j^)  represents  the  probability  that  the  random  variable 
takes  the  value  xk,  and  f(xk)  is  the  density  probability  distribution 
of  X  at  Xj^.  This  can  also  be  expressed  through  the  cumulative  distribution 
function  defined  as 


P(Xlxk)  *  F(xk)  -rf(u) 


du 


(2.2) 
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When  success  is  defined  by  X  taking  all  values  greater  than  x^, 
reliability  takes  the  form  of 


R  =  P(X  >  xfc)  =  1  -  F(xk)  =  1  -  /  f(u)  du 


(2.3) 


It  can  be  seen  that  reliability  can  be  directly  computed  from  the 
probability  distribution  of  the  random  variable,  directly  from  its 
cumulative  probability  distribution,  or  from  integration  of  its  density 
probability  distribution. 


Extension  to  Multi  Variate  Expressions 


The  case  of  one  random  variable  is  the  simplest  solution.  However, 
in  practice,  two  or  more  random  variables  are  involved  in  the  process. 

If  one  examines  the  case  of  two  random  variables,  the  joint  distribution 
function  of  the  random  variables  X  and  Y  is  defined  as: 


P(x  -  Xy  Y  *  yk)  =  f( Xy  yk)  (2.4) 

with  f(x..,  yk)  dx_.  dyR  =  1 

The  above  probability  expresses  the  occurence  of  two  events  which 
can  be  either  dependent  or  independent.  In  the  case  of  independent 
variables,  the  "product  rule"  applies,  i.e. 

P(X  -  x.,  Y  ■  yk)  -  [P(X  =  x.))  •  [  P(Y  =  yk)]  (2.5) 
or  f  (x,  y)  *  [fj(x)3  *[  f2  (y)l 

where  fj(x)  and  f^(y)  are  called  the  marginal  probabilities  of  X  and  Y. 
The  extension  of  the  one  random  variable  to  two  or  more  random  variables 
is  easily  done: 
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(2.6) 
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Change  of  Variables 


If  X  is  a  continuous  random  variable  with  a  probability  density 
function  f(x),  U  *  #  (X)  and  X  *  f(U);  then  the  probability  density 
of  U  is  given  by  g(u)  where 

g(u)  *f  j  =  f(V(u))|  *’(u)|  (2.11) 

If  X  and  Y  are  continuous  random  variables  having  a  joint  density 
probability  function  f(x,y),  U  “^(X.Y),  V  =  *2(X,Y),  X  =*j(U,V)  and 
Y  *  #2(U,V);  then  the  joint  density  function  of  U  and  V  is  given  by 
g(u,v)  where 

g(u,v)  =  f(x,y)  =  ft^jCu.v),  *2(u,v))‘  | J | 

where 

J  *  the  Jacobian  of  the  transformation  (2.12) 


Probability  Distributions  of  Functions  of  Random  Variables 

If  X  and  Y  are  continuous  random  variables  and  U  *  ♦ j(X,Y),  then 
the  density  function  of  U  is  the  marginal  density  obtained  from  the 
joint  density  of  U  and  V  (while  V  is  chosen  arbitrarily,  V  *  X  or  V  *  Y). 

If  f(x,y)  is  the  joint  density  for  X  and  Y,  then  the  density 
function  g(u)  of  the  random  variable  U  »#j(X,Y)  is  found  by  differ¬ 
entiating  with  respect  to  u  the  distribution  function  given  by 

G(u)  *  Pl*j(x#y)£u]  ■  /  /R  f(x,y)  dx  dy  (2.13) 

where  R  is  the  region  for  which#  j(x,y)£u 
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The  density  function  of  the  sun  of  two  continuous  random  variables 
X  and  Y,  i.e.  of  U  *  X+Y,  having  joint  density  function  f(x,y)  is 
given  by: 

g(u)  *  f.m  f(x,  u-x)  dx  (2.14) 

For  the  special  case  where  X  and  Y  are  independent, 


f(x,y)  *  t  f^x)]  ‘[f2(y)] 

4*oo 

and  g(u)  »  /  [fjtx)]* [f2(u-x)  dx]  =  fj  *  f2 

This  is  called  the  convolution  of  fj  and  f2# 


(2.15) 


Expectation 

For  a  continuous  random  variable  X  having  a  density  function  f(x), 
the  expectation  of  X  is  defined  as 

*  CD 

E[X]  *  /  xf(x)  dx  (2.16) 

-  00 

The  expection  of  X  is  very  often  called  the  mean  of  X  (denoted 
by  u  x)  and  is  a  measure  of  central  tendency. 

If  X  is  a  random  variable  then  Y  =  g(X),  a  function  of  X  is  also 
a  random  variable.  The  expectation  of  Y  is 

E[Y]  -  E[g(X)J  -  /* "  g(x)  *  f(x)  dx  (2.17) 

oo 

For  a  function  of  two  random  variables  g(X,Y) 

*  m  ♦  eo 

E[g(X,Y)]  •  /  /  g(x,y)  •  f(x.y)  dx  dy  (2.18) 

*  •  -  00 

The  variance  of  a  random  variable  is  defined  by: 


VAR[X] 


E[(x  -u)2] 


(2.19) 


2 

where  g(X)  =  (x- u)  is  a  function  of  X.  According  to  the  above: 

Efe(X)]  =  E((X-y)2]  =  X-m  (x-u)2f(x)  dx  (2.20) 

The  variance  (or  its  square  root,  the  standard  deviation)  is  a  measure 
of  the  dispersion  or  scatter  of  the  values  of  the  random  variable. 

The  covariance  of  two  random  variables  is  defined  similarly  as: 

Cov  [  X,Y ]  =  E[(X-vx)  (Y-v  )] 

*  /  "  /  (x-y  )(y-P  )  f(x,y)  dx  dy  (2.21) 

-  OD  -00  * 

which  is  equal  to  zero  when  the  random  variables  are  independent.  When 
X  and  Y  are  dependent,  the  following  relation  between  the  covariances, 
the  variances  and  the  correlation  coefficient  can  be  found  by: 


Cov  [  X,Y  ] 


xy 


Var  [X  ] 
o 

y 


Var  [  Y  ] 


where  p  is  the  correlation  coefficient. 


(2.22) 


Higher  orders  of  expectation  of  the  random  variable  are  defined 
in  the  literature,  to  describe  mathematically  the  probability  density 
distributions  of  the  random  variable  (skewness,  kurtosis,  etc.) 
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Chapter  3 


PROBABILISTIC  THEORIES  FOR  PAVEMENT  DESIGN 
General  Background 

The  reliability  of  a  system  composed  of  several  units  depends 
upon  the  reliability  of  these  units  and  upon  the  way  they  are  intercon¬ 
nected  (in  series  or  in  parallel  as  defined  in  the  reliability  theory). 

The  reliability  of  each  unit  can  be  determined  by  experiment  or  assessed 
from  the  stress- strength  method. 

At  present,  the  rigid  pavement  design  framework  is  based  on  slab 
failure  only,  without  any  consideration  of  subgrade  or  joint  failure. 

It  is  therefore  a  one  unit  system  with  a  large  number  of  variables 
(traffic,  material  characteristics  and  slab  geometry).  A  basic  approach 
to  the  problem  would  be  that  of  the  stress- strength  analytical  methodology. 
The  second  one  would  be  based  on  simulation,  to  bypass  the  complexity 
of  the  basic  approach,  and  the  third  one  would  be  an  approximate  one  based 
on  normal  distribution  assumption  and  Taylor  series  expansion  of  the 
functions.  These  approaches  are  discussed  in  detail. 

Stress- Strength  Approach  (Interference  Theory) 

The  unit  is  assumed  to  fail  where  the  "stress" (Y)  reaches  its 
ultimate  value,  its  "strength"  (X).  The  probability  density  distribution 
of  both  the  stress  and  the  strength  must  be  known.  If  it  is  assumed 
that  the  stress  and  the  strength  are  independent,  the  probability  of 
failure  is  expressed  as  the  probability  of  X  being  lower  or  equal  to  a 
given  value  of  Y,  i.e., 

y 

P  (X<y/  Y-y)  »  /  f.(x)  dx  (3.1) 

at  ^ 
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Multiplying  by  the  probability  that  Y  is  in  the  neighborhood  of  y, 
one  obtains  the  joint  probability  function 

y 

P  (x<y,  y<Y<y  ♦  dy)  =  f  fj(x)  f2(y)  dx  (3.2) 

and  extending  this  to  the  whole  range  of  Y, 

♦do  y 

P  (X1Y)  =  /  /  f.(x)  •  f,  (y)  dx  dy  (3.3) 

The  reliability  of  the  unit  is  given  by 

+oo  y 

R  *  1  -  P(X£Y)  =  1  -  /  /  f  (x)  f,(y)  dx  dy 

“00  -oo  ‘ 

+  00 

*  1  -  /  F(y)  f,(y)  dy  (3.4) 

where  F(y)  is  the  cumulative  probability  function  of  x  evaluated  at  y. 

It  is  seen  that  the  reliability  can  be  computed  if  the  density 
probability  functions  of  the  stress  and  the  strength  are  known.  However, 
only  for  very  simple  cases  can  the  above  integral  be  evaluated  in  a 
closed  form  model.  Yet  the  stress  itself  can  be  a  function  of  random 
variables,  and  its  density  probability  function  can  be  expressed  by  the 
above  theorem  by  differentiating 

G(u)  -  //R  f(x,y)  dx  dy  (3.5) 

with  respect  to  u. 

The  problem  involving  a  number  of  variables  is  found  to  be  very 
complex,  first  to  define  the  probability  density  distribution  of  the 
'’stress’'  and  of  the  "strength”,  and  second,  to  evaluate  the  joint 


probability  (or  the  probability  of  failure).  It  should  be  noted  that 
for  the  case  of  the  normal  density  distribution  of  the  stress  and  the 


strength,  the  probability  of  failure  P(X_Y)  is  also  normally  distributed  (2). 

In  passing,  it  should  also  be  noted  that  log-normal  distributions  or  other 
statistical  distributions  (e.g.  Beta)  can  be  utilized  in  the  above  noted  procedure. 

Simulation  Technique 

The  technique  of  simulation  has  become  widely  used  for  estimating  dis¬ 
tribution  parameters  of  the  variable  and  for  providing  insight  into  the  cause 
and  effect  relationship  within  a  system.  In  simulation,  the  system  is  divided 
into  elements  whose  behavior  can  be  predicted,  at  least  in  terms  of  probability 
distributions.  They  are  then  combined  in  their  natural  order,  allowing  the 
computer  to  present  the  effect  of  their  interaction  on  each  other.  After 
constructing  the  model,  it  is  activated  by  generating  appropriate  input  data 
to  simulate  the  actual  aggregate  behavior  of  the  system  over  time. 

The  technique  of  simulation  has  not  been  applied  widely  in  estimating 
pavement  performance.  Few  studies  (e.g.,  Ullitz  (3)have  advanced  reliability 
based  pavement  performance  schemes.  Simulation  can  be  applied  for  the  rigid 
pavement  case  to  estimate  the  probability  distribution  of  the  overall  behav¬ 
ior  (performance)  of  the  pavement,  given  the  probability  distributions  of 
the  individual  variables  and  to  estimate  the  effect  of  each  variable  on  the 
overall  (aggregate)  behavior.  Formulation  of  the  model  is  quite  simple 
and  requires  the  following  steps: 

(1)  Define  the  random  variables  and  their  probability  distributions; 

(2)  Express  the  rules  (equations)  that  link  the  different  variables 
to  the  performance  variable  -  the  one  whose  average  value 

and  probability  distribution  is  of  interest. 

Activation  of  the  model  is  performed  by  generating  input  of  the 
random  variables  from  their  probability  distribution  and  deriving  the 
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performance  variable.  Repeating  the  process  a  sufficient  number  of  times 
gives  the  description  of  the  probability  distribution  of  the  output 
variable.  The  appropriate  flow  chaTt  for  simulating  rigid  pavement 
performance  is  shown  in  Fig.  3.1. 

For  the  case  of  several  variables,  and  in  order  to  estimate  the 
cause  and  effect  relationships,  a  factorial  simulation  experiment  must 
be  conducted.  The  factorial  experiment  must  be  designed  in  such  a  way 
to  allow  to  extract  the  individual  effects  of  the  random  variables. 

It  should  be  noted  the  number  of  solutions  required  (N  -  in  Fig.  3.1) 
may  be  quite  large  in  OTder  to  accurately  reproduce  probability  distri¬ 
butions  of  the  various  variables.  Because  of  this,  the  simulation 
technique  should  be  restricted  to  situations  that  are  too  complex  to 
be  handled  by  closed  form  formulation  and/or  for  cases  where  there  is 
strong  evidence  that  probability  distribution  of  variables  are  far  from 
the  normal  approximation. 

Approximate  Closed  Form  Probabilistic  Formulation 

According  to  the  probability  concepts  briefly  summarized  above, 
the  expected  value  (mean)  of  a  function  of  variables  is  expressed  as: 

E  [g(X.)]  «  /  g(xi)  •  f(x1,...,x_)  dxr..dxn  (3.6) 

*00 

where  f(Xj,  . ...x^)  is  the  joint  probability  distribution  of  all  Tandom 
variables  X..  This  is  ’•ather  difficult  to  evaluate.  Instead  of  evalu- 

l 

ating  integrals,  the  function  may  be  expressed  in  terms  of  the  Taylor 
series  expansion  evaluated  near  the  variable  means  and  its  expected 
value  derived: 
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Fig.  3.1  SIMULATION  SCHEME 
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where  R  is  the  residual. 

In  the  case  of  the  first-order  Taylor  series  approximation  (corres¬ 
ponding  to  linear  approximation),  the  following  can  be  found. 


EtgCXj)] 


Var 


g(xlf  x2, 


...  v 


I.CV1 


Var  [X.]  ♦ 


(3.8) 


n  n 

i=l  5 

ki»i 

where: 

g(X1,X2f...Xn)  is  the  function  g(X^)  evaluated  at  *  mean 
values  of  the  variables  X. 


3  gCXj) 

T~*r~1xi 

VartX.] 

Cov  [Xv  xg 


is  the  derivative  of  g(X^)  with  respect  to 
X^,  evaluated  at  X..  *  X^ 

is  the  variance  of  variable  X^ 

is  the  covariance  of  variables  Xj  and  X^. 
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In  the  case  of  two  variables,  it  can  be  shown  that: 


£[g(x,y)]  -  g(x,  y) 

2 

Var(g(x,y)]  =  (|f| -)2  Var[xJ  +  (f£|  J  Var[y]  + 

ax'  x  ayi  y  (3.9) 

2  <HI  =•  5  <t?U  ’  c°v  <’■* 

*  y 

It  should  be  noted  that  the  first  order  approximations  of  the  Taylor 
series  is  in  most  engineering  cases  a  good  approximation.  Taking  the 
second-order  approximation  leads  to  the  following  mean: 


EtgfX^]  *  g(Xj,  X2>. 


n 


n 


1 

2 


32g(X^ 


n  a2g(x.) 

L 

m 


Var[X^]  + 


Cov[Xi> 


xkl  -  R 


(3.10) 


Since  VarfX.^  and  CovfXj,  are  relatively  small,  the  approximation 
is  a  good  one  where  the  second  order  derivatives  of  g(X^)  evaluated  at 
are  also  relatively  small  so  that  their  product  will  be  negligible.  The 
variance  term  of  the  second- order  approximation  involves  higher  order  of 
measures  of  central  tendency  (the  skewness  and  kurtosis)  and  multiplication 
of  the  variances  and  covariances  of  the  variables.  As  such,  the  precise 


derivation  becomes  rather  complex.  The  linear  approximation  is  improved 
in  symmetric  probability  distributions  when  the  skewness  coefficient  is 
nil.  Therefore,  the  approximate  closed  form  formulation  should  be  adopted 
only  when  the  output  variable  has  a  symmetric  density  distribution.  For 
the  closed-form  approximation  to  be  appropriate,  it  should  satisfy  the 
following  conditions:  (a)  the  output  variable  is  normally  distributed 
and  (b)  the  second  order  terms  of  the  Taylor  series  expansion  can  be  ne¬ 
glected,  in  comparison  of  the  first  order  form. 
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Chapter  4 


SUMMARY 

Three  basic  approaches  for  expressing  rigid  pavement  design  in 
probabilistic  and  reliability  terms  were  summarized.  From  the  above, 
it  is  noted  that 

(a)  in  the  stress-strength  approach,  definition  and  characteriza¬ 
tion  of  the  stress  component  need  to  be  clarified.  In  the  current  design 
method,  initial  failure  is  defined  as  the  "point  at  which  the  crack, 
which  originates  at  the  bottom  of  the  slab  and  migrates  upward  to  the 
surface,  commences  to  spall  and  ravel,  which  produces  debris  on  the 
surface  (4).  The  "stress"  that  would  cause  such  a  deterioration 
(cracking  initiation  and  propagation,  spalling  and  ravelling)  cannot 
at  this  stage  be  given  a  mechanistic  interpretation. 

The  stress-strength  approach  assumes  that  the  variables  are  inde¬ 
pendent.  Moreover,  it  is  based  on  comparing  the  values  of  the  random 
variables,  X  -  the  stress,  and  Y  -  the  strength,  which  is  equivalent  to 
defining  the  new  variable  Z  =  X-Y.  The  approach  currently  followed 
in  the  rigid  pavement  design  is  quite  different.  The  failure  is  defined 
through  the  Design  Factor  which  is  the  ratio  of  the  strength  (the  modulus 
of  rupture)  and  the  stress  (further  reduced  by  the  load  transfer  coeffi¬ 
cient).  For  this  case,  the  probability  distribution  of  the  design  factor 
remains  unknown  and  quite  difficult  to  determine. 

While  the  discussion  has  denoted  the  random  variables  in  authentic 
fasion,  it  should  be  understood  that  other  distributions  (i.e.,  log  normal. 
Beta,  Gamma)  can  be  assumed  to  represent  the  PDF  of  each  variable.  Such 
assumptions  would  cause  fur4-’  °r  mathematical  problems  in  the  "Strength- 
Stress"  approach  to  the  reliability  solution. 
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(b)  Simulation  techniques  are  powerful  for  solving  complex  problems, 
which  cannot  be  solved  in  closed  form.  However,  full  factorial  experiment 
analysis  say  be  very  computer  time  consuming,  and  hence  costly.  It 
could  be  combined  with  the  approximate  closed  form  solution  which  could 
be  used  for  evaluating  the  effect  .of  each  variable  on  the  overall 

aggregate  performance. 

(c)  Where  analytical  expressions  are  available,  or  can  be  developed 
by  statistical  regression  predictive  techniques,  the  direct  application 
of  the  Taylor  Series  expansion  (approximate  closed-form)  to  a  probabil¬ 
istic  formulation  of  reliability  analysis  appears  warranted.  It  is 
without  question  a  very  powerful  mathematical  tool  that  has  great  potential 
in  reliability  studies  of  pavement  performance. 

Because  of  these  general  considerations,  it  is  noted  that  the 
primary  method  of  developing  the  probabilistic  solutions  to  the  rigid 
pavement  performance  analysis  of  the  USACE  for  airfield  pavements  was 
the  Approximate  Closed-Form  solution  (Taylor  Series).  This  approach 
was  used  to  formulate  both  the  probabilistic  Westergaard-USACE  procedure 
shown  in  Volume  III  as  well  as  the  Multi-Layered  Elastic  Design 
(developed  by  the  USACE)  presented  in  Volume  IV.  In  the  latter  case, 
a  simulation  scheme  is  also  presented  to  clearly  demonstrate  that  the 
two  approaches  are,  from  an  engineering  viewpoint,  equivalent.  Such 
a  conclusion  only  enhances  the  utilization  of  the  Taylor  Series  approach 
to  the  reliability  problem.  Inherent  within  this  suggested  probabilistic 
approach,  is  the  necessity  to  develop  closed-form  solutions  and  models 
by  statistical  regression  techniques  for  various  variables  important 
in  the  design/analysis  solution.  The  development  of  these  particular 
expressions  are  clearly  presented  in  the  appropriate  report  Volumes. 
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Chapter  1 

INTRODUCTION 

General 

The  design  of  rigid  pavements  (plain  and  reinforced)  for  military 
purposes  is  currently  based  upon  the  classical  Westergaard  free  edge 
stress  slab  theory.  Present  U.S.  Army  Technical  Manuals  and  U.S.  Air 
Force  Manuals  use  this  theory  as  the  basis  for  their  design  methods.  The 
Westergaard  model  consists  of  a  linear  elastic  slab  supported  by  a  dense 
liquid  (equivalent  to  independent  linear  springs).  While  Westergaard  de¬ 
veloped  stress  solutions  for  loads  placed  at  various  slab  locations 
(i.e.,  edge,  comer  and  interior),  the  military  design  approach  is  based 
upon  a  (modified)  free  edge  stress  condition.  (1) 

At  present,  the  design  method  approach  is  deterministic,  i.e.,  a  unique 
pavement  system  is  designed  for  the  set  of  input  variables  which  are  also 
unique.  The  design  corresponds  to  an  acceptable  performance  level  (first 
crack)  with  an  unknown  margin  of  safety.  The  deterministic  design  approach, 
indirectly  and  only  qualitatively,  may  account  for  the  effect  of  material 
variability  on  pavement  performance  by  judicious  selection  of  the  design 
input  values.  However,  a  quantification  of  these  effects  is  possible  and 
will  improve  the  design  procedure  by  showing  the  partial  effect  of  each 
variable. 


Study  Objective 

The  aim  of  this  project  was  to  include  the  design  parameter  variabi  lity 
in  the  current  USACE  design  procedure  for  unreinforced  concrete  pavements, 
and  to  evaluate  (quantitatively)  their  effect  on  the  final  design.  The 
final  rigid  pavement  design  is  then  expressed  in  reliability  terms.  The 
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mathematical  and  analytical  development  of  the  probabilistic  approach 
used  in  this  volume  is  presented  in  Volume  II  of  the  study. 

It  was  also  the  aim  of  this  project  to  expand  the  current  design  pro¬ 
cedure  so  that  systems  with  more  than  one  foundation  layer  could  be  evalu¬ 
ated.  This  was  accomplished  through  the  use  of  the  composite  modulus  of 
reaction  which  considers  all  layers  below  the  pavement  to  respond  as  an 
"equivalent",  "lumped"  or  "composite"  foundation.  The  layer  parameters 
of  this  composite  section  are  then  used  as  the  foundation  material  respons 
within  the  Westergaard  model. 

Report  Organization 

This  volume  has  been  subdivided  into  two  major  chapters.  They  are: 
Chapter  2:  Composite  Modulus  of  Subgrade  Reaction 
Chapter  3:  Probabilistic  Analysis  of  PCC  Airfield  Design 

Chapter  2  will  discuss  the  background  of  the  composite  modulus  of 
subgrade  reaction  as  well  as  the  development  of  the  equations  which  are 
used  within  the  reliability  based  design  method. 

Chapter  3  deals  with  the  development  of  the  equation  used  to  predict 
the  maximum  stress  computations.  These  equations,  in  turn,  are  used  to 
evaluate  the  variability  of  the  pavement  design.  The  probabilistic  solu¬ 
tion  used  in  the  Westergaard  analysis  is  based  upon  the  Approximate  Closed 
form  Probabilistic  Formulation  (Taylor  Series)  discussed  in  Volume  II. 

The  detailed  mathematical  derivations  used  to  develop  this  solution  are 
contained  in  this  volume.  A  computer  program  has  been  developed  to  solve 
this  probabilistic  methodology.  A  user's  guide  and  a  program  listing 
have  been  included  in  the  Appendices. 


Chapter  2 

COMPOSITE  MODULUS  OF  SUBGRADE  REACTION 

Introduction 

The  composite  modulus  of  subgrade  reaction  was  introduced  to  take 
into  account  the  combined  effect  of  a  layered  system  beneath  the  rigid 
pavement.  This  composite  modulus  is  a  very  useful  and  necessary  tool  in 
the  evaluation  of  layered  systems.  The  Westergaard  theory  is  based  upon 
a  slab  resting  on  a  foundation  material.  As  a  consequence,  multiple 
foundation  layers  cannot  be  solved  for  directly.  However,  the  composite 
modulus  is  a  means  of  making  the  layered  system  into  a  single  equivalent 
layer  which  may  be  used  within  the  Westergaard  model.  The  composite  modulus 
is  defined  as  an  equivalent  modulus  that  will  lead  to  the  same  response 
from  the  original  layered  system  (in  this  case:  equal  deflections  of  the 
subsystem  below  the  slab).  The  major  variables  in  determining  the  composite 
modulus  are  the  layer  thickness  and  modulus.  This  chapter  presents  the  eval¬ 
uation  of  the  composite  modulus  of  subgrade  reaction  for  the  case  of  a 
layered  system.  All  of  the  composite  modulus  relationships  used  in  this 
solution  have  been  based  upon  those  contained  in  TM  5-824-3  (AFM  88-6, 

Chap.  3),  Rigid  Pavements  for  Airfields  Other  Than  Army,  August  1979. 

Current  USACE  Method  for  Determining  the  Composite  Modulus 
of  Subgrade  Reaction  (k£  Value) 

In  the  current  USACE  method  (August  1979),  several  relationships  for 
composite  modulus  (kc)  are  presented  for  different  materials  (2J) .  Figure 
2.1  shows  the  k£  relationship  for  a  well-graded  crushed  material  and 
Figure  2.2  shows  the  relationship  for  natural  sands  and  gravels  (PI <8). 

The  kc  is  developed  as  a  function  of  the  layer  thickness,  the  modulus  of 
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subgrade  reaction  and  the  material  type.  Figure  2.3  shows  the  kc  relation¬ 
ship  for  stabilized  subbase  materials.  In  this  case,  the  kc  value  is  depen¬ 
dent  on  the  modulus  of  elasticity  and  layer  thickness  of  the  stabilized  layer 
as  well  as  the  modulus  of  subgrade  reaction. 

Figures  2.1  and  2.2  are  based  on  field  plate-loading  tests  while 
Figure  2.3  is  based  on  computations  using  the  elastic  layered  theory  with 
corrections  (_3 ) .  These  k£  values  are  based  on  an  equivalent  deflection  cri¬ 
terion.  The  deflection  of  the  composite  section  is  equal  to  the  deflection 
of  the  layered  system  under  the  same  30  inch  diameter  plate. 

Formulation  of  Composite  Modulus  Equations 

The  approximate  closed-form  probabilistic  approach  presented  in  Volume  II 
and  implemented  in  this  volume  requires  the  variables  to  be  expressed  in  equation 
form.  Therefore,  the  results  of  the  composite  modulus  from  figures  2.1  to 
2.3  were  analyzed  to  derive  regression  equations.  This  was  accomplished  by 
the  use  of  the  Statistical  Package  for  the  Social  Sciences  (SPSS).  A  multiple 
regression  analysis  was  performed  on  the  information  taken  from  these  figures 
to  generate  the  equations.  (The  SPSS  outputs  may  be  found  in  Appendix  I). 

The  final  equations  selected  resulted  in  excellent  correlation  coefficients 
and  agreement  between  actual  and  predicted  values.  The  SPSS  generated 
equations  are: 

For  well-graded  crushed  materials: 

log  k  =-1.251182  +  2.219732  log  k  -  0.2949522  (log  k  )2 
c  sg  sg 

♦  0.08901252h  -  0.0004425194  h2  (2  ^ 

-  0.02901488  h  log  k 

sg 

2 

R  =  0,998,  Standard  Error  of  Estimation  SEE  *  3% 


EFFECTIVE  b  AT  SUBFA 
OF  SUBBASE,  PSt, 


THICKNESS  OF  BASE  OB  SUBBASE.  IN. 


Figure  2. 1  COMPOSITE  SUBGRADE  MODULUS  OF  REACTION 
FOR  WELL-GRADED  CRUSHED  MATERIAL 
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COMPOSITE  MODULUS  OP  SOIL 


FOR  DIFFERENT  BASE- SUBBASE  MATERIALS 


where:  h  =  layer  thickness,  in. 


k  =  modulus  of  subgrade  reaction,  pci 
sg 


For  natural  sands  and  gravels  (PI<8): 

log  k  =  -1.296084  +  2.263407  log  k  -  0.3013741  (log  k  V 
C  6  sg  6  Sg" 

+  0. 08554373h  -  0.0002574619  h2 
-  0.03050173  h  log  k 


(2.2) 


R  =  0.997 


sg 

SSE  =  3.6“ 


For  stabilized  layers: 

log  k  =-0.1578667  +  1.02813  log  k  +  0.0S44761h 
c  sg 

-  0. 8473852x10" 3h2  +  0. 7254749xl0'6E  (2.3) 

-  0.1937293xl0"12E2  -  0. 4409096xl0"2h  log  k 

&  sg 

-0. 4601653x10” 2E  (logk  )  -0. 2465638xl0_8Eh 

sg" 

R2  =  0.997  SSE  =  2.3% 


where:  E  =  modulus  of  elasticity  of  the  stabilized  layer,  psi. 

It  should  be  stressed  that  equation  2.3  is  only  valid  within  the  range 
of  E  *  105  to  2  x  10^  psi. 

The  approximate  closed  form  solution  was  then  used  with  these  equations 

2 

in  order  to  generate  an  equation  for  the  CV  (k)  value  to  be  used  in  the 
probabilistic  analysis.  (This  value  will  be  discussed  in  Chapter  3.) 

It  should  be  noted  that  the  maximum  number  of  layers  for  which  the 
program  may  be  used  is  four  (including  the  slab  and  the  subgrade).  The 
computation  of  the  k£  for  a  single  base/subbase  layer  is  a  straightforward 
solution.  The  required  values  (subgrade  modulus,  layer  thickness  and  the 
modulus  of  elasticity,  if  applicable)  are  input  into  the  equation  for  the 
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composite  modulus.  If  there  are  two  base/subbase  layers,  the  procedure 
used  is  as  follows' 

a.  The  single  layer  procedure  is  followed  for  the  subgrade 

t 

and  the  layer  directly  above  it  to  obtain  a  k£  value 

I 

b.  This  k  value  is  then  used  as  the  foundation  material 

c 

of  the  upper  base/subbase  layer 

c.  Step  (a)  is  repeated  to  find  a  kc  for  the  entire  subsystem. 

Since  extrapolated  values  were  used  in  the  generation  of  equation  2.3, 
it  is  possible  to  get  a  composite  modulus  greater  than  the  maximum  k£  for 
design  of  S00  pci.  In  the  program,  however,  a  maximum  value  of  500  pci 
is  used.  The  actual  calculated  kc  value  is  printed  as  output  along  with 
an  indication  that  the  maximum  design  k  is  being  used.  These  equations 
have  been  set  up  for  a  pavement  system  which  has  an  increasing  modulus  from 
the  subgrade.  Therefore,  when  analyzing  two  base/ subbase  layers,  the 
stiffer  layer  must  be  on  top. 


There  are  two  approaches  available  for  determining  the  composite  modulus 
of  subgrade  reaction  for  stabilized  layers.  One  method  accounts  for  the  base 


material  by  increasing  the  modulus  of  the  subgrade  while  the  other  method 
accounts  for  the  base  by  using  a  section  of  PCC  with  an  increased  thickness. 

For  the  first  method,  the  composite  k-value  chart  (fig.  2.3)  was  based 
on  calculations  using  the  elastic  layer  theory.  The  values  which  were  ob¬ 
tained  by  using  the  elastic  layer  theory  were  then  corrected  based  on  field 
data  and  experience.  The  value  of  kc  was  determined  by  using  equivalent  de¬ 
flections.  The  k  value  was  raised  until  the  deflection  of  the  "improved" 
subgrade  equalled  that  of  the  base-subgrade  system. 

The  other  method  of  accounting  for  a  base  material  is  the  use  of  an  in¬ 
creased  PCC  thickness.  This  increased  thickness  is  determined  by  the  partially 
bonded  rigid  overlay  pavement  design  equation: 

hdoc  '  1-4'/hd<;1-4  -  (0.0063  V1' 

Eh3 

This  equation  is  based  upon  equivalent  stiffnesses  (D  =  - ■=— )  i.e,,  the 

12(1-0 

stiffness  of  a  PCC  layer  for  a  given  modulus  of  elasticity  and  thickness  is 
equivalent  to  the  stiffness  of  a  base  layer  for  a  given  modulus  of  elasticity 
and  thickness.  From  this,  the  increased  thickness  of  the  PCC  layer  may  be 
calculated. 

In  this  probabilistic  design  methodology,  the  method  of  the  composite 
modulus  of  subgrade  reaction  has  been  used.  It  has  been  chosen  because  it 
yields  the  best  estimate  of  values  to  be  used  in  a  composite  pavement 
section. 
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Chapter  3 

PROBABILISTIC  ANALYSIS  OF  PCC  AIRFIELD  DESIGN 
USING  THE  WESTERGAARD  FREE  EDGE  THEORY 

Introduction 

The  analysis  presented  in  this  chapter  is  based  on  the  approximate 
closed-form  probabilistic  approach  which  was  formulated  in  Volume  II.  The 
use  of  this  closed-form  approach  requires  stress  computations  which  can  be 
made  with  the  H-51  computer  program.  While  the  use  of  the  H-Sl  is  justified 
for  the  pavement  design  for  a  particular  aircraft,  it  can  become  uneconomical 
for  the  probabilistic  approach  due  to  the  large  amount  of  computer  time  needed 
The  analysis  was  therefore  confined  to  the  controlling  aircraft  of  each  USAF 
ACI  group.  Solutions  for  other  aircraft  may  be  developed  by  using  the  user 
defined  option  in  the  program  to  develop  the  regression  constants  to  be  used 
in  the  maximum  tensile  stress  equation. 

The  probabilistic  approach  presented  in  the  following  paragraphs  includes 

(1)  Stress  computations  and  derivation  of  an  equation  to  predict 
the  maximum  tensile  stress  for  USAF  AGI  1  to  USAF  AGI  15; 

(2)  Derivation  of  the  relationship  between  variances  of  the  de¬ 
pendent  and  independent  variables  for  the  approximate  closed- form  approach. 

The  linear  or  first  order  Taylor  series  expansion  is  assumed  as  presented  in 
Volume  II, 

Stress  Computation 

The  Corps  of  Engineers  uses  the  H-51  computer  program  (£)  to  compute 
the  bending  stress  at  the  edge  of  a  slab  supported  by  a  dense  liquid. 

The  program,  developed  by  General  Dynamics  Co.,  numerically  integrates 
the  number  of  blocks  of  the  influence  chart  covered  by  the  load.  The  wheel 


load  is  assumed  to  be  uniformly  distributed  over  an  elliptical  area  shape 
with  any  desired  ratio  of  the  axes.  The  results  of  the  computation  are 
generally  within  2-3  percent  accuracy. 

The  numerical  computation  of  the  stress  constitutes  a  limitation 
for  deriving  a  closed  form  probabilistic  model.  It  was  therefore  necessary 
to  develop  an  equation  for  the  stress  at  the  edge  for  the  load  of  interest. 
The  form  of  the  equation  was  taken  from  the  original  Westergaard's  work. 
Westergaard  C£,5»)  developed  an  approximate  general  formula  for  stresses 
for  the  interior  and  edge  cases  for  any  loaded  area.  The  formula  for  the 
edge  case  reads  as  follows  (6): 


o^  =  tensile  stress  at  the  bottom  of  the  slab  along  the 
edge  or  joint  at  x=y=0. 

P  =  load 

y  =  Poisson's  ratio  of  the  concrete 

h  =  thickness  of  the  slab 

K  and  S  *  area  coefficients  defined  by  the  equations  noted 
A  =  area  of  contact 

r,6  *  radius  and  angle  in  polar  coordinates  describing  the 

distance  of  the  infinitesimal  area  from  the  origin 
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£  =  radius  of  relative  stiffness  given  by  the  equation  noted 

k  =  modulus  of  subgrade  reaction 

y  =  the  distance  from  the  edge  or  joint  to  the  center  of 
gravity  of  the  load. 

According  to  Westergaard's  results.  Equation  3.1  can  be  written  as 
follows,  for  the  elliptical  contact  area: 

o  =  lOnOP  [4  infc  *  f (x,y)  +  1.18(1  +  2V)  £  ]  (3.2) 

6  »(3+p)JT 

where  f(x,y)  is  the  function  of  the  load  configuration  and  geometry.  It 
should  be  observed  that  this  function  is  constant  for  a  given  loading 
(aircraft)  condition. 

It  should  be  stressed  that  equation  (3.1)  was  developed  using  a  finite 
polynomial  series  and  is  restricted  to  the  case  of  a  "small  distance"  of 
the  loaded  area  to  the  origin  at  which  the  stress  is  computed.  Investiga¬ 
tion  of  the  extent  of  this  "small  distance"  for  the  particular  case  of  a 
rectangular  loaded  area  showed  that  the  small  distance  can  be  as  large 

as  about  1.5  t.  Without  this  restriction,  equation  (3.2)  is  a  powerful 
result  which  simplified  the  analysis.  It  is  interesting  to  note  that, 
excluding  the  load  geometry  conditions,  only  two  pavement  variables:  h  and  £, 
determine  the  edge  stress. 

With  the  above  restriction  in  mind,  the  following  general  equation  form 
for  expressing  the  edge  stress  as  function  of  h  and  £  was  adopted: 

o  =  -  (a  +  a  In  t  +  a  /  l  ]  (3.3) 

e  o  1  2 

where  aQ,  a^,  *2  ~  coefficient  dependent  upon  load,  load  configuration 

and  geometry. 
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The  steps  involved  in  computing  the  above  coefficients  for  any  aircraft 
type  are  as  follows: 

(1)  Compute  the  maximum  stress  at  the  edge  induced  by  load  for  different 
l  values.  This  can  easily  be  achieved  using  the  H51  computer  program  by 
choosing  one  set  of  values  of  k,  E  (modulus  of  reaction  of  subgrade  and 
modulus  of  elasticity  of  concrete)  and  varying  h  -  the  slab  thickness. 

(2)  Find  the  coefficients  of  aQ,  a^,  by  the  least  square  error 
method,  over  the  stress  values  computed  in  step  (1). 

The  methodology  is  illustrated  for  the  aircraft  designated  AGI  1  to 
AGI  15.  Their  characteristics  are  shown  in  Table  3.1.  The  results  of  these 
computations  are  shown  in  Appendix  II.  The  stresses  computed  using  eq.  (3.3) 
are  also  shown  in  Appendix  II.  It  can  be  observed  that  they  are  almost  iden¬ 
tical  to  the  original  values  computed  from  the  H-51  computer  program.  In 
general,  correlation  coefficients  greater  than  0.999  were  obtained.  The 
errors  in  computing  the  stresses  with  the  regression  equation  have  been  tab¬ 
ulated  as  a  percentage  of  the  original  stress  and  are  also  presented  in 
Appendix  II.  (The  maximum  error  is-0.8lof  1%).  These  results  clearly  support 
the  use  of  equation  (3.3)  for  a  closed-form  solution  of  the  problem. 

If  the  analysis  is  to  be  used  for  a  user  defined  aircraft  (other  than 
AGI  type),  it  is  necessary  to  input  at  least  five  slab  thickness  levels. 

This  is  necessary  so  that  the  least  square  error  regression  has  enough  data 
points  to  predict  coefficients  which  are  truly  representative  of  the  data. 


Summary  of  USAF  Evaluation 
Characteristics  for  Aircraft  Groups 
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Probabilistic  Closed-form  Approach 
General  Formulation 

The  rigid  pavement  design  is  based  on  computation  of  the  design  factor 
DF  defined  as: 

DF  =  —  (3.4) 

aa 

where  MR  =  modulus  of  rupture  of  concrete 

a  =  load  transfer  coefficient  (=0.75  in  the 
Corps  of  Engineers  procedure) 
a  =  free  edge  stress  induced  by  load  (Westergaard 
solution) 

The  design  factor  is  related  to  the  number  of  coverages  on  the  basis 
of  USACE  test  section  results  and  experience.  The  relationships  are  shown 
in  Figs.  3-la  to  Fig.  3-lc  for  three  levels  of  failure  mode:  (i)  initial 
failure,  (ii)  shattered  slab  condition  and  (iii)  complete  failure.  There¬ 
fore,  knowledge  of  the  design  factor  and  its  probability  distribution 
leads  to  a  complete  description  of  the  design  in  probabilistic  and  relia¬ 
bility  terms.  In  the  following,  the  closed  form  approach  is  described. 

In  equation  (3.4)  the  stress  is  expressed  as  function  of  the  material 
and  slab  geometry  variable,  (E, u  ,  k  and  h)  using  equation  (3.3).  In  the 
analysis,  a  constant  y  of  0.15  was  used.  It  has  been  held  as  a  constant  since 
there  is  a  relatively  insignificant  effect  of  y  on  o.  Also,  H-51  uses  a 
equal  to  0.15  in  order  to  generate  the  stress  predictor  coefficients.  There¬ 
fore,  u  is  equal  to  a  constant  value  of  0.15  instead  of  being  stoichastic. 
Therefore,  the  average  design  factor  is: 

MR _ 

— —  S  ■  «  —  1  —  — 

DF  n  .  o(I,  k,  F) 
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OES'GN  r*CTOP 


and  its  variance  is: 


n 

E 

jx)2  Var[X.] 

i=l 

n 

n 

r 

i=l 

L 

i=i 

^lxk>  c°v(VV  '5-6> 

k?<l 

where  X^  -  denotes  the  different  variables,  MR,  a,  E,  k,  h 

-  Var[X.1  -  variance  of  variable  X. 

iJ  1 

Covar  [X.,Xkl  -  covariance  of  variables  X^  and  X^  which  is  zero 
when  these  variables  are  independent,  and  equals 

p  /"varfx.  ]  "Var  "  when  the  variables  are  dependent, 
with  a  correlation  coefficient  p  . 

The  derivatives  of  DF  with  respect  to  all  variables  can  be  found 
from  eq.  3.3  and  3.4  using  the  chain  rule.  For  example,  the  derivative 
of  DF  with  respect  to  E  is  given  by 


and 


3  DF 

.  3DF 

3o  _  MR 

3  E 

3o 

*  3E  '  2 

ao 

3o 

3o 

.  3l  =  al-a2/t 

5E 

3* 

3  E  h2^ 

3o 
3  E 


4E 


leading  to 


Similarly: 


3  DF 

MR 

Te 

“  "  2 

ao 

3  DF 

MR 

TT 

ao“? 

3  DF 

_  MR 

3  h 

~  y 

ao  4 

(1  -  *2/1.-) 
4h2E 


(al-^2/i) 

4h2E 


2o  .  3Cai-a2/X) 

^  4? 


.  h 


(3.7) 
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r 


TT 


and: 


3DF  =  _J_ 
3MR  ~  ac 

3DF  _  MR 
oo  2 


a  a 

For  the  following,  only  MR  and  E  are  assumed  to  be  dependent.  If  the  co¬ 
efficient  of  variation,  CV,  is  defined  as  the  ratio  of  the  standard  deviation 
to  the  mean  value  of  the  variable,  then  substituting  eq  (3.7)  into  eq.  (3.6) 
and  rearranging  the  terms  lead  to: 


CV2  [DF]  =  CV2[MR]  ♦  CV2[a]  ♦ 


1  ai-a2/t  1 

l  4ch2  J 

[-  *  l£t 


(CV2[E]  ♦  CV2[k])  ♦ 


2/1) 


CV^h] 


(3.8) 


'  al-a2/t  | 

.  4<jh2  J 


p [MR, E]  .  CV[MR]  .  CV[E] 


where  CV[X.]  is  the  coefficient  of  variation  of  variable  X^.  In  this 

2  2 
equation,  the  value  of  CV  [k]  is  either  the  input  value  of  CV  [k]  or  it 

2 

is  the  calculated  value  of  CV  [k]  if  the  composite  modulus  of  subgrade 

2 

reaction  is  calculated.  These  CV  [k]  values  may  be  calculated  for  the 
materials  found  in  equations  2.1  to  2.3  by  using  the  first-order  Taylor 
series  approximation  (from  Vol.  II). 

For  well-graded  crushed  materials: 

Var [k  ] 


CV2[kc] 


F  2 

C 


5.3018981  Jh2 (0,08901252  -  0. 0008850388h 

-  0.02901488  log  k  )2CV2fh] 

(0.9640174  -  0.2561922  log  k 

-  0.012601  h)2  CV2[k  ]] 

Sg  J 


sg 


(3.9) 


where: 


CV(i] 


h 


k 


sg 


coefficient  of  variation  of  variable  i 
composite  modulus  of  subgrade  reaction,  pci 
layer  thickness,  in. 
modulus  of  subgrade  reaction,  pci 
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For  natural  sand  and  gravel: 


CV2[k  ] 
1  c 


*  5.3018981  £h2(0. 08554373  -  ,0005149238h  -  0.03050173  logksg)2CV2[h] 

+  (0.9829852  -  0.2617702  log  k  -  0.0132467h)2CV2[k  ]]  (3.10) 


For  stabilized  layers: 

CV2[k  ]  =  5.3018931  [h2(0. 0544761  -  1. 6947704x10" 3h  -  0.4409096xl0'2log  k 

c  ll  ^8 

-0. 2465638x10 CV2{h]  ♦  E2t0.7254749xl0-6  -  0. 3874568xl0*12E 

-  0.4601653xl0"7log  k  -  0.2465638xl0"8h) 2  CV2[El 

s?  (3. 

+  (0.4465112-0. 1914846x10" 2h  -  0. 1998463xl0'7E)2  CV2[k  jj 

sg  J 

where  E  =  modulus  of  elasticity  of  the  stabilized  layer,  psi 


Equations  3.5  and  3.8  are  used  to  compute  the  average  DF  and  its  co¬ 
efficient  of  variation,  using  (i)  equation  3.3  and  (ii)  material  characteris¬ 
tics. 

These  equation  can  easily  be  numerically  transformed  into  different 
reliability-number  of  coverage  -  slab  thickness  relationships.  Two  possi¬ 
bilities  of  implementing  the  above  formulation  are  discussed  in  the  next 
section. 


Reliability-Number  of  Coverages  Relationship 
for  Different  Slab  Thicknesses 

With  the  design  factor  assumed  normally  distributed,  the  number  of 
coverages  corresponding  to  DF(1  k.CV[DFJ)  can  be  computed,  and  the 
probability  of  Df£dF(1  *  kCV[DF])  is  taken  from  the  normal  distribution. 

In  order  to  do  this,  equations  had  to  be  developed  for  the  DF-Number 
of  Coverages  curve  (Figure  3.1a,  since  initial  failure  is  assumed). 
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Since  each  curve  is  bilinear,  two  equations  were  derived  for  each. 

The  first  equation  defines  the  curve  below  5000  coverages  and  the  second 
equation  defines  the  rest  of  the  curve.  Since  these  equations  were  de¬ 
veloped  independent  of  the  modulus  of  subgrade  reaction,  ranges  had  to 
be  assumed  for  each  curve.  These  ranges  are: 

Range  of  k 

Curve  Assigned  to  Curve 


25  -  200 

25  - 

250 

300 

250  - 

350 

400 

350  - 

450 

500 

450  - 

500 

This  part  of  the  design  procedure  yields  the  reliabilities  for  a 
given  slab  thickness  in  terms  of  the  number  of  coverages  it  will  take  to 
produce  that  reliability. 

Reliability-Slab  Thickness  Relationship 
for  Different  Number  of  Coverages 

This  part  of  the  design  procedure  gives  the  reliability  as  a  function 
of  the  slab  thickness  for  a  given  number  of  coverages.  This  can  be  done 
by  basically  working  the  previous  procedure  in  the  reverse  order.  If  the 
desired  number  of  coverages  is  known,  then  the  DF  may  also  be  calculated 
by  using  the  equations  described  in  the  previous  section.  The  probability 
is  then  taken  from  the  normal  distribution.  The  information  obtained 
from  this  analysis  can  be  very  useful.  If  a  certain  reliability  is  de¬ 
sired  for  a  given  coverage  level,  then  the  necessary  slab  thickness  can 
be  found. 


22 


Study  of  Effects  of  the  Variation  of  the  Design  Parameters 


The  closed  form  solution  has  the  advantage  to  bring  an  insight  of  the 
effect  of  the  design  parameters.  In  eq.  (3.8),  the  coefficient  of  variation 
of  DF  and  hence  the  reliability  of  the  design,  is  expressed  as  a  function 
of  the  coefficient  of  variation  of  the  design  parameters.  The  effect  of 
these  parameters  is  enhanced  if  equation  (3.8)  is  written  as  follows: 

CV2(DF]  =  Z  vi  .  CV2[X.]  + 

1  ?i  *  plXi*  h]  *  CV[XiJ  *  CV  ^k1 

(3.12) 

al  ”  a2^ i  2 

where:  w.,£.  -  weighting  factors  (w.  =  ( - - - ) 

11  1  4oh 

for  the  variable  E,  w^  =  1  for  the  variable  MR,  etc.) 


The  effect  of  each  parameter  is  either  increased  or  reduced  by  the 
corresponding  weighting  factor  which  comes  from  the  response  function.  It 
should  be  noted  that  the  weighting  factor  is  dependent  upon  load  and 
material  characteristics  and  slab  thickness.  An  evaluation  of  the  weight¬ 
ing  factors  for  the  thirteen  Aircraft  Group  Indices  is  presented  in 
Appendix  II,  over  the  practical  range  of  l  values.  It  is  seen  that  the 
weighting  factors  of  E  and  k  (the  modulus  of  elasticity  of  the  concrete 
and  the  subgrade  modulus  of  reaction)  are  quite  small,  resulting  in  a  minor 
effect  of  their  variability  on  the  reliability  of  the  design.  The  weighting 
factor  of  the  slab  thickness  is  the  largest  one,  with  the  result  of  ampli¬ 
fying  the  effect  of  thickness  variability  on  DF  variability.  However,  it 
should  be  remembered  that  thickness  variability  is  relatively  small.  It 
seems  that  the  leading  parameters  controlling  DF  variability  are  thus  the 
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modulus  of  rupture  and  load  transfer  variabilities.  With  the  increase 
of  number  of  wheels  (from  light  to  heavy- load  design),  the  contribution 
of  the  variability  of  E  and  k  increases. 

Run  Example 

The  new  approach  was  used  to  develop  reliability-number  of  coverages 
curves  for  different  slab  thicknesses  (11,  12,  13,  14,  IS  and  16  inch) 
for  AGI  13.  Reliability-slab  thickness  curves  were  then  developed 
for  different  coverage  levels  (1000,  5000,  10,000,  20,000,  50,000,  and 
100,000  coverages).  The  average  values  of  the  variables  and  their  coef¬ 
ficients  of  variation  are  given  in  Figure  3,2a.  Output  may  be  found  in 
Figure  3.2b.  Graphical  results  (plots  of  Reliability-Coverages-Slab 
Thickness  relationships)  may  be  found  in  Figure  3.3  and  3.4. 
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Reliability-Log  Coverages  Curves  for 
Different  Slab  Thicknesses 


N=100 


Reliability  -  Slab  Thickness  Curves 
for  Different  Coverage  Levels 


Chapter  4 


SUMMARY 

The  current  U.S.  Army  and  U.S.  Air  Force  Design  Manuals  use  the 
Westergaard  free  edge  stress  slab  theory  as  the  basis  for  their  rigid 
pavement  design  methodologies.  The  design  procedure  has  been  expanded 
to  provide  a  solution  expressed  in  probabilistic  and  reliability  terms. 
Further  developments  were  required  in  the  original  procedure  to  make  the 
analysis  more  practicable.  Two  major  investigations  were:  (1)  The  use 
of  the  composite  modulus  of  subgrade  reaction  to  expand  the  procedure  to 
solve  problems  having  multiple  subbase  layers,  and  (2)  The  evaluation 
of  a  general  equation  form  used  to  predict  maximum  tensile  stresses  at 
the  bottom  of  the  concrete  slab  for  aircraft  types  designated  by  USAF 
AG I  1-15. 

SPSS  Multiple  Regression  runs  were  made  to  generate  the  composite 
modulus  of  subgrade  reaction  equations  for  three  different  material 
types  currently  used  in  military  manuals.  The  correlations  between  the 
curve  data  and  the  equations  that  fit  this  data  are  very  close  to  1 
(R2  =  0.998,  0.997,  and  0.997). 

A  regression  equation  was  developed  to  predict  the  maximum  tensile 
stress  at  the  bottom  of  the  concrete  layer.  Regression  constants  were 
then  determined  for  each  of  the  1  USAF  Aircraft  Group  Index  controlling 
aircraft.  The  stresses  predicted  by  this  equation  are  extremely  close 
to  those  predicted  by  the  H-51  and  well  within  the  bounds  of  an  acceptable 
margin  of  error. 


The  derivations  above  have  been  included  in  a  computer  program  for 
the  probabilistic/reliability  analysis  of  rigid  pavements.  The  approximate 


closed-form  probabilistic  approach  (Taylor  Series)  has  been  utilized. 

The  computer  program  can  be  used: 

(1)  in  the  analysis  in  probabilistic/reliability  terms  for  a  given 
pavement  system  and  loading  aircraft.  The  means  and  coefficients  of 
variation  of  the  design  parameters  serve  as  the  input.  The  computer  pro¬ 
gram  produces  values  of  the  number  of  coverages  and  their  reliability 
levels. 

(2)  in  the  design  of  a  rigid  pavement  for  a  given  reliability. 

Here,  the  number  of  coverages  is  also  an  input  so  that  the  slab  thickness 
may  be  determined  for  the  given  reliability  and  coverage  level. 
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SPSS  Outputs  for  the  Composite  Modulus  Equations 
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igure  1-3  -  SPSS  Multiple  Regression  Output  for  th 
Subgradc  Reaction  of  Stabilized  Layers 
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Regression  Constants  for  Free  Edge  Stress  Equations 
(AGI  1  to  AG I  13) 
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Table  11.1  -  Summary  of  Stress  and  Weighting  Factor  Computations 
for  AGI-1/C-12S 


h 

(in.) 

t, 

(in.) 

0  - 

(psi) 

Weighing 

Factors 

Error  as 
%  of  H-51 

H-51 

EQ 

am 

warn 

5E,MR 

■HBjfl 

26.47 

1212.38 

1212.56 

1.0000 

0.0351 

2.0681 

-0.3746 

♦0.01 

29.72 

968.94 

968.38 

1 .0000 

0.0300 

2.1916 

-0  .  3464 

-0.06 

8.0 

32. 8S 

791.94 

792.98 

1 .0000 

0.0265 

2.2840 

-0.3258 

♦  0.15 

9.0 

35 . 88 

663.21 

662.59 

1.0000 

0.0238 

2.3636 

-0.3084 

-0.09 

10.0 

38.83 

563.20 

562.93 

1.0000 

0.0218 

2.4258 

-0.2950 

-0.05 

11.0 

41.71 

484.50 

484.90 

1.0000 

0.0202 

2.4765 

-0.2842 

♦  0.08 

12.0 

44. S2 

422.35 

422.55 

1.0000 

0.0188 

2.5230 

-0.2744 

♦  0.05 

13.0 

47.27 

372.17 

371.89 

1.0000 

0.0177 

2.5642 

-0.2658 

-0.07 

14.0 

49.98 

330.12 

330.18 

1.0000 

0.0167 

2.5979 

-0.2588 

♦  0.02 

a  *  -  70026  R  *  1.000 

o 

ax  =  34  374  R2  »  1.000 

a  *  44166 

2 

k  *  150  pci 
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Table  11.2 


Suanary  of  Stress  and  Weighting  Factor  Computations 
for  AGI-2/F-4 


26.47 

29.72 

32.85 
3S.88 

38.85 
41.71 
44.  S2 
47.27 
49.98 


0  -  (psi) 

Neighing 

Factors 

H-Sl 

EQ 

w 

ll 

CE,MR 

1676.00 

1673.95 

1 .0000 

0.0221 

2.4165 

-0.2970 

1311.20 

1312.76 

1.0000 

0.0179 

2.5565 

-0.2674 

10S6.24 

1057.95 

1.0000 

0.0151 

2.6611 

-0.2458 

872.14 

871.61 

1 .0000 

0.0151 

2.7450 

-0.2288 

731.12 

731.24 

1.0000 

0.0117 

2.8090 

1 

-0.2160 

623.84 

622.81 

1 .0000 

0.0105 

2.8646 

-0.2050 

537.68 

537.22 

1.0000 

0.0096 

2.9084 

-0.1964 

468.09 

468. 4S 

1.0000 

0.0089 

2.9453 

-0.1892 

412.12 

412.37 

1 .0000 

0.0084 

2.9784 

-0.1828 

*2  "  * 
k  *  15 


Table  11.3  -  Summary  of  Stress  and  Weighting  Factor  Computations 
for  AGI-3/F-111 


h 

t 

0  • 

(P*i) 

Weighing  Factors 

Error  as 
t  of  H-&1 

(in.) 

(in.) 

H-Sl 

EQ 

esa 

*E,MR 

20.47 

2122.98 

2121.58 

1.0000 

0.0318 

2.146S 

-0.3566 

-0.07 

IF^B 

29.72 

1685.46 

1687.93 

1.0000 

0.0276 

2.2560 

-0.3520 

♦0.15 

8.0 

32.85 

1378.84 

1373.34 

1 .0000 

0.0244 

2.3461 

-0.5122 

-0.04 

9.0 

35.88 

1150.03 

1149.21 

1.0000 

0.0221 

2.4165 

-0.2970 

-0.07 

10.0 

38.83 

974.10 

974.65 

1 .0000 

0.0203 

2.4728 

-0.2850 

♦  0.06 

11.0 

41.71 

838.07 

838.36 

1.0000 

0.0188 

2.5230 

-0.2744 

♦  0.03 

12.0 

44.52 

730.33 

729.69 

1.0000 

0.0176 

2.5670 

-0.2652 

-0.09 

13.0 

47.27 

641.42 

641.55 

1.0000 

0.0166 

2.6027 

-0.2578 

♦  0.02 

14.0 

49.98 

S68.94 

569.08 

1 .0000 

0.0158 

2.6358 

-0.2510 

♦0.03 

*o 

* 

-115728 

R  - 

1.000 

*1 

* 

57671 

R2  * 

1.000 

*2 

m 

84016 

k 

u 

150  pci 
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Table  11.4  -  Suanary  of  Stress  and  Weighting  Factor  Confutations 
for  AG1-4/C-130 


h 

t 

0  • 

(psi) 

Weighing 

Factors 

Error  as 
\  of  H-Sl 

(in.) 

(in.) 

H-51 

EQ 

W  *W 
m  "a 

"EL-Wk 

N, 

h 

CE,MR 

7.0 

26.15 

762.33 

763.86 

1.0000 

0.0463 

1.8344 

-0.4304 

♦  0.20 

8.0 

28.91 

641.88 

639.96 

1.0000 

0.0537 

1.7028 

-0.4O54 

-0.30 

9.0 

•31.58 

549.96 

549.82 

1.0000 

0.0578 

1.6353 

-0.4808 

-0.03 

10.0 

34.17 

481.24 

480.68 

1 .0000 

0.0590 

1.6162 

-0.4858 

-0.12 

11.0 

36.71 

425.55 

425.94 

1.0000 

0.0590 

1.6162 

-0.4858 

♦  0.09 

12.0 

39.18 

380.20 

381.17 

1 .0000 

0.0582 

1.6292 

-0.4824 

♦  0.2b 

13.0 

41.61 

343.93 

344.03 

1.0000 

0.0565 

1.6569 

-0.4752 

♦  0.03 

14.0 

43.98 

312.90 

312.56 

1 .0000 

0.0547 

1.6864 

-0.4676 

-0.11 

15.0 

46.32 

285. 88 

285.71 

1.0000 

0.0530 

1.7145 

-0.4604 

-0.06 

a  «  -332671  R  =  1,000 

o 

ax  *  94339  R2  -  1.000 

a2  *  1626321 


k  »  230  pci 
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Table  II. S  -  Summary  of  Stress  and  Weighting  Factor  Confutations 
for  AGI-5/C-  9 


h 

(in.) 

l 

(in.) 

0  - 

(psi) 

Weighing  Factors 

Error  as 
%  of  H-Sl 

H-51 

EQ 

W  =W 

Ml  a 

N  *W 

EL  k 

£e,mr 

7.0 

26.15 

1240.10 

1238.26 

1.0000 

0.0496 

1.7748 

-0.4452 

1 

o 

M 

cn 

8.0 

28.91 

1033.77 

1033.08 

1 .0000 

0.0422 

1.9149 

-0.4108 

-0.07 

9.0 

51.58 

873.87 

875.66 

1 .0000 

0.0372 

2.0209 

-0.3856 

♦  0.21 

10.0 

34.17 

749.80 

752  .  37 

1 .0000 

0.0333 

2.1098 

-0.365  0 

♦0*  34 

12.0 

39.18 

575.56 

574.63 

1 .0000 

0.0275 

2.2569 

-0.3318 

-0.16 

13.0 

41.61 

509. 9S 

509.23 

1 .0000 

0.0256 

2.3113 

-0.3198 

-0.14 

14.0 

43.98 

455.69 

454.67 

1.0000 

0.0239 

2.3608 

-0.3090 

-0.22 

16.0 

48.62 

369.76 

369.77 

1.0000 

0.0214 

2.4380 

-0.2924 

♦0.004 

17.0 

50.88 

335.62 

336 . 26 

1.0000 

0.0204 

2.4690 

-0.2858 

♦  0.19 

aQ  *  -127S82  R  =  1.000 

=  56841  R2  =  1.000 

a2  =  71556 

k  »  250pci 
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Table  II. C  -  Summary  of  Stress  and  Weighting  Factor  Computations 
for  AGI-6/T-43 


26.  IS 
28.91 
31.  S8 

34.17 

39.18 

41.61 
43.98 

48.62 
SO.  88 

-183723 
69928 
434391 
250  pci 


o  -  (psi) 


1247.20 

■1041.67 

880.58 

759.96 

580.79 

518.12 

463.69 

378.62 


1247.41 
1039.90 
882.23 
759.30 
582.44 
517. 3S 
462.97 
378.18 


R“  =  1.000 


Weighing  Factors 


U 

m  a 

*EL-"k 

"h 

€e,mr 

1 .0000 

0.0476 

1.8109 

-0.4  3b2 

1.0000 

0.0424 

1.9108 

-0.4118 

1 .0000 

0.0388 

1.9861 

-0.3938 

1.0000 

0.0354 

2.0604 

-0.3764 

1.0000 

0.0309 

2.1677 

-0.3518 

1.0000 

0.0288 

2.2219 

-0.  339(> 

1 .0000 

0.0273 

2.2632 

-0. 5304 

1.0000 

0.0247 

2.3351 

-0. 5140 

i  nnnn 

0.0238 

2  3 1*22 

-0-3086 

Error  as 
*»  of  H-Sl 


♦0.02 

-0.17 

♦0.19 

-0.09 

♦0.28 

-0.15 

-0.16 

-0.12 


Table  11.7  -  Summary  of  Stress  and  Weighting  Factor  Confutations 
for  AG1-7/B-727 


h 

(in.) 

t 

(in.) 

0  - 

(psi) 

Neighing  Factors 

Error  as 
t  of  H-Sl 

H-51 

EQ 

mm 

^E.MR 

mmm 

20.  S8 

1415.27 

1414.53 

1.0000 

0.0517 

1.7366 

-0.4S4S 

-0.05 

29.03 

1207.37 

1209.18 

1.0000 

0.04  80 

1 . 8028 

-0.4332 

♦  0.15 

11.0 

33.75 

920.14 

919.40 

1.0000 

0.0415 

1.9282 

-0.407b 

-0.08 

13.0 

38.25 

728.09 

726.88 

1 .0000 

0.0368 

2.0286 

-0.3838 

-0.17 

15.0 

42.58 

589.67 

591.66 

1 .0000 

0.0356 

2.1037 

-0, 3664 

♦0.34 

17.0 

46.77 

493.45 

492.57 

1.0000 

0.0304 

2.1818 

-0  .  34  8b 

-0. 18 

19.0 

50.84 

417.62 

417.52 

1.0000 

0.0282 

2.2389 

-0.33SR 

-0.02 

20.0 

52.84 

386.65 

386.62 

1.0000 

0.0272 

2.2659 

-o.s:9s 

-0.01 

21.0 

54.81 

359.03 

3S9.17 

1 .0000 

0.0263 

2.2904 

-0.3244 

♦  0.04 

a  =  -348726  R  *  1.000 

o 

ax  *  121 87S  R2  >  1.000 

a2  *  1049552 
k  *  350  pci 
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Table  II. 8  -  Summary  of  Stress  and  Weighting  Factor  Computations 
for  AG1-8/E-3 


h 

lin.) 

1 

Un.) 

0  - 

(psi) 

Weighing  Factors 

H-Sl 

EQ 

EES9 

CE,MR 

8.0 

26.58 

1006.29 

1003.38 

1.0000 

0.0525- 

1.7224 

-0.4584 

9.0 

29.03 

869.08 

863.61 

1.0000 

0.0668 

1.5001 

-0.5168 

11.0 

33.75 

687.93 

685.33 

1.0000 

0.0789 

1.3393 

-0.5618 

13.0 

38. 2S 

565.64 

565.66 

1 .0000 

0.0807 

1.3179 

-0.5680 

15.0 

42. S8 

478.17 

479.41 

1 .0000 

0.0778 

1.3526 

-0.5580 

17.0 

46.77 

413.21 

413.74 

1.0000 

0.0731 

1.4152 

-0.5408 

19.0 

50.84 

362.05 

361.93 

1 .0000 

0.0682 

1 .4  796 

-0.5224 

20.0 

52.84 

340.47 

339.97 

1.0000 

0.0658 

1.5141 

-0.5130 

21.0 

S4.81 

320.04 

320.08 

1.0000 

0.0639 

1.5423 

-0.5054 

a  »  -823985  R  =  1.000 

o 

ax  =  221392  R2  =  1.000 

a2  =  4314471 
k  =  350  pci 
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Table  11.9  -  Summary  of  Stress  and  Weighting  Factor  Confutations 
for  AGI-9/C-141 


h 

(in.) 

t 

(in.) 

0  - 

(psi) 

Weighing  Factors 

H-Sl 

EQ 

E59 

IMB1E 

"h 

5e,mr 

8.0 

2o.S8 

1027.58 

1029.53 

1.0000 

r 

0.0783 

1.34b3 

-0.5598 

9.0 

29.03 

902.06 

900.45 

1 .0000 

0.0866 

1.2479 

-0.5886 

11.0 

33. 7S 

722. b4 

721.84 

1 .0000 

0.0897 

1.2133 

-0.5990 

13.0 

38.25 

598.93 

599.54 

1 .0000 

0.0854 

1 .2614 

-0.584b 

1S.0 

42.58 

509.12 

509.24 

1 .0000 

0.0789 

1.3393 

-0,5618 

17.0 

46.77 

439.47 

439.61 

1  .0000 

0.0726 

1.4197 

-0.5390 

la.o 

50.84 

384 . 38 

384.35 

1.0000 

0.0669 

1.4987 

-0.5172 

20.0 

52.84 

360.99 

360.87 

1.0000 

0.0643 

1.5364 

-0.5070 

21.0 

54.81 

339.60 

339.58 

1.0000 

0.0620 

1.5703 

-0.4978 

a  * 
0 

-802693 

R 

=  1.000 

ai  * 

220139 

R2 

=  1.000 

a2  * 

3893752 

k  * 

350  pci 
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Table  11.10  Summary  of  Stress  and  Weighting  Factor  Con>utations 
for  AGI-1U/C-SA 


h 

1 

0  - 

Cpsi) 

Weighing  Factors 

Error  as 
\  of  H-Sl 

(in.) 

(in.) 

H-Sl 

EQ 

HEM 

PCTI 

mm 

^E,MR 

8.0 

26.58 

838. IS 

842.22 

♦0.49 

9.0 

29.03 

719.03 

715.91 

-0.43 

11.0 

33.75 

556.58 

554.84 

Weighting  factors  not  given 
due  to  the  non-symraetric 

-0.31 

13.0 

38.25 

454.00 

4S3.11 

nature  of  the  C-5A 
gear  configuration 

-0.20 

15.0 

42.58 

379.20 

381.48 

*■0.60 

17.0 

46.77 

327.63 

327.78 

♦0.05 

19.0 

50.84 

286.07 

285.88 

-0.07 

70.0 

52.84 

268.25 

268.23 

-0.01 

21.0 

54.81 

2S2.59 

252.29 

-o.i: 

a  *  -653058 
o 

aj  =>  174658 

*  3563166 

k  *  350  pci 
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Table  II. 11  -  Summary  of  Stress  and  Weighting  Factor  Confutations 
for  AGI-ll/KC-10 


h 

t 

0  - 

Ipsi) 

Weighin’ 

Factors 

Error  as 
\  Of  H-Sl 

(in.) 

(in.) 

H-51 

HQ 

iza 

EIPH 

"h 

^E.MR 

8.0 

2o.58 

1089.90 

1094.80 

1 .0000 

0.0344 

2.0846  j 

-0.3708 

♦0.45 

9.0 

29.03 

935  .  34 

930.85 

1.0000 

0.0517 

1.7574 

-0.4546 

-0.48 

11. U 

33.75 

720. b4 

725.47 

1.0000 

0.0727 

1.4182 

-0.5394 

-0.16 

13.0 

33. 2S 

590.92 

596.81 

1 .0000 

0.0794 

1.3531 

-0.5636 

-0.02 

15.0 

42.58 

505.14 

S05.94 

1 .0000 

0.0794 

1.3338 

-0.5634 

♦0.16 

17.0 

40.77 

436.00 

437.35 

1  1.0000 

0.0763 

1.3722 

-0.5524 

♦0.17 

19.0 

50.84 

383.50 

383.43 

1.0000 

0.0720 

1.4275 

-0.5568 

-0.02 

20.0 

52.84 

300.37 

3b0 . 5  8 

1 .0000 

0.0701 

1 .4542 

-0.5294 

♦  0.06 

21.0 

54.81 

340.47  i 

339.86 

1-0000 

BWffWH 

.  ,-P.U— . 

a  = 

-968715 

R  -  1.000 

0 

ai = 

254788 

R2  =  1.000 

a2  ' 

5396732 

k  = 

350  pci 
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Table  11.12  -  Summary  of  Stress  and  Weighting  Factor  Computations 
for  AGI-12/E-4 


h 

t 

o 

1 

(psi) 

Weighing  Factors 

Error  as 
\  of  H-51 

(in.) 

(in.) 

H-51 

EQ 

***« 

E8BW 

"h 

4, MR 

9.0 

26.55 

833.52 

839.32 

1.0000 

0.0377 

2.0090 

-0.3884 

♦0.70 

11.0 

30.87 

650.71 

646.26 

1.0000 

0.0683 

1.4789 

-0.5226 

-0.68 

13.0 

34.99 

536.25 

S31.93 

1.0000 

0.0823 

1.2987 

-0.5736  | 

-0.81 

1S.0 

38.95 

452.55 

453.09 

1.0000 

0.0875 

1.2379 

-0.5916 

♦  0.12 

17.0 

42.78 

392.51 

394.00 

1.0000 

0.0866 

1.2479 

-0.5G86 

♦  0.  38 

19.0 

46.51 

345.70 

347.59 

1 .0000 

0.0834 

1.2850 

-0.577b 

♦  0.53 

21.0 

50.13 

308.94 

309.80 

1.0000 

0.0788 

1.3407 

-0.5614 

♦  0.28 

23.0 

53.67 

278.86 

278.49 

1.0000 

0.0739 

1.4026 

-0.5438 

-0.13 

24.0 

55.41 

265.82 

264.74 

1.0000 

0.0715 

1.4  34  7 

-0.5343 

-0.41 

aQ  -  -101891$  R  =  1.000 

aj  «  266253  R2  *  1.000 

a^  *  5077764 

k  *>  500  pci 
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Table  11.13  -  Summary  of  Stress  and  Weighting  Factor  Computations 
for  AG1-13/B-52 


Table  II  -  14  -  Summary  of  Stress  and  Weighting  Factor  Computations 

For  AG I  -  14/0V-1 


h 

(in.) 

l 

(in.) 

0  - 

(psi) 

Weighting  Factors 
_ 

Error  as  | 
\  of  ll-Sl! 

i 

i 

H-51 

EQ 

ma 

m 

ee,mr 

26.47 

561.83 

563.00 

1.0000 

0.0147 

2.6758 

-0.2420 

i 

+0.21 

1 

29.72 

437.56 

436.77 

1.0000 

0.0130 

2.7S10 

-0.2276 

-0.18 

8.0 

32.85 

350.14 

3S1.02 

1.0000 

0.0117 

2.8060 

-0.2166 

♦0.25 

9.0 

3S.88 

287.58 

286.81 

1.0000 

0.0108 

2.8514 

-0.2076 

-0.27 

10.0 

38.83 

239.66 

239.94 

1.0000 

0.0101 

2.8839 

-0.2012 

■ 

♦0.12 

11.0 

41.71 

203.41 

204.00 

1.0000 

0.0095 

2.9135 

-0.1954 

♦0.29 

12.0 

44.52 

175.35 

175.76 

1.0000 

0.0090 

2.9412 

-0.1900 

♦0.24 

13.0 

47.27 

153.13 

153.13 

1 .0000 

0.0086 

2.9670 

-0.1850 

♦0.02 

14.0 

49.98 

135.17 

134.78 

1 .0000 

0.0081 

2.9908 

-0.1804 

-0.29 

ao  “ 

-9544 . 

R  «  0.999 

*1  * 

9268. 

R2=  1.000 

a2  * 

14558. 

k  « 

150  psi 
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Table  II  -  IS  -  Summary  of  Stress  and  Weighting  Factor  Computations 

For  AGI  -  1S/C-S4 


h 

(in.) 

9 

(in.) 

0  - 

(psi) 

Weighting 

Factors 

B 

H-Sl 

EQ 

B5S 

■sa 

n 

EC  .MR 

7.0 

26.15 

635.12 

635.23 

1.0000 

0.0301 

2.1880 

-0.3472 

♦0.02 

8.0 

28.91 

520.93 

520.87 

1.0000 

0.0283 

2.2362 

-0.3364 

-0.01 

9.0 

31.58 

436.59 

436.38 

1.0000 

0.0266 

2.2813 

-0.3264 

-0.05 

10.0 

34.17 

372.55 

371.86 

1.0000 

0.0252 

2.3232 

-0.3172 

-0.18 

11.0 

36.71 

320.77 

321.48 

1 .0000 

0.0241 

2.3553 

-0.3102 

♦0.22 

12.0 

39.18 

280.57 

281.12 

1 .0000 

0.0230 

2.3886 

-0.3030 

♦0.20 

13.0 

41.61 

248.55 

248.31 

1.0000 

0.0218 

2.4239 

-0.2954 

-0.10 

14.0 

43.98 

221.79 

221.15 

1.0000 

0.0209 

2.4549 

-0.2888 

-0.29 

15.0 

46.32 

198. IS 

198.46 

1 .0000 

0.0202 

2.4756 

-0.2844 

♦0.16 

a„  ■  76186.  R  -  1.000 

o 

Ij  ■  30235.  R2.  1.000 

a2  -  225665. 
k  >  250  psi 
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APPENDIX  III 


USER'S  GUIDE 


III 


USER'S  GUIDE 


Card  1  (12,  13,  4F10.4) 


1-2 

IG 

3-5 

NOLBP 

6-15 

XK 

16-25 

CVK 

26-35 

ETEMP 

36-45 

XMUTEM 

-  USAF  Aircraft  Group  Index  (=0  for  user  defined  aircraft) 

-  Number  of  base  (subbase)  layers  (excludes  subgrade,  maximum  of  2) 

-  Modulus  of  subgrade  reaction,  pci 

-  Coefficient  of  variation  of  XK 

-  Modulus  of  elasticity  of  the  pavement,  psi 

-  Poisson's  ratio  of  the  pavement 


Card  2  (12 
1-2 


3-12 

13-22 

23-32 

33-42 


,  4F10.4)  Repeat  NOLBP  times  *See  notes 
MATOPT  -  Material  of  the  base  (subbase)  layer 
=  1  for  well-graded  crushed  materials 
=  2  for  natural  sands  and  gravels  (PI  8) 

=  3  for  stabilized  materials 
HB  -  Thickness  of  base  (subbase)  layer,  in. 

CVHB  -  Coefficient  of  variation  of  HB 

EB  -  Modulus  of  elasticity  of  base  (subbase)  layer,  psi  (only 

required  when  MATOPT  *  3) 

CVEB  -  Coefficient  of  variation  of  EB  (only  required  when  MATOPT  =  3) 


Card  3  (10A6) 
1-60 


Required  if  IG  =  0  *See  notes 

-  Comments  (Only  input  if  aircraft  is  user  defined) 


Card  4  (4F10.0) 
1-10  G 
11-20  P 
21-30  A 
31-40  XNOG 


Required  if  IG  *  0 

-  Load  on  gear,  lb. 

-  Tire  inflation  pressure,  psi 

-  Contact  area  of  one  tire,  sq.  in. 

-  Number  of  Y  (GAMMA) 


III-l 


Card  5 

(6F10.0) 

Required  if  IG  *  0 

1-10 

XLA 

11-20 

XLB 

Gear  spacing  parameters  a,  b,  c,  d,  in. 

21-30 

XLC 

31-40 

XLD 

41-50 

XNOD  - 

Number  of  A's  (DELTA) 

51-60 

NXOSG 

-Number  of  (o,G)»s  (ASIG,  AG) 

Card  6 

(5F10.0) 

Required  if  IG  =  0 

1-10 

XNA 

11-20 

XNB 

21-30 

XNC 

Gear  spacing  parameters  n&>  n^,  nc,  n^ 

31-40 

XND 

41-50 

PHIE  - 

Ratio  of  length  to  width  of  the  tire  print  (optional) 

Card  7 

(5F10.0) 

Required  if  IG  *  0 

1-10 

B 

Number  of  tire  approximation  points 

11-20 

XOP3  - 

Print  option.  Enter  1.0  to  print  subtotals  for  each  wheel 

21-30 

BIGX  - 

Location  of  stress  calculation  point  relative  to  gear,  in. 

31-40 

BIGY  - 

Location  of  stress  calculation  point  relative  to  gear,  in. 

41-50 

X0P6  - 

Coordinate  option.  Enter  1.0  when  using  optional  H-51  input  method 

Card  8 

(10F6.0) 

Required  if  IG=0 

1-6 

GAMMA(l) 

-  Position  of  gear  (angular),  deg. 

7-12  GAMMA (2) 


GAWA(XNOG) 


III-2 


1 


Card  9  (10F6.0)  Required  if  IG  =  0 

1-6  DELTA(l) 

7-12  DELTA (2) 

*  • 

DELTA (XNOD) 

Card  10  (8F10.0)  Repeat  XNOSG  times,  Required  if  IG  =  0 

I- 10  ASIG(I) 

Input  Pairs 

11-20  AG (I) 

Card  11  (3F10.2) 

1-10  HMIN  -  Minimum  slab  thickness  to  be  evaluated,  in. 

II- 20  IWAX  -  Maximum  slab  thickness  to  be  evaluated,  in. 

21-30  HINCR  -  Increment  by  which  the  slab  thickness  is  increased 

from  HMIN  to  f#<AX,  in. 


Card  12 

(7F10.0) 

1-10 

AMR 

Modulus  of  rupture  of  concrete,  psi 

11-20 

CUMR 

Coefficient 

of  variation  of  AMR 

21-30 

ALPHA  - 

Load  transfer  coefficient 

31-40 

CVALP  - 

Coefficient 

of  variation  of  ALPHA 

41-50 

CVE 

Coefficient 

of  variation  of  ETBAP 

51-60 

CVH 

Coefficient 

of  variation  of  slab  thickness 

61-70 

R0EMR  - 

Correlation 

coefficient  between  AMR  and  ETEMP 

III-3 


Card  13 

(3IZ) 

1-2 

NON 

Number  of  N's  (coverage 

levels)  to  be  input 

2-3 

JOPPLI  - 

Plot  option.  Enter  1  to 

plot  Reliability  vs 

log  coverages 

3-4 

JOPPLZ  - 

Plot  option.  Enter  1  to 

plot  Reliability  vs. 

.  Slab  thickness 

Card  14 

(10F8.0) 

1-8 

BN(1)  - 

Number  of  coverages 

9-16 

BN  (2) 

BN (NON) 

Notes 

1.  Card  2  should  be  omitted  if  no  composite  modulus  of  subgrade  reaction  is 
to  be  calculated. 

2.  If  NOLBP  =  2,  the  data  for  the  layer  directly  above  the  subgrade  should 
be  input  first. 

3.  In  the  output.  Layer  2  will  always  be  labeled  the  "BASE",  and  Layer  3 
will  always  be  labeled  the  "SUBBASE". 

4.  Cards  4-11  should  only  be  included  if  the  aircraft  is  user  defined  (i.e  IG=0) 

5.  The  parameters  listed  on  cards  3-10  are  defined  in  the  same  manner  as  in 


the  H-51  program.  Any  additional  information  concerning  these  parameters 
may  be  found  in  the  H-51  user’s  manual (4). 
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APPENDIX  IV 


PROGRAM  LISTING 


IV 


yv.'M 


MSI  *  MAIM  MOdl* 


MSI  COMPUTES  BIN# IMS  S1AESS  IN  tONCAEM 

COMMON  /ABIC/  A.AP<10),AM(10).AA(tQ),tSI6«13>.A6MC> 
COMMON  ten.lt  BUI. PICT, B 

COMMON  /  #  B  l  k  /  DT(31).(PI,tuU(t0),lilTAI(10) 

COMMON  /IBll/  l,LA(<(35> 
common  item  1 J < i o > 

COMMON  /CPic/  C.CAMMAMO)  ,6P*lNt(100,10) 

COMMON  /MBit/  M 1 1 0 ) 

COMMON  /me/  )CMl)0>,lLN(IPI,IT,U».tSG,IN0C  ,ABP 
COMMON  /NPt* /  N06  .NON ,NP1 ,N* ,NT ,N0#,N0S6 
COMMON  /  PBl  C  /  P.P1  ,MN"1S1  (lOI.PNlt 
common  /CBl</  AP# , A  1 1  A  p 

COMMON  /  SBl  a  /  $16*4(10). S.SA 


COMMON  /TBit/  1.1A  0 

Common  teem  m#,ma?  6 

COMMON  turn  lCtALA,l|.B.ALC,ALBlALAA.nBA,ALCAfILOAlIlA.ALA;1  fi 

1  *"U, ANA ,INB. ANC .lNB.INOC ,AMOM. ANOD , ANOS6 . «NT .TOM ,  0 

2  AOPJ.*oMj,*op*,»o*5,*oP6,jpnOo,io)  .iPiuowi,  o 

^common  /TBit/  TMN*?6  >iwPM00.10>  ,TPH*0*7>,T||AO.T|IAOM  8 

COMMON  fHiUll  CCANIJJ  _  p 

COMMON  /Bt  A  1 /  NW , N ! , C V ( 700 >  01 

COMMON  mem  UNUMMOilOI.UMAl  (10)  p 

COMMON  teen  AO , A 1,A?,AIPMA ,AMA,BN( SO) , C VBB  V ( 10) .  EVA , 1C, NON ,  0 

1  iOPPlf  ./0PPI.2  0! 

BINfNSlON  AAMS.S)  t  0 

CAT*  AA  /  -7C6Z6..U6.  ,-115728.  , -S3?t?  t .  12JU* .  16  .*77  J  . .  0 

8ijt  | S{  5::  t 

. .  s, 

CONTINUE  0 


ODOCCOSO 

00000040 

C6f£c6ip 

ooo?Iffc5 

CfODl;;7D 

POOOSDPO 

B52OD1O0 

009C„’i* 

00C.C01/O 

focr-i5P 

poooom 

OOCOuiSC 

pcoccuo 

OOOOOiMO 

55000700 

Q5Q03?iC 

00003770 

OOOOC750 

8E83Ei;8 

00000760 

00003770 

IpEt'SE 

8S3EEIT8 

OOCOOI7C 


10  CONTINUE 

ISO  A | AN ( 1 1  ■  APAOC 

AtANl?)  ■  A  N  AM  | 

AEAN(S)  ■  1 
BE  *M (4  >  ■  6 
BE  AN ( S  )  ■  0 

200  A? !#(s!^i!f N#I»^»> ic.noiop.ia.cva.etemp.amutim 
IE  EETEMP.CT.O. )E>ET(HP 
If  ((TEMP. GT. 0.160  TO  A00 
{>4000000, 

A  30  JMlMUIr  M.6T.0.)ANU»AMUTEM 
IMAMut{m.6T.0.)CO  TO  S03 
««U*0.1S 
SOO  CONTINUE 

If (nOiB*.Eb,OICO  TO  1003 

,  C  A 1 1  ACOMM 
loop  CONTINUE 
1010  11(16)7000,2000.2010 
2030  C  A 1 1  PA 33  A  D 


UODDO ) S3 

88ESeU8 

888S8U8 


^3000470 

procitjo 

ooocottf 

ppoootjc 

ooccoieo 


1010  A(  A#  (  S  ,71  Q  I  MM  I  N, 
210  f OAmatIjM0.2) 
N(1>»mmin 


NMAI.MINCA 


Oocoo; jc 

B8B8EHS 

ogcoosB* 

min* 

ccBoEIo? 

55330*10 


NOM.MOM.1 

N(N0n1«N(N0M-1)*MINCA 
60  10  273 
300  CONTINUE 

„  -  ll«l6)210f .2100,2110 


J|‘l6)210f. 

2110  AD*  a  a (16.1) 

A1  *AA (16,7) 
A2  *A  A  (  16  •  )  ) 
60  TO  ipOO 
2100  00  100  I  IN* 
C»tt  6f ON 
C a 1 1  Obit 


H6J? 
640 

6sc 

6*8 


1.N0N 


cob.jiP 


03033»49 

O0C307S9 


IV~1 


r 


w 

n 

It 

u 

*9? 

92 

93 
ft 
*3 

*a 

99 

m 

m 

& 

m 


in 

12/ 

8 

i! 

139 

Vs 

1*3 

1*4 

1*3 

14* 

1*7 

1*8 

1*9 

150 

151 


Ml 


If(*OP6.GT.0.5 I  CO  TO  fto 
umam jlh)*uhu"M ,1) 

00  101  IU2*1,N0B 
90  101  144>1,N0$G 
IMUMK  |LH>.GT.UNUN(IU7.1  JA71 
U"A*  HIM  )  «UNUN  (IUJ  ,1  J»  > 

101  CONTINUE 
60  TO  133 

UN*I()LM)«PVHTST(1) 

CONTINUE 

*0*0. 

•  1*0. 

*2*0. 

3000  CALL  *IL1»B 
CALL  S10IAL 
CO  TO  13 
9999  STOP 
IN9 
C 

c 

a'OR'lS  9ELIAB 

SUBROUTINE  REL IAB 
C 

c 


CO  TO  131 


COMMON  / ABLR /  A.  AP <  10 >  .  Ai< <  1 0 >  •  A<  (  1 07  •  A  S I C  C 10  > •  AC  ( 1 0> 
COMMON  /BBC  A /  BiCl.BICr.B 

COMMON  /OBI  A  /  0P<31>.BP*,DELT»M0).0ELTAR(10» 

COMMON  /E8L */  E.E*OOD> 

COMMON  /f  BLR /  >4110) 

COMMON  /CPI*/  C.CANNA(10>,CPBINT(100,101 
COMMON  / NbL  * /  M«10> 

COMMON  / 1 BL  K /  1 C  M(| 0  >  1 1 


iILH*|PT»IT,JLD.IS6,1 NOG  .*BP 


COMMON  /NBL*/  MOC .NON , NPT . NI ,N T , NOB, NO SG 
COMMON  /PBL*/  P.pJ.PMMTSTIIOI.PNIE 
COMMON  /RBL* /  RPD.RIERO 
COMMON  /SSL*/  SICMA  C 10 ) »S  «SN 
COMMON  /1BL* /  1.1 R 
COMMON  / NBL  * /  MR »MR2 

COMMON  /1BL* /  X*.XLA.UB.*LClXL9tRLAR.ILB*.KLCR,IL9RiILR,XL*2. 

1  XMU,*N*,*N8,*Nt ,*N6 ,*NOG .XNOH.INOO.INOSC.IHT.IOBI  , 

2  XOP2f*OP3.*0P*,*OPS,XOP6,lP(lfl0,t0)  ,*PI(*097>, 

3  RtERO.xIEROP 

COMMON  /  VSL*  /  TNN(10),WP(1OOi1O>  ,»P|  (*0971  .HERO.  T/EROP 
COMMON  /N5 1 A  1 /  RERNE5) 

COMMON  /BL* 1 /  NM.NZ.CV(7Q0>  .... 

COMMON  /UBLK /  UNUM(IO.IO) .UMAIC101 

COMMON  /MPJ/  *0. Al. A? .ALPNA .AMP.BNC 101 *C  M06M 101  §  CM* » 1C . NON , 

1  JOPPLi  .J0PPL2 

BI MENS  ION  AC].3T,*A<J,4>tBBf2,S>,BBri)3>.PRt3)>.«N(33>.9f3>./CI3» 
8  *1(10). «(3) 

BATA  RPROC/3MMST/.KNAME/4NN51A/ 

*  *  1  • 

*2*1. 

BO  1  1*1  ,3 

8(11*0. 

BO  1  J-1,3 
A(I.4>*3. 

CONTINUE 

RE  AO (5 , 1 000 )*"*. CRN*, ALPHA, C»AlPtC*E,CVM,»OEN» 

RE *9(5(70) NON, JOPPLI •J0PPL2 
EOAMAMJI?) 

EORMAT (7E10.01 


1008 

.  NRJTf (A. 10017 

1001  EORNATMmI) 


!»‘l 


8.07CO  TO  2000 


BO  w  |lm*1 .non 
AX*C*M(11M>**3./11.73/XR 
ALC  !LH)»«x**0.25 

UN*(I!lm)*A0**1*AL0C(AL(ILM>I*A;/Al(1LM) 
2001  CONTINUE 
CO  TO  210 

2000  60  10  ILM*1,N0M 

AI*E*H(Ilh)**3./11.73/IK 

AL(ILmJ*ai«»0.25 

I1*AL0C(ALMLN>> 

*2*1. /AHUM) 


IV- 2 


A(1,1)*i(1,1)*1. 

All 

A(t,3)«A( l,3)«>2 

A(t,2J*i<2,2)«ii*ii 
A(2,3)*A(2,3 >*Rl*vJ 
I(3)<M3)m1<UMI  ( luo 
A(3«5)*A{3.S)*R2*R2 
Srr*srv*o*Ai(i  lmwunav  cum) 

10  CONTINUE 

A(2,))*A(1,2) 

A(J»1>*A(1,3) 

A(3,2)*A(?.3) 

•0  Jo  1*1,3 

AA(I  .0*8(1) 

BO  20  J* 1 , 3 
20  AA(|,J)-a(i,J> 

TAWBOWMI ,1 ) 

VC1)«4. 

CALL  6JR(RA,4,3,3,4,8100tiC,V) 

A0«4A(1 ,4) 

A1*AA(2,4) 

A2* A A( 3. 4 ) 

60  TO  ISO 

100  gSMM‘12?’.}’ 

BB(I,3)*B(l) 

DO  130  J*1,2 
130  BB({.f)*»Cl.)» 

CALL  CJR(BB,S,2,2,3,S999,/C,V> 

A2«0. 

ISO  CONl 1 NOT 

V2EST*A(1,1)*A0*AQ*A(2,2)*A1*Al*A(3,3)*A2*A2«2.« 
I  2.  •40*A2*i<f,s)«2.*Al*A2*M2,3) 

VAVN*VAV*A(1,1) 

»  2AVN*TAVN*VAV 

AJ»(T2IST-T2*VN>/CSTT-T2AVN) 

A*SBBT (R2 ) 


A0*A1*A(1,2)« 


11(16)200,200.210 

210  HAIU(A,;flllt 

211  FORAATi/.SV.’uSAF  AIRCRAFT  6R0U*  INBEV :  '.12,/) 

il?g  rsi;5«*,Sta> 

1 F  OOPPL 1.LT?fPio*TO*?l 

ti  lsfcVt:jirw,€*w . 

?o*° 

S1C*U«A1(1)/H(1)/N(1) 

MCT1»(A1>A?/*L(1 ) 1/4 . /UNR* (1 ) 

FACT2*-2.*3.*FACT1 
B(*ANR /ALPHA /SIC 

CVBF*CHNR*C  NNR*C9ALP*C9ALP*F AC  T 1*  FACT  1 *(CVE*CVE*CVK*CVR)* 

I  FACT2*FACT2*CWH*CVH-2.*FACT1*R0HA*C9E*CVHR 

C9BF  *SOR  T  (C  VBF  ) 

CVOBMl)*Cy(F 

HAITI (6. 101 2)  1,M(1) .SIC.BFjCVBF.FACTIjFACT? 

1011  FORNATi/.ll.'NO.'.SR.'TNlC*  ,4l,*SI4NA  ,51,  0F',Sl,'C*8F',5V, 
1FACT1',4V.»  At T  2  '  ) 

1012  F0RNAT(14,Jf9.2,4F9.4) 

■  ■*3. 

HAITI (A, 1013) 

BO  40  J«1 ,16 

ZV*E  VP(-ll*l V/2.)/2.5066283 


AIRCRAFT  CROUP  INBEV:  ', 
A1.R2.A2. A 

A1  A2 


TI*1./(1. *0.2316419*11) 

P1*1.-ZV*T»* (5.31 9)81 53*TR*(-0.3}S5SS782* 
TR* (1.781477937*11* (>1.121233978 


C  T»* (1.78147 79 37*TI*(«1 

^BF ( ^)*^F*(1.~CVBF*I1) 

BBF(32-| )*0F«(1.*CVBF*11) 


8*11*1 .330274429)))) 


00002270 


PR  ( 32- J > 

ion-5.; 


•  1.-M 


503 


11*7 

•0  39  K»1.2 

IMH.U.  O.)  CO  TO  31 
II  (IK.LT.330.)  CO  TO  33 
1MIK.LT.L30.)  CO  TO  33 
11(001(11). CT. 0.91)  CO  TO  37 
M(|))>iO.«((tOMU)>O.S(01)/0.1) 

37  AM(I1)«T0.*«(»»MII)70.?46> 

CO  TO  31 

35  11(001(11). CT. 1.99)  CQ  TO  36 

*0(1 1>»10.«* ((001 (II >-0.6413)70.1186) 

CO  TO  38 

3*  *0(11)»JO. •*((001(11 1-0.0563) 70.2767) 

33  11  (00M1 1  >.CT,T  .21  )  CO  TO  34  _ 

*0(1 1)« 10.*« (<001(11 ) -0.7029)7 0.1 371) 

CO  TO  38 

ANCIDMO.** MOOT (11 >-0.06(3)70.31 51) 

CO  TO  38 

11(001(11). CT. 1.3)  CO  TO  32 
«O(ll)«T0.**((OOMll)-0.7653)7O.1645) 

CO  TO  39 

32  «0(i|)MO.**((OOMl  l>-0.069)70  >3382) 

39  COOT  1  OUT 
CONTINUE 

00 ]T((6,1016>  ((*»(K>l00MK)l«0(K)>fK«T .31) 

•0  62  7*1,31 
AN(71«*10C10(6N(7>> 

11(60(7).  LT.Q.QQOT)  4N(7>*0,0001 
IT (60 ( 2 ) . C  T.6 • 999 )  «*(7)*6.999 

62  COOT 10UC 


36 

31 


39 

40 


1M70TPL1.LT.DC0  TO  72 
C*ll  PL3TUM6N.PI, I) 

COOT 10UC 
COOT 10UC 
CO  TO  998 
M01TM6.1013) 

1 0*0  6  T ( *  RELI6 

TOKO  61  (2112. 4, MS.  0) 

1013  f 000  6T ( '  ILL  06TK1K  •• 

ilTUKN 
100 


TOT  3 
1014 
1013 

998 


or 


o') 


TOO, 13  110130 


SUBROUTINE  TlNlSH 


COUPON  76BLK7 
COPOOO  7BBLK7 
C0OO0N  7BAIK / 
C0OO0O  7I9LK7 
COPPOO  7TBLK7 
COPOOO  7CBLK7 
COPOOO  70BLK7 


*t *P (10), *0 (10), *K(10>. 6313(10), *3(10) 
BICK.BIGy.B 

•  P(31|^PR,0ELT4(T0>  ,»ELT*R(10) 


it\lO) 

4. 46 PP 6(10)  ,  4PR1PT (100,  1 0) 

_  H(T0) 

CONNOR  71BLK7  >C0(>0), UN, JPT, IT, III. 133, 1003,  KBP 

COPPOO  7NBLK7  003, OOP, OPT ,01 ,09 , POD ,00 S3 

COPPOO  7PBLK7  P.PI .PVPTST (10), Poll 

COPPOO  7RBLK7  RPR , R I E  R  0 

COPPOO  7SBLK7  3 ICBA ( TO ) ,3 , 3* 

COPPOO  7TBLK7  T.TR 
COPPOO  7MBLK7  N*,tf*2 

COPNON  7KBLR7  KK,X16, 118, 11C,XL»,K16R,IIBR, KICK, 11RR,XIR,XI*2, 

KPU.X0A.K08.KNC .XNR.XNOt .XOOP.XOOO .XOOSC .KOT.XOPT , 
10P?,lOP3,l0P4,XOP3,X0P6,fM10C,1O)  ,XPI (4097) , 
xriRO.K/(*op 

COPPOO  7T8LK7  TPO(TO).MP(100,10)  »TPM 4097)  , Tl 8*0 , If  E ROP 
COPPOO  7031*17  K|R0(3)  ’ 


! 


CONNOR  7BL 
UR  HE 


BLKT7  Otf.Ol .CV(700) 
.T0>_(1CP(I),1«T,101 


10  IO*P*T(tp^10x;iCi*6) 


oo  m 

“2320 

ooojiip 

0002130 
"■  S60 
*.-.370 
0002380 
2390 
'  *30 
.410 
00002(20 
2430 
244ft 

002470 

8BIU8 

00250ft 
00251 0 
0002520 


mm 


J257( 
3002580 
002390 
302600 
32610 
32620 
--32630 
000026(6 

83881*18 

00002670 
302680 

00002710 

~m 

02750 


3002790 

stef 

00002830 

0000284ft 

mhii 

02890 
1002900 


9*0 

0002968 
0002970 

0088^998 

00003000 

00003010 

SBBB18I8 


IV-4 


MMM  *  N  f 1 L N ) 

WAITE  (6,20)  HIIIH),  IK.  E,  KNU.  R2EA0,  SUTMIlt) 

20  FORNAT  (1m0.10k,22hRuRWAT  CHARACTERISTICS  /  tlK,3MH  >tF8.3, 

1  4h<1R).5k  ,IHK  »  ,  F  8*1 ,11H(LBF /  1 N  •  •  3  ) ,  5  >  ,  3  H  E  >.110.1  / 

2  11X.4HHU  »,F6.3,5X,4HAO  >,»3.2.*M<IN)  ,5x,7M»f  IT*  «,F7.3  > 

WRITE  (6,30)  KIR. ILB. UC. ILfi 

30  FDRNAT(1h0,10k  ,  4hsea«  1 

1  In  ,10k, 3nr  >,(7.2, Bn (in)  B  >  ,F 7. 2 ,8n (! n)  C  «,F7.2, 

2  BM(1N)  6  ■  ,  F7 ,2  ,4M  (  IN  )  ) 

WRITE  (6,10)  KNR ,KNB,KNC ,KNO 

40  FORRATCIh  ,10k  ,4NNR  >,15.1  ,6K,4NNB  •  ,  F  5 . 1 , 6  K  ,  4  HR  C  ■  ,  F  5 . 1 , 6X  ,4NN»  > 
I  ,  F  5  •  I  ) 


)i  i]  i  if  &  fir 


0003070 


NRITf  (6,30)  R.RMli 

50  F0A"AT(lHO.1OK,?6HCpNTACT  AREA  OF  ORE  TIRE  «  ,F 9. 2  ,8N ( S S . I N  .  ) , 

1  10K,  6N AN ( E  »  ,F5,3)  * 

WRITE  (6,60)  P 

60  FORNAT ( 1  HQ,  10K.20NINFIAT ION  PRESSURE  * , F 7 . 1 ,5N (PS  I )  > 

WRITE  (6,70)  4,6 

70  FORNATOhO,  tSi.IIWCEAR  LOAD  >,F10.1.5N(IBF)(10N  B  >,FB.1//> 

BO  500  I  NOG  ■  1.N0S 
SAHNAMNOS)  ■  SANNA( 1N06) *BPR 

WRITE  (6,90)  1  NOS .6 AN" A ( 1  NOS > , PVNT S T ( 1  NOS > .THRU  NOS) ,S16NA(lN0t> 

90  FORNATIlN  ,10i,6nSANNA(,I2«3N)  » , F 8 . 3 ,9H (BEG >  N  >,F7.2,3n  NN  *, 

1  F9.2.8M( 1N.IB.)  t 

2  IN  .101  ,17HPAWENENT  STRESS  » ,F T 1 . 3 , 5" ( PSI > / / ) 

10  UIKOPS  -  0.5)  160.  160.  110 
10  WRITE  (6.120) 

20  FORNAT  (fHO,Ul.1N!,14X,lHN,10t.8HSUBTOTAL,10K .5N6ANNA) 

BO  150  J«1 ,N06 
BO  150  1*1,  NPT 
B1STY«RIE*0  •  IP(1,J) 

B1ST|*R7ER0  *  WP(I.J) 

STRESS  •  (R7ERO*R7ERO/lOOOO,)*P*SPRINT (1 , J ) 

STRESS  *  6.*STRE$S/(HMN*HMN) 

WRITE(6.130)  DISTT.BISTK,  STRESS  .SANNAUI 
130  FORNAT  {(ION  ,2<M0.4,5K),2(F10.3,5K>>) 

160  CONTINUE 

BO  600  1 NOS  *  1,  NOS 
600  SANNA(INOS)  *  6AN"A(JNOE)*RPO 
RETURN 


•FOR. IS  BLOCK 


BLOCK  BATA 


(ONNON^BIK,/  NN.N^CWETOO) 

i  OgtOgtO^tOgtU* *0  • *0 • «Ob«0 i«Or iOo *  0  • • 0 • *0 • • 0 • • 0 • «0 • *0  • |0a i0b»0« t 


D  jvts^.30i;yg.  ir  'OoTVicrTiof  > v«o  i  «  irr.o  Tf 

BATA  (Cwlj), )»176, 325)  /OstT93. 18, 343. 27, 4J8. 41, 545.37, 609. 60,  „ 

J  655.41 ,7o’.64, 712. 96, 693. 75, 656.43, 607. 82, 552.89, 495. 29,4:7. 69, 
K  361. 96 *329. 41, 280. 86,? 36. 75, 197. 31, 162. 52, 132. 22,106. 17, B4. 03, 

L  0* *215.61,38? ,94 ,51 5,16 ,616.59 ,692.58 .747. 58.B3 8 ,45 ,821 , 26 , 

N  86l.}i.7i9. 75,7fi4. 47, 641 .42, 5T4, 95, 508.25. 4(3. 57, 382. 47,325.96, 
N  *74. SB. 228. 61, 168.05 ,152. 73, 122.37196. 58, 96. 58, 

0  0. ,246. 13, 446. 42, S06. 88, 732. 99, 829. 74. 931. 53, 964. 96, 1007. 91, 


,*77.98,3)*. 10,519 
tit. J*S?6,*7S)  10. 
5. 1590. *7. 1*50 
,43?.?8,7lS 
.137.78 


91g 


mm 

00003930 

00003960 

8B83{»8 

0000 J99n 
0000*006 


.10.319. 99,733. 3?.? 00. 77.153. *5 ,133  . *5/ 

7  5)  /0.  ,}1 J  .0? .571.73.799.85,981 .97 .1178.36, 

*30.1*. 1**5. 38, 1389. 98, 1303. *0.1196. 05, 1077. 76 
.29. 603. 9?. 303. 79. *13. 79. 336. 77. 76  7. 17. 707. 8*. 

37.78. 

*  ,  59 J . 95  , 873. 03. 1017. 11. 11 6*. 98, 178 5. 3*. 14*7.1*. ISO  7. 91, 

.1*51.0*.  1367. *5.^ 175 1.35,1 17  7.56,998. 78  .870.73,7*7.60, 

, 1*98. 00 .1*07. 9 5,1 79 3.9* ,1166. 76.10 3 7. 90 .899.9*. 7? 1.*2, 

<*1.«S,*4 ?. 78. S3*.45, 778.57.7 1S.J4, 138.4 7 .158.* 7. 

.551. 79. *<6. 81, 351. 74. 779. 59, 717. 31.15J. 76. 135. 74.  0000*056 

57,671 .03 ,86  3.0* . 106*. S3 ,17?6.*9 .1359.0*,  15  37 .61 ,1609.7*.  060 0*0*0 
1,1559.70,1*67. 67.13*9. 63. 1716. *7.1975, 35, 936. *9. 801. IS,  ODD  0*05  0 
.357. ?9,* 51. 95, 359.0*. 7 78. *3, 739. *9, 131. 65, 151. 65,  q335*063 

>6,673.71 ,869. 96, 107 3. 37, 1?59. 30, 1371. $■ ,1j*8. 36, 16? 7. *9,  66op*Q70 
7,1578. 58, 1*85  . 87,1  366.3* ,1731  . 16, 1088.06,9*6. *7, 8  01.51,  0000*086 

1  8$88:?88 

7. 1559. *7. 16 39. 9 8, 16**. 37, 1591  .69.1*98.2  7,1377.61,17*0.86  " 
7,957.30, 8 17. 36.680.71 ,559. 69, *50. 81 ,3 5*. 6*. ?71.78. 

So  .6  51 .30 ,8  78  .  ??,  1084. 15,1  ?5?.6?,1  38  7.11 ,1 56 7 .04 , 1 64 8 . 4 3 # 

r\  linn  t  >  «  (  73,  li  till  9  4  il,  *  i  •  9  «M6  «  /  nt  8  ft,  At!  «  « 


>.1605^70.1511 .13 .1388 ,64 ,1?49 . 56.110? ,84 .955*49, PI?. 71 , 
^5*. 77, **3.6*, 3*5. 57,760. 6®, 188. *4,1? 8. 7f, 178. 71, 

7, 635. 77,681.80, 1038. 81, 1758, 76,1393. 67,157*. 90, 1657. 03, 
.1608.35 ,1513. *3.1 390.74.1 255.78,1 107, 69, 95i.77.610.3l . 


.IS 

«  10  6.1.7  E..0*  MO.,  JOB  HO..  (1C.  OM  6.. /OK  ON-LIMI. 
<T  DM1NIH0N 
101  NUMBi ■ • 

1(6  KOIMI  1NTE6E.  61  MOT  MOK  1M6N  3  61*113. 

0 1 H I  MS 10M.L  *.i*T  COM1 A  t  Ml  MG  TNI  iolloming. 

(l)--.AOCf DU»£  COCt. 

3  NOLLI.IIH  CM*. ACKKS.  Ltd  ADJUST!  D. 

I2>--M*"t  0!  SUSiOUHN! 

6  MOL  L !  .  1 1M  CMA.ACT!  MS  0*  KSS 
(  ?  >  — or  »-l  IMI,  ON-LINE  mint  o.iion 
N?t3)»1  ..INI  0M-L19E  0911 
N?IS>»?  ..191  ON-LINf  ONL 1 
N?  I  3  )  *  3  ..INI  OM  —UNI  AND  ON-LINE 
[O--EA.0  NO.  O.IION 

N?  t  *  )  •  0  USE  (Ml  NO.  19  I  SC  *  > 

NJt*)  NOT  *0  SET  iSI*)«N?f()  AFTER  CALL  STATUS 


BiwB 
*‘)S8 
tm 

000*738 

030*1*1 

*76? 
0*?70 
0*780 
*  7  90 
*300 


0000*380 

0000*390 

0000**70 

«**3n 

*4*6 

**50 
**60 
0**70 
**BC 
**90 
*500 
*510 
*57C 

0000*550 


’-6 


fNTAT  COUNTED. 

*2(5)  IS  INC«l«E*l1f o  BT  1  EACH  Til 
...  ..  N2(|)  BUST  BE  INITIALIZED  BT  CALLING  AAOWAAM, 

NS— STAENENT  ND.  OA  SON!  OThEA  NUAEAIC  TO  INDICATE  STATE"ENT 
(■ 0?  HM1(N  EBN.MAS  CALICO  IN  CALLING  PBOSAAN.  NS  NUS I  BE  A 


pms 

00OO45SQ 

00064590 


MUD  POINT  1NTE6EB 
SUBBOUT INE  CBN  (NT •  N?  ,  N3> 
IS<7>.  N2C5) 


TOT  NONE  THAN  S  OIGITS. 


DIMENSION 


BATA  K2  ,  AS  /DHL 


6HP 


88 

too 

ns 


N?  (S>  ■  N21 5)  ♦  T 
IE  (N2 (A 1 1  90.9S.90 
IS ( 4 )  ■  N2 ( 4 ) 

If  CISC?)  -  2)  TOO.T  SO.TjO 
XT  *  *2 
60  TO  T90 
XT  ■  AS 

1»  (N2«S»  -  2)  200.503 .200 
^®OjNA  I  T£  (6,900)  1SC4),  ISIS),  1S(6>.  !S(7), 

..  I»  (N2(S)  -  T)  400.to0.500’ 

Mg  BETUAN910*  ,S<1>‘  M  '  *A2 Cl ) 

UJ?H 

•  TO 


J0t?465n 
.5004660 
50004670 
“304680 
304690 
304700 
'£>4710 
34720 

000^4748 
00004750 
00004760 
00004770 

88884  W8 


.1)N.  P»SE  NO.  .I3/6H  DATE  .Aft, 
PUJINC  .A6.15H  STATEMENT  NO.  ,15) 
IOAHAT  (T8h  EAAOA  IN  JOB  NO.  ,16, 


JOB  NO. 
SEBUENCE  NO. 
, A6 . 1 Sh 


«  .  .  .r.  .  BECK  TT*E 
.15,14m.  LINE  CTA 
MOCCOUBC  ,A  3 «  1  5  H  . 


t.^is.tih. 

c 

arox.is  neon 
c 

t 


PAOCEOUAE  .as ) 


6 »  T  SH •  DECK  TTPE  .AT.T2H. 


-.4810 
„  004820 
AT,  “  - 

.11 _ _  . 

SUB*000048sf 
OOOO4P6O 
30004  R  70 


ECX  TTPE  .AT. 00004850 
UNI  CTA.  ,1  500004840 
S.TSH.  SUB *00004850 


SUBAOUTINE  SEON 


COMMON  (ABLK /  A,AP(TO),AH(TO),AK(TO),AS16(10>.A((TO) 
COMMON  / B9 L X /  Bltl  .11(1,8 

COMMON  /DBLX/  D*<31),D*D,DELTAfTO)«DELTAN(10) 
common  /(BIX/  E.EAMOO) 
common  / e blx /  r C  c  to> 

COMMON  /SSL  I /  «.SAMMACT0»,SPBINT(100,T0) 

COMMON  /HBLK/  H<1Q> 

522???  il 5.L5< 


,IN06  , X  BP 


COMMON  / NBL  X /  MOS.NOH.NPT .Nl. NT, t OD.no  St 

COMMON  /NBLX /  «A,UX2  uuuu ]udu 

COMMON  / xblx /  »x,xla,*lb,xlc,*ld,*lab,*lb*,*lca,*lo*.*l*,»l*2,  B8bbso?8 

I  XMy.XNA.XNB.XNC.IND.ANOt.ANOM.ANOD.ANOSG.AllT.TOPT.  00005100 

l  *0*2 .TOPS ,*0P4 ,XOpS,xOP6,2P(i3o.13)  ,1*2(4095),  OOOO5IIO 


DECK  nooc 

4880 
48gu 
.4960 
...04910 
00004920 
4930 
4940 

4950 

'M 

0005020 
‘050J0 
05040 

£W 

005070 
5C  - 


*2( BO, 17 (BO* 

^s^/snilir000-10’  .*««*o*7».f»i.0,T,i*0* 

common  /BLXT/  nw . NI , C V ( 700 ) 

BfEAO  «SBAT(  SOAfl  ( E *H ( 1 LH) •• 3)  /  (  1 2 . • C T ,-XMU» • 2) • XX  )  )  ) 
XL  A  B  •  XL*  /  A2EA0 


*LBA  ■  TLB  /  A  2  E  A  0 
XICA  ■  XLC  /  A2EA0 
ALDA  •  1LD  /  A 2f BO 
NB  •  ( XNA * T . )  •  ILAB 
*LB«  (XNC-T.)  •  XLC A 


NX  •INT((1NA*1NB«T.) 

NT  ■1NT((1NC*KM0*T.) 

N*>.  NA/2. 

CBC2)  •  XLAA-XLBA 
*1*2  •  XL  A / 2 . 

IMS)  •  XLCA-1LDA 
21  •1NTIXNB  ♦  0.0001) 
BO  TO  I  ■  T ,NX 
E  *  ( 3)  •  1  -  T 


«  ANA* CXNB-1 . )»TLB* 
•  INC*(XN0*1,)*XLDA 

‘  ‘  0.0001) 

0.3031) 


!  5:8 : 


mm 

30005150 
30005160 
00305170 
3005180 
3005190 

itel!8 

00p05?20 


00005260 

00005270 

mm 

00005310 


IV- 7 


ER  (  I  ♦ SO)  *  HR2  -  ER(3)*XLBR  -  £R(2)  •  FLOAT  ((I-1)/N1> 
10  CONTINUE 

m2  xintjxno  ♦  0.0001) 

DO  20  4  »  1,NT 
|R(6>  «  4  -  1 

C*(4«7>  •  XLR2  -  EI(6)*XLDR  -  1*15)  •  FUtt  ((4-D/N21 
20  CONTINUE 
Nil  «  N I 
Nil  «  NT 

IF  (  ABS(EA(N1«50>)  -  1.E-07  )  10,30,40 
SO  £R(NX*53)  *  0. 

Nil  *  NX  -  1 
40  NPT  •  0 

DO  60  4  «  1,NT 
t)  SO  1  •  1,N* 

NPT  •  NP1  ♦  1 
XPZ(NPT)  x  £fl( 1*30) 

TPZ(NPT)  *  £R ( 4*  7  I 
50  CONTINUE 
60  CONTINUE 

IT  (  «8S (f»(NT«7>)  -  1.E-07  >  70,70,75 
70  ER(NT*»)  *  0. 

NT  1  «nV-1 

75  IF  (  N T 1  )  105,105,80 
Bo  DO  100  4  *  1 , N T 1 
DO  90  I  *  1.NX 
NPT  *  NPT  ♦  1 
XPZ(NPT)  s  E  R ( 1*50) 

TPICNPT)  «-ER(4*7) 

90  CONTINUE 
100  CONTINUE 

105  If  (Nil)  123,123,107 
107  DO  123  4  ■  1  ,NT 
DO  IlD  1  *  1,N*1 
NPT  *  NP1  ♦  1 
IPZ(NPT)  •  -EKI*50> 

TP7 (NPT )  «  EH  (4*7 ) 

110  CONTINUE 
120  CONTINUE 

123  IF  (Nil)  145,145,125 
125  IF  (NT1)  145,145,128 
128  DO  140  4  «  1.NT1 
00  133  1  «  1.NM 
NPT  «  NPT  ♦  1 
IPZ(NPT)  *  -|» ( 1*50) 

TPZ  ( NPT )  »  -£R (4*  7) 

130  CONTINUE 
140  CONTINUE 

145  IF  (OP(1 )*0.5)  170.170,150 

15Q  HO  1 T E  (6.160)  RZER0.*l*R,IL8R,ILCR,ILBR,HR.ILR 

160  F00N4T  ( 1  HO  ,2l,32MRZ£RO,XLAR ,  llBR ,1LC  R ,XLDR,WA,XL*  / 

1  In  ,10(12.4) 

170  IF  (DP(2)-Q.5>  200,200.100 

18Q  HRITE  (6,185)  NX  ,N T ,NP T,NX1,NT1,ILH 

185  FORNRT  ( 1H0, 35M  NX  NT  NPT  NX  1  NTl  IlM  /  1) 

_  HRITE  (6,190)  ER  ,  < XPZ ( I ) ,TPI ( 1 > , I «1 ,NP T > 

190  fORNRT(1MO,2l,37HER(1)-ER(100) ,(1PZ(1) ,TPZ( 11,1x1 ,NPT) 
..  1  <1H  ,  10112.4  )  ) 

200  CONTINUE 
RETURN 
END 


1  ,TPZ (l).lxl , NPT )  / 


OFOR.IS  OUTLNE 


SUBROUTINE  OUTLNE 


CONNON  / ABL X /  A,  AP  (  1  0  )  ,  AM  (  1 0 1  ,  A*  (  1 0>  ,  A  5  1 1  (  1 0 1 ,  AS  ( 1 0) 
CONNON  /BBLX /  Bit*  ,B  I  G  T ,B 

CONNON  / DBL X /  DP ( 3 1 )  .DPR , BELT M 1 01 ,D EL T AR ( 1 01 
CONNON  /  ( BL  *  /  E.ER(IOO) 

CONNON  /FOLK/  FG ( 1 0  > 

CONNON  / t BL* /  t,tANNA(101,(PRINT(100,l01 
CONNON  /MNLXZ  Ml 1 0 ) 

CONNON  / 1 nL X 4  ]CN(10>,llH,lPT,lT,)LD,I$t,IN06  ,XBP 


g88^5J!8 

n5340 
5350 
00005360 
0000}J70 

ass®? 

5400 

00005*10 

00005420 

DQ0OS43O 

00005440 

00005450 

00005460 

00005470 

00005480 

00005490 

OOOP55QP 

»5510 
5520 


00005560 

mite 

pshoS 

00005610 

m5620 
5630 


000Q564Q 

00005650 

00003660 

E  05470 
0568Q 
05690 

pm 

mull 

60005750 

00003760 

00DQ5Z70 


00005610 

OD005820 

00005830 


PQQ05880 

mm 

0000591 n 
00005920 
00005930 

§888$95n 

00005960 

ggm?«7§ 

00005990 

00006000 

MNtfH 

88888828 

00006050 

W&W 


Man  'MU',  i;ii;iw5\rtir.wi°'M“ 

COMMON  / IBIK /  MP6 . M I E  M  0 
COMMON  / SBl A /  S1SNA<  lO.S.SM 
COMMON  fmi/  1.11 
COMMON  /NBK/  MM  .  MB? 

COMMON  /  MBLK/  l(,IltlUI|llC.U»iIlU,U«l,UCMlH.UIllli;, 

1  MNU.MNA.XNB  .INC  ,1*0  ,»N06  ,IN0M,  «NOC>  ,  XNOSS  ,INT  .I0M  , 

{  MOP2,MOP3,mOP4,MOP5,MOP6,IP(100,19)  .MPI<4097), 

3  Mil MO. 17 (MOP 

COMMON  /me/  TNN(10>|UP(100»10>  .  TP  1 <  *097  )  ,t|  EM  0  ,T  2t  MOP 
COMMON  /H5 1 M 1 /  KEAN<5) 

COMMON  /BL  M 1 /  NM,N7.CV<700) 

COMMON  / UBLM /  UNU"«1p,10l,UM»*l10) 

S4M(PM>"<1 .<<PHIE*1.)<Pl<<5)<(S1N<Pl*PM/2.)>/PM!E 
S  *  SQMT(A/<PI<PMIE>> 

T  ■  $0MT<4<PMIC/P1 I 
SM«  S/M  I ( MO 
IM*  T/M2EMO 
1LM  «  1 

If  (NOD  .(•.  0>  SO  TO  a 
7  DELTAMdlB)  ■  DEL1A(  IID)/MIEM0 
a  DO  50  M  *  1.10 
PVH1SKK)  •  n«0 
DO  51  i  »  1  .100 
SPMIN1U.K)  «  0.0 
51  CONTINUE 
50  CONTINUE 

DO  305  I  NOS" 1 .NOS 

If  (I0P6  .SI. 0.5)  SO  TO  10 

If («BS(S*mm«(inOS>).LT. 0.00011  60  10  70 

If  (  ABSESANMACINOS)  -  9q.<MPD)  ,11.  1.0E-7  >  SO  10  30 

SC1MN  ■  SIN (GANNA (IM06))/C0S(GAMMA(1M0G>) 

MIEAOP  •  -<MP2<l)<SA<4ETAN/SOAT((GETAN>.*2 
1  •  PMIE""2>  •  D|L!AM(IlD)"COS(GANMA(INOS)n 
TIE  MOP  ■  -(TPI<11«l*>PHtE/SBMl((Sf 1AN)«*2 


Df  LIMB (I ID 1  •  SIN ( G AMP  A ( INOS) ) ) 


20 

30 

SO 


IfKNDi 


BISI/BIEMO 

B1ST/M2IBO 

•(DEUAKIlD)  ♦ 
•I1P7I1I  ♦  1B» 


-(MPIfl)  •  SR)  . 
(DELIAM(ILD)  *  » PllNIl) 
0.999999991 


1  PHIE""21 
SO  TO  40 
10  1 7  E  MOP 
NlfcMOP 
SO  TO  40 
1 1 C  MOP 
f  2  S  M  OP 
SO  TO  40 
MIEMOP 
NIIMOP 

M  "INTCS/4.0  ♦ 

ABP"  4"M 

STAMT  LOOP  ON  TIMf  POINTS 

im 

DO  250  1PT  •  1.NPT 
IIEBOOPKIPT)  ♦  SM  -  MIEMOP 
T7E»0*TPI C 1PT)  -  NIIMOP 

Ilf NO<MIENO<COS< SAMMAd NOG)) -TIER 0*S IN (GAMMA (I  NOG) ) 
NIEBO"BIEBO" SIN (GAMMA  I INOS)) «NIIMO<COS (GAMMA dNOG)) 
STAMT  LOOP  ON  TIME  PEM1METIM 

DO  200  A*T.>6P 

BN  "TLOAKM) 

OJ  • f LOAT ( A ) 

1M1.LE.J  .AND.  J.LE.M)  SO  TO  60 
1MN«1.1E.J  .AND.  i.LEt?*N)  SO  TO  79 
If (2"M«1 ,LE. J  .AND.  J.IE.3*M)  60  TO  69 
PM  ■  (4.<DM-DJ)/DM 

MI"MPI(|PT)-MTEM0P«SM<S8MT(1.-SAM(PM)<<2) 

NI"TPI(IPT)-NIEMOP-TM<SAM(PM> 

SO  TO  90 
60  PM  ■  (Dj/DM) 

MI"MP»(1PT)-MIEM0P<SM<  SO  M  T<  1  •  “ S AM( P* )  •  •  2  > 

N1-TP7 (I PT)-TZ EM0P<TR<SAM<PM ) 

SO  TO  90 

70  PM  ■  (2.<0"-D))/DM 

MI«MPI<|PT)-M/EM0P-SM<S0MT(1.-S*M(PM)<*2) 

TI"NP7(|PT)-T2EMOP<TR»SAM(PM) 


6 
6 

.  _  61  . 

000*110 

“H8 

6130 

OOogii^o 

00006700 

06730 
. .067*0 
0006230 
jB006?60 
00006270 

mim 

Si 

5006340 
$006330 
5006360 

!  006370 
006380 
006390 
006*00 
006410 
30064  20 
50006430 
00006440 
50006430 

56A$8 

00006480 

56510 

m 

?6550 
00006560 

wi?8 

Mil*0 

300661 
30066 
500066  ^  - 
50006640 
106650 
36660 
36670 
5*680 

30669 


VU 

II 


5292970 

3000671 . 

mtw 

3888^50 

M 

888S8S88 

88888(18 


IV- 9 


i 

I 


60  TO  00 

so  pa  ■ 

ii>apm  ipii-aie aop-sa« 


to 


TI*»PI(!Pt)-mAOP-TA*S*N<PA) 
tl'll»C6S<6ANNA( 1  NOG) 1 -Tl* SIN ( GANNA C 1100) » 


H)*11*S1N(6ANNA(|NOG))*T1*(OS( 6 ANN A f j  N  OG ) ) 
7JN1  •  ill  ♦  IjfAO)  /  7-0 
IT  dim  .IT.  0.0)  60  TO  110 
tAtt  T  KB  |N1  (TINT.  TNA,  HI,  K|.  NW,  (V,  )M0> 
If  (1NB.IT.7I  60  TO  10(5 
IT  <ASS(*I>  .IT.  1.0E-05 


-05)  60  TO  153 


35)  60  TO  TOO 

It  T ■ It  1*1  „ 

IMIE1.6T.7)  60  TO  K00 
HtITf  (6.199)  INS,  1 1 N T .  HI 
IBB  TORNAT  1 1  NO  ,32MDA1 K  IS  BuTSlSE  UNITS  Of  TABU.  / 
1  ,  5A ,6 HI  INI  «  ,(15.7,31 ,*MWI  * , 1 1 5 . 7) 

400  CONTINUE 
TNK  «  0.0 

100  CALL  TABINT  (TINT,  TNB.  H7ER0,  HI .  HH ,  CM,  1N6) 

If  UNO  .(1.  f)  60  TO  1)0 
IT  (A8S(tfH»0)  .IT.  i.or-c 
1E2'IE2*1 
1MIE2.GT.2)  60  TO  500 
HAITI  (6.199)  1KB  a  TINT.  HltAO 
169  tOANKT  (1M0 .S2N0ATA  IS  OUTSIDE  UNITS  Of  TABLE.  / 
T  ,5a, 6HZ INI  *,tlS.7,S(,7NH2CRO  ',115.7) 

500  CONTINUE 
TnB'O.O 
60  TO  130 
110  AUNT  •  -TINT 

CALL  TABINT  (AUNT,  TSA,  HI, 

INK  ■  -ISA 

IT  UNO.  LI.  2)  60  TO  120 
IT  (ABS(Hl)  .LT.  1.0E-05)  60  TO 
IE3'IE3»1 


TON 


|N0  >,12 


TON 


INO 


6950 
_  696C 
006970 
id006°BP 
0006990 
7PQ0 

S8#S!!<j 

Oc88?36j; 

00007070 

-*l28888?35B 

00007100 

‘"im 


NT,  NH,  CV,  1ND) 


120 


1M1E3.6T.2)  60  TO  BOO 
HAITI  (6.191)  INS.  All  NT ,  HI 

191  TOANAT  <1*0,)2MDa1a  IS  OUTSIDE  LINUS  Of  TABLE.  /  TON 
A  1  .SN.ZnAIINT  • ,f I  5, 7 , 5  A , ANH1  ',tT5.7) 

600  CONTINUE 


1HB  »  ,1  2 


m 

7240 

00pj?60 

0000727n 

vmw 


00007100 
00007)10 
307)20 


THA  *  o.O 

120  call  TABINT  (AUNT,  ISA,  HZEAO,  NT,  NH,  CV,  INS) 

TNB  «  -T  SB 

IF  (INS  .LT.  71  60  TO  t)0 
If  (ABS(HItAO)  .it.  1 .OE -05)  60  TO  T50 
IE*«IE*«1 

IMIEA.6T.2)  60  TO  700 
HA1TE  (6,192)  INS,  AUNT,  H2CA0 
192  TOANAT  I1mO,)?M6ATA  15  OUTSIDE  UNITS  Of  TABLE,  f 

1  41.5N1N6  «(12,5k, 7NAI1NT  >,E15.7,Sa,7hwT(A0  ',115.71 
700  CONTINUE 
TNB  '  0.0 
130  TANS'TNA-TNB 

IT  (BA (5 1  *  0,5)  T50.  150.  f 40 

140  HAITI  (6,95)  II.  TI,  21,  HI,  A  TE AO,  11110,  11(10,  Hit  AO, 

T  TINT,  6iNNA(|N06),  TNA.  TNB 

95  TOANAT  (fH0,6i. 2hai,13a,2HT], 131.2mm. ISA. 2HHtf  4(tt5.7)7/ 

1  5A(5MAtE«6,TQA,5MTliAj,1OA,5MTUAO.T0A,5MHZEAO/  *(f15.7)// 

6A  ,*HI1N1  ,6A,11N6ANN»(InOG).9A,3MTNA,1?A, 3mTNB  f  M(15.7I) 

6PAINTMPT  ,IN06)  ■  GPAINlUPf.lNOG)  ♦  TANS 
AllAO'Il 
TIEAO'TI 
UEAO'U 
HTEAO'Hl 
CONTINUE 

PVNTSTlINOG)  ■  PNNTSTI1N06)  ♦  6*A  1  Nl  (  I  P|  ,  I  Not  > 

IT  (lOAS  -  0.5)  250.250,230 

230  ZPUPT.lNf'G)'  COS(6ANNA(1NOG))«(IPZ(IPT)-AZEAOP)  -  S  IN  (6  ANNA  ( 1 N06  ) _ 

1  )«(YP|(lPT)-TZ|kOP>  00007520 

HP ( 1 PT  ,  I NOG 1  *  S1N(6ANNA(|N0G)>'(AP2(1PT)-AZE AOP)  ♦  COS I6ANNA ( I NOGIOfioO 7 5 30 
1  )*(TP/(|PT)-TTf AOP)  000075*0 

250  CONTINUE  00007S50 

TNN(lN0G)>A7EAQ*A?f AO/10000. O'P'PNNTST  (1N06) 
SlGNA<|N0G>'6,Q'««N(IN0G)/H(lLM)"2 
IF  (NOSO  .(B.  0)  60  TO  335 
SO  2B0  1 SG *  1  ,N0S6 


150 

197 


200 


07330 
7340 
7350 
7)60 
0007370 

tmmi 

“gglifg 

007*20 
007*30 
007**0 
. 007*50 
0007*60 
,0007*70 
00007*80 

KW 

7510 


IV-10 


•  MISS)*  *6CI$6)/<(NT*41  „„„„„ 

4N ClS6>«tZER0*SaaT(ABSCC6.0* APE !SG)*PVNTSTCIN06))/ (10000. -AS16 
1  CISC)))) 

A*ElSi>  *  C3*.1/»NCIS6>>-C10.07*HCIS6>/AZEAO>*-4 
280  CONTINUE 
290  CALL  CU*VE 
305  CONTINUE 

CALL  I1NISH 
•0  330  ] NOS* 1 .NOt 

330  UNUNCILB.INOGI-PVNTSTCINOG) 

II  <N0»  .(8.  0)  (0  TO  320 
ILO  •  IL6  ♦  1 
II  (ILO  .IE.  NOO)  SO  TO  7 
320  BE1UAN 
(NO 


E8?b 

7620 

3888^58 


aioR.is  curve 


007660 

QQ7670 

007680 

007690 

007700 

mw 

0077  JO 

BB&B 


subroutine  CURyE 


CONNON  / ABL(/  A.AP 1 1 0 > , *H t 1 0 > , A( « 1 0> . » SI S ( 1 0 > , AS (1 0) 
CONNON  /BBL(/  BIS1.B1ST.B 

CONNON  / 6BL( /  6PC3l>.9PR,6ILTAC1C>»D£lTARM0) 

CONNON  /{  BL* /  E,t»(*JO) 


1007780 

1007790 


CONNON  / C  BL*  /  E.ERII33) 

CONNON  / 1  BLN /  I6C10) 

CONNON  /6BL( /  S.SANN6M0)  ,6NR INT <  TOO ,  1 0> 

CONNON  /MBLK/  W/f/>> 

CONNON  /I  BLN  /  1CNTiQ),1LN,IP1.1T.ILD,IS6,1N06,  (BN 
CONNON  /NBL( /  NOS ,N0H,NP1 ,NI .NT .NOO.NOSG 
CONNON  17811/  P.Pl.PVNTST  MO)  ,PM1  £ 

*/  RNB.RIEAO 


CONNON  /RBI  ...... 

CONNON  /SBLN /  S1SN A ( 1 0 > , S , SR 
CONNON  / TBL( /  T.TR 


0007B60 

0007870 

000788Q 

0007890 


LVKDUK  / IBLN#  III! 

CONNON  / VBL(/  UR.HR2 

CONNON  /XBL( /  I(.(L*,ILB,ILC,(L»,IL*R,(LB*,(lCt,(L0».(L*,ILR2, 

1  (NU.INA.INB. INC ,RN6,IN0G,XN0M.XN0D.IN03G.XNT.X0P1 , 

2  IOP2,*OP3,XOP*,XOP5,*OP6,ZPC100,10>  .  IP  7  (A  09  7  )  . 

3  I7E A0.I7ER0P 

CONNON  /TBLK/  TNN C 10 > . WPC100 • 10)  , TPM *097 ) , T7 E» 0 , Tj E • OP 
CONNON  / H3 1  A 1 /  (ERNES) 

CONNON  / BL  * 1 /  NH • N I . C  ¥ ( 700 ) 

6ANNAE1N0G)  *  SANNACINOG)  •  OPR 
WRITE  (6,10)  E.  INU 

10  TORN AT  C1H1.41,  SHE  -,110.1,5*,  SMXNU  -, 1 6.3  > 

.  WRITE  16,30)  (LA|(LB|(LC,(LO 


5NXNU  I  6.3  > 


40  EO 


..Ml  «o,  VJl 

30  IORNAT (1  NO  ,101,  (nsear  / 

1  IN  .101. 3NA  -  , 1 7. 2 • BNE 1 N)  B  * , 1 7.2 , 6h < i N)  C  -.17.2, 

2  she fw)  i  -.17.2  ) 

40  IORNATEl64l2lO**46NA#l*l5tl ,6X,4NNB  « , I S .1 ,6* ,t HNC  « ,15 .1 ,6*,4HNB 

1  ^|So|  )  ^  ^  ^  ^ 

30  I g@N AT M N3?l  0* *26hC0N1 AC T  AREA  01  ONI  TIRE  • , 1 9 . 2 ,8NI SB . I N . ) , 

1  10*.  6hpm 1 E  -.13.3) 

WRITE  <6. 60)  HMtN),  SANNACINOG),  BELTAE1L9).  PVNTSTEINOG) 

60  IORNAT  MH0,4X,  JHH  *,18.3.5*.  7HS ANNA  -,78.3,51,  7H6ELTA  «,I8.3, 
1  5*.  3MN-  .*11.3) 

00  1O0  ISC-5, N0S6 

WRITE  E6.70)  ASISEISS).  ASEISS).  APEISG),  ANEIjS),  A(EISG) 

70  IORNAT  E 1  HO  *2  * ,  7N5 1  SNA  -.1  6.3.2*,  |NS  -,7  8.1,21,  JNP  -17.1,2*. 


00007910 

00007920 

8888^28 

00007950 

00007960 

8888^^8 

00007990 

Q000800D 

88882818 

00008030 

OOQQ83‘Q 


70  I3RNAT  ( f HO ,2  * ,  7N5 1  SNA  *.f 
1  3HN  -78.3.2*  ,3N(  »,*8.1) 
100  CONTINUE 


»NS  -,7  s.1,21,  3NP  *17.1,2*. 


8888il|8 

00006140 

000081S0 

sons 


RETURN 

EN6 

C 

bior.is  tabint 

SUBROUTINE  TABINT  C(  ,  V  ,1 ,N( ,NI  ,C V , INB > 
BINE  NS  ION  CVE700) 

INB  «  1 

BO  100  1-26.30 
1IEI.LT. CVE26))60  TO  200 
ITEZ.LT.C9EI))  60  TO  300 
100  CONTl-NUE 

IN6  -  5 
IS  ■  4« 


>00822n 

>008230 


irninz 

00008270 

00008280 

50006290 

00008300 

00008310 

tains 

8888*350 


IV-ll 


ik  To'JjO 

200  1ND  «  4 
IS  •  2* 

Ik  To?2S0 

500  is  :  J-, 

2S0  90  103  j  *  1,25 

IMa.lt.CVIll)  60  TO  900 
1M1.L1.CVU)>  60  TO  1000 
•00  CONTINUE 
190  •  J 
JS  •  2* 

JL  •  ?5 
60  TO  9S0 
900  JS  *  1 
J  t  -  2 
190*2 
60  TO  9S0 
TODO  JL  *  J 

JS  *  j-1 
9S0  COWTlNUl 
11  • 1 S  *2  S 
19  *JL-25 

Rl«30*  (U-1)«2S«JS 
V1$«CVt(1> 
i2«at*1 
Vtl*CVC*2  > 

K3*SO«C10-1).25«JS 
T2S»CV<* J) 

T2L*CV1K4> 

1  f  (  1  NO  .61.1)  60  TO  1100 

1010  11*<(<1U*MS>/<CVOl>-CVMS)>)*(l-CV(JS>})«MS 
T2«M(T21*T2S)/(CV(JL)'CV(JS)))«(1'CVCJS)MM2S 
111! NO .61. 1)  60  10  10S0 

1020  T*(((T2-T1)MCV<IL)-CV(|Sm*<2*CV(tS)))«Tl 
1MIN0.61.1)  60  TO  10)0 
63  TO  20]3 

1100  1MX.LT.CV(1))60  TO  1010 
1M1.6T.CV(2S>)60  TO  1040 

60  ii  ifiio 

1040  T 1  *  U<m-T1S  >MCVC  Jl  >-tVUS)>)«(X-CVlJL>>)«VlL 
12*  (<<V2l*V2S>/<<V(JL)-CV<Jsn>*(l>CVUL>))*T2L 
10S0  1M1.LT.CV<?6))  60  TO  1020 

- -  -  JO  10  1060 


10S0  1M2.L1.CVC26))  60  TO  1020 
If (2 . 6T.C V( 50) )  60  TO  1Q60 
60  TO  1020 

106Q  T  »  tl(V2-T1)/iCVtlL>-CVUS)>>Ml-CVtIL>))«T2 


loot)  t  «  cut; 

)!ss  isansi 


■ITU9N 
•  NO 


•TOt.lS  flOIIB 


SUBROUTINE  pROBRB 


0000858 

0000851 


n*20 

430 

OOOO644O 

QQQQBiSO 


00008460 

00008470 

P?tt8 

00008500 

00008510 

00008520 

m\\h 

OOOOB550 

pi*»s 

00008580 


00008600 

0000661Q 

8?88Stij 

00006640 


00008660 

0000869n 

00008700 

0QQQ871O 

Q0Q08720 
000087 Jn 
00008740 
00008750 
00058760 

n 00008770 
00008780 
000067Vn 
00008800 
00008810 
820 
8  JQ 
00000840 


CONNON  /4BLK /  », 88110) ,RM 1 101 , A* 1 10> , * S 1 6  *  10 > , AG110) 

COMMON  /BBLK /  B16I.816V.8 

COMMON  /D9L1  /  OP151>.6P*,6ELT»110>,OELTAR110> 

COMMON  /(  BL  K  /  l  ,(  8  11  00  ) 

COMMON  /TBL*/  T611D) 

COMMON  / 6BL*  /  6.6ANMA1 10) t6PR!NT1100» 10) 

COMMON  /MBLK/  NlfO) 

COMMON  /  18L  K /  1CM110)»ILH,IMT,1T,ILB»IS6 • 1 HOC  ,187 

COMMON  /NBLK/  NOG.NON.NPT .NI.NT.NOO.NOSG 

COMMON  /PBL * /  P.PI ,PV"TS1110l.MNIE 

COMMON  /IBL* /  APD.RlllO 

COMMON  /SPIK/  SI6M8 1 10 ) .$ .SR 

COMMON  / TBL* /  T.T* 

COMMON  /ttPLK/  NR.VR2 

COMMON  / 1 9L* /  1*  ,1L4,ILB,11C |ILB,IL4R,1LBR.ILCR,1L0R1ILR,(LR2. 

1  lHU.lNA.XNB.INC.XNO.INOG .(NON , KNOB , 1N0S6 ,*N T , KOP T . 


05008860 

00008870 

Sim 

s:?( 

md 

00008950 

00008960 

sssslils 

00008990 


00009010 

0000*020 


000090 jo 
00009040 

mim 

OQODvOTO 

oooo9ofcr 

8888W8 

00009110 


»  *OP2.SOP5|SOP4,IOP5«XOPS,7P<100.10)  ,*PI<4097), 

COPPON  / 181* /  1PN < id) .9P< 1 00,10)  ,1PI<4097),1IC*O,1Z<ROP 
COPPON  /NS  1 A 1 /  *ERN<5) 

tS55o3  /PP 5/'  AO J Al • A2  I  ALPHA , APR ,BN< 1 0) #CVBOf<10)«CVK,lS» 

l  40PPtl  .iOPPL? 

DATA  KNAPC/6HPR0BR0/ 

KE  AN  1 2)  *  *P»»l 
PI  »  1-1*15927 
OPR*  llO.O/Pl 


265  I?  IsNOSS?  273,  i£s*,iZS 
270  (>U  IRP  <08,  «{»R,  2*01 


00  10  6990 
275  IE  <RNOSG  -  10.)  265.  265,  280 

l#0  ^Vo'S^0’’  *b6> 


0*1A  KNAPE/6HPR0BRD/ 

KERP121  *  *PA"E 
PI  *  5-1*15927 
OPR*  llO.O/Pl 

AEAD<5{io!eND*999)<ICN<I),I*1,10> 

10  P0RPATI1DA6) 

0P(1  1*0. 

20  f 0ANA7<lHl°)l9X,15HP  ?  9  *  m  1  “  P  U  T  #  *  1  *  * 

JO  EORPAT  <1  HO .  11*.  10 A6  1 
15  iCAO<5,*0,END«99$)G,P.A,*N06 

17  EOPTINUE  „  „„ 

40  EOANATUf  10.01 
211  if  ( SNOG  1  220,  220,  225 
220  CALL  CRN  <05,  KfRN,  2201 

225  If  ESNOG^  io.)  250.  250,  215 
215  CALL  CRN  <0*.  KERN,  2151 

250  If  EP.LT.p.OOOl .ANO.A.lT .0.0001 1  SO  ro  240 

45  fORPAI  (!h&.5*,11HK<LB/CU  1N),S*,5H6<LB),8*,6HP<PS1)»6S,8HA<S0 
1,41, 12HN0.  SAPPAC 1 1  1 

50  fORPAl'llp  ,5s*^7.l,5s,flo.1,5S,f7.1,5s,f9.1,5S,f*.1) 

SO  10  245 

240  CALL  E RN(05,  KERN,  240) 

SO  TO  6990 

2*1  MAD*  <5  ,41 SLA,  SL8 «  ICC,  *L0,  RNOD,  SN0S6 

250  ^’We^'k^:  »0» 

255  If  <SN009-  10.1  265,  265,  260 
260  CALL  CRN  <0 7,  KERN,  2601* 


KinS 

00009150 

00009160 

miw 

00009190 

00009200 

«9210 
9770 
9  i  50 
00009240 

K9250 
9760 
9?7n 
00009780 

sm&B 

00009510 


00009540 

00009550 


ggggstfg 

*  >0009470 

!888Si*8 

{{885118 

00009470 

Kill 

500 

00009510 


50009550 

S8888U8 


!!5  mm  . . .  <n«>  •«*>■ 

1N1.22H  NO.  OEtTAll)  NO.IS.S)  1 
NOS  «INT<SN06*  0.0001) 

NOD*IN1<SNOD*O.OOpl) 

65  rOMAT  (IN  #5l  , *7.2 *3(6* f  *  7*2 > #7*  •  *4 •  • 

4  BCAD<5!45)SNa1kNB,«NC,SND,PH1I 

Iff<PNII. 11. 0.0001)  PHIC  ■  *»/3* 

70  f ORPAT*< Im815S.5PSNA |9R,5HSNB,Ix,5AsNC ,9S ,SHSND.9S,6PPNI<E>7 

,J?.«4!j;44s]^‘l;A«:'5.,G.,  BIST.  S0P6 

75  fDANArSlNolas  J1h8,6*.9HPr!n1  9P!l|i,j59|9  "S^fl'l")6  f*6** 
IBmcOOR  0P1  /1h  ,ss  ,f 7. 1  ,5 S ,f 4. 1 ,2 «5S,f 7.2 ) ,7S,f 4 .1  > 

42  roRPAinofA.oi 

READ  <5,421  <SAPPA  < 1 1 ,  1*1, NOS) 

•PRITE  <S,§5>  fSAPNA(l),  1*1,996). 

85  fORPAI  <1hQ,27S,14H6ApAa<I )  IDEG17 
1 1N  .21, 607.2,41)  > 

If  <N0D.E9.01  60  TO  97 


B<  1 N  )  SNAIL  C  < I N  )  SPALL 


vm 

00009650 

Hums 


00009710 

00009720 

00009750 

» 00009740 
00009750 
00009760 

9760 

9790 

OODD9SOO 

00QQ9810 


00009950 

300 ?9 640 


IV-13 


■  I  AD  <5.4  2>  OUIMI),  1*1,M0D> 

Ml|1(  <6.90)  (DELTAd),  1*1  .NOD  ) 

Do  fOANAT  <1H0,30A,BNDELTA(1 >  / 
iin  (2i.6<f7.2,4i)  > 

97  If  (NOSt  .10.  0)  «0  TO  99 

MAD  (5.92)  fASIt(I).  At  f  1  )  ,  1*1  ,N0S6) 
92  f 0ANA1 I ft  9 10. 0> 

If  A  ]  T  £  (6.93)  (  ASlt(l)  .AGO).  1*1  .NOSG  > 
95  fOANAT  MH0.6X ,7HAS16( I) ,9x,5HA6< I)  / 

99  AO  300  1  *  1.N0G 

NANNIE  1  #  •  6ANN4(I)*AAD 
300  CONTINUE 

If  (DA(1|-p.5)  340.  340.  325 

H3  nss.i'i  ifi.sr''- 

HAITI  (6.333)  (D A ( I )  .  1*1,31) 

333  FOANAT  ( 1  HO .  5HDUNAS,  31F4.1) 

340  *NT*l|JjA*IHBJ*NtJ*HJ  ^  4Q3>  4Q[) 

400  (ALL  (AN  (10,  At AN ,  400) 

60  TO  >990 

410  If  (A)  490.  500,  490 
490  A  «  t/(4< INT) 

GO  TO  StO 

5?8  G  *  A*APi!?T> 

°  60 *T 0  S4q 
>990  CONTINUE 
GO  TO  15 
999  CALL  (All 
340  A( THAN 
(ND 


>f OA , I S  ACONA 

SUBAOUTINE  ACONA 


2(3A,f10.1)/) 


CONNON  /ABLA/  Ik, ALA, HB,XLC,ALD,1LAA, (LBA, 1LCA. ALDA, ALA, ILA2 
1  ANU, ANA, ANB. ANC,AND,ANnG,AKOH. A NO  6. INOiG ,*NT  .10 

|  *OA2,XOA3jXOA4,AOAS,XOP6,IA(100,1D>  ,XA2 (409>) , 

* CONNON  /NAJ /  A&lAl!A^ALAHA,ANA,BN(lO),CVDDf ( 1 0 ) , C VK , I 6 , NON , 
1  JPAALI . JOAAL2 

1004  FOANAT (1  Hi ,161 ,  'AULT I-LATI AED  SYSTEN  ANALYSIS'//) 

TOO^fOANATjlAA^NODULUs'OF  SUBGAADE  AE  AC  It  ON  •  ',f7.2,“  ACI',5A, 

IIP  F0ANA1(14a*^1YA*:  HE  LL-GAADE  D  CAUSHED  NATEAIAlj'/) 

1009  fOANAT (14aJ'1Y At:  NATUAAL  SAND  AND  6AAVEL  (A)<j>'/> 

1011  fOANATC14l,'TYAE:  $1 A9 I L 1 1 ID '/ )  ... 

1012  fOANAKlAA, 'HEIGHT  -  ,75.2,  INCH!  S '.BA  . 'C  Y  *'f  7.4/1 

1013  f 0ANAT(14A,'N0DULUS  *  *«f 10.2,  ASI  ,5X,  CV  *  ',f7.4 //) 

1014  fOANATClH0.2A,  CONpnjITE :') 

1015  fOANAT(3A,*LAY( A  J*l 

it:  6  FOANAT (51,  'Subbase  :  ) 

7Q1  fQANAI (12.4F10.4)  ....  .... 


CONNON  /NAJ / 

)4  fOANAT(1Hl.1« 
15  fDAN*T(4A,*SL 


701  r0ANAT(}2.4f10.4) 

431  f 0ANA1(1h3,2OA ,  A  NAXINUN  DESIGN  I 
HAITI (6,10645 
HA  1  T( (6, lpp3 ) 

HAlTE(6|lp06)AK,CVA 

AE  AD (^ ,701 JnA^OAT ,H0 ,( VHB ,EB,(VEB 
If  (NA10AT.E9.DG0  TO  101 
1MNATOA1.C0.2)GD  10  107 
If (NAT0A1. (9.3)60  10  103 
V«ALQ61Q(IX) 

12* (0.08901 2 52* ( 2, •( “0.0004475 194 )*hb) *(-0.07901 4>>*T) ) 
Y3»(HB**2.)»(Y2**2.>*(},}6l8981)*(CYHB**2.) 

VC ><0.964 31 74) *<-5.256 1922«V I *(-6.01 2601 *H>) 

Y5»(14**2. >• (5,3018981 )><CYK*>2.) 

CV«C*Y3*  T5 

Y?5-??25 I ?>!♦ (2.21 9737* Y)* (-0.29495 22 *(Y»*2.)) 

I  ♦ <0. 08901 2 52 *HB)*< -0.0 0044251 94* (HO* *2.)) 


A  NAXINUN  DESIGN  I  Of  500  AC  I  IS  USED') 


00009920 

000099)0 

»99<n 

9956 

9960 

PQQQ997Q 

00309980 

00009990 


66515626 

00010030 

00010040 

OQQl807p 
00010080 
DDp1p090 
00013100 
000 1 01 10 
PPP'Ol 20 
00010130 
00010140 
000 1 Ql 50 

888i6i70 

000 lQl 80 
000 101 9n 
00010200 
00010210 
10220 
10230 

?8?58 
10260 

\m 

?8IS8 

m 

10  340 
10350 

?8?*8 
10380 
10390 

18798 
18718 

10440 

ms§ 

1Q470 
104BP 
10490 

!  » 
MX 

10550 

18578 

\m 
10600 

MX 

10630 


2m){9nf)4B8lll>iV> 


1J1 

it! 


11 3* 


SO  TO  19} 

102  i*Aiosn<iK> 

22«<0.0S55*373*<2.* (-0.00021 7*61 9) • MB) *<-0, 03050173*7) > 
IS*<Hfl**2.)«<f2**2.>*<5.30189B1)* ICVH8.*2 .1 
I* *0.9829852 *<-0.261  7702* 7  > ♦ <-0 .01 32*6 *hb» 


15*<14»*2,>*<5.3108961>*<c»«**2.> 
CVKC  *73*15 


l¥*-f.2963B4«<2. 263407*2)  *<-0.391 3741*0* 

1*<0.08 334373* MB l*<-0. 000 2S74619*<NB**2.)) 

2* <”0.030 50 173* MB* 2) 

KK*10.**2l 
SO  TO  195 
103  M»AlO610<*t> 

M2 >0.054 4 761 *<”0.001 6948 *NB)«( ”0.44090961 -2* U) 

,M3*?mI**J^*I<mI**5!  )  *<3.301 8981  )*<CVM8**7.) 

M* *<0.725474 92 -6 )*<-0. 38 7*S86<“'2*IB>* (-0. 4601 6532 -7*m> 

tM5«?2l*2fli*<Mi*5?!)><5.30l89S1 >• <t»28**2.» 

Mft«0. 4465 1 12 *<”0. 19148462  - 2*MB >♦ <-0.1 9984732-7*20) 

M2* <M6* *2.1* <5. 3018981 »*<CWK**2.» 


<7*  *2 • )) 


CV(C*M3*M3*W7 

CVK-S0BT<CVKC> 


I 


Ml >-Q.1 578667* <1. 0281 3*M I* <0.0544761 *MB>«  <-0 .84738522-3* < MB* *2.) > 
♦ <0. *2547492-6*28) *<-0.19372932 -12* <28. *2. )> 

♦<-0.44090962 -2*MB*M)« <-0.4601 65 32 -7* 2 B*M» 

♦<-3. 24656382-1* - 


9.24656382-fl*2B*NB) 
8K*1Q.**M1 
SO  TO  195 
195  CONTINUE 

IMNOLBP. <0.1)60  10 
ir<NOLB«>. 20.2)60  TO 

201  MB !TC<6.10p7) 
IIOUT0M.t9.1)60  TO 
IT  < NO  I  OPT .29.7)60  10 
IT <M»TOBT. 20.3)60  TO 

203  MBIT2<6.1008) 

60  TO  206 

H".S'l5i00” 

2 05  MBIT! (6.1011) 

.  60  TO  206 

206  CONTINUE 

MB  I T2  <6.1 01 7)  MB.  C  VMB 
If <M6T0PT.LT .2.3)60  TO  207 
_  MB1TC<6.1013)2B.C«2B 

207  CONTINUC 
MBJT2 (6.101* ) 

MB 1T2  <6. 101 5 ) 

MR112<6.1 006 ) IK .C VK 

(f <8K.l!.500.>60  TO  400 

■6*300. 

MB 1TE  <6.401) 

400  irmu 

202  |f<l. [0.1)60  TO  300 
...  If <1 . IB. 2)60  TO  301 
300  MBITC(6. 1016) 

If (MATOPT. [0.1)60  TO  302 
If <MATOPT.I«*2 >60  TO  303 
If <MAT0PT.2«.3)60  TO  304 

302  MBITE<6.T008) 

303  MB1T2<6*1009) 

60  TO  305 

304  MBIT2<6,I011> 

69  TO  3)5 

305  CONTlNUl  ... 

•  MB  I T2  <6 . 1012 )MB. C9MS 
If (MATOAT.LT. 2. 5) EO  TO  306 
M«1T2<6, l0l3)fB. CVEB 

306  CONT1NU2 
60  TO  700 


mi 

110730 

mitt 

1U7B0 
..  J10790 

J0010800 

1Q1Q810 
Ju1y82Q 
3010830 
00010840 
07010850 
00010860 
00010870 

888i8§§ 

0001090 

08818978 

mn 


OOOln’SO 

00010960 

«970 
980 
990 
OOO11OOO 

pita 

50011030 

00l i0*0 

0011050 

0011060 

001167b 

wm 

moo 

-11110 

H1158 

turn 

OooiiiSn 

>11220 
30011230 
300112*0 

22011?SQ 

Bo8iif$o 

30011280 
5001129n 
501 1300 
01131P 
,.011320 
00011330 
000113*0 
00011350 
00011360 
00011370 
00011 J8C 
00011390 


8 


IV-15 


JUS 

m 


\\ 


!** 

1144 

11*7 

11*8 

11*9 

1150 

iili 

1155 

Vl!67 

1158 

1159 

1160 
1141 
1162 

mi 

1165 

1166 

m; 

1169 

w? 

m 

v,n 

im 

1178 

1179 

1180 
1181 
1182 
1181 

jjh 

it* 

1188 

im 

191 

ji 

.19* 

1195 

US* 

1198 

1199 

\m 
m 


n 

1! 


I 

i«f 


1  tl  <N*10&1°?0  .1  160  10  507 
1MNATOP1.E0.2)GO  10  3r« 

If  (NAtOP1.E9.3JGO  to  jfl? 
507  UR!TE(6.1u08) 

60  10  llO 

308  6*17116* 1009  ) 

60  TO  113 

309  MR1TE (6,1011) 

60  10  110 

310  CONTINUE 

BR1TE<6, 13121MB. CVM9 
17<NAT0PT. IT. 2.5)60  TO  511 
«*)TE(6,1013)|8,CNEB 

311  CONTINUE 
MR1TE16, 101*1 
UNITE <6,1015 > 

U*jTE(6, 1006)11, (VI 
17(11. L(. 500)60  TO  401 
11*503. 

U* 1 TE (6 ,401 > 

*03  CONTINUE 

?99  CONTINUE 
00  CONTINUE 
>E  TURN 
END 

70*  ,1$  S 1  OR  PL 


Subroutine  siorpl 


common  / »bl* / 

COMMON  7BBL1/ 
COMMON  /DBL* / 
COMMON  / E  BL 1 / 
COMMON  /7BL1 / 
COMMON  /6BLK/ 
COMMON  /MBL1/ 


A,AP(10),*M<10),Al(lQ),ASf6(10),A6(10> 

BI61.BICt,B 

BM(11).D»*«,0EE1*(10>.DELT*«{10) 

IrtWl”* 


,1BP 


111 

1215 


6.6*MM*( 10) ,6**1  NT  1 100, 10) 

......  M(io> 

COMMON  /IBll /  1  C  7*  *  I  0), iim.ipt.it  ,ILD,I  sg , i nog 
COMMON  /NBL1 /  NO 6 .NON ,NPT . Ml ,N T ,NOD,NOSG 
COMMON  /PBL1/  P,PI,P*PTST( 10),PNIE 
COMMON  /RBL1/  RPD  «R I E»0 
COMMON  /SML1/  SI6M* 1 10 ) ,5 , S* 

COMMON  / 1BL1 /  1,1* 

COMMON  /HBL1/  U* , 6*2 

COMMON  / I BL 1 /  l*,*l*,lLB1lVC.*tO,*L»».*LB*,lLC*,*lD»1*LR,*l»2. 

1  IMU, IN A, INS, INC ,>ND,1NOC,1NOM.1NOO,1NOSG.1NT.10PT  , 

j  *OPj.*OP3,l0P*,XOP5,*OP6,7Pll00,1O)  ,1P7 (*097) , 

COMMON  / VBL1 7  TNN* ifi > *  UP < 1 00 ,10 >  ,TPI  (*097),T7E*0 , T{ E  ROP 
COMMON  /N5 1*1 7  KE*N<5) 

COMMON  /Bill/  NU,NI,CU(7Q0)  .... 

COMMON  /UBL1 /  UNUM (lO.IOl.UNAlOO)  . ..  , 
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Chapter  1 


INTRODUCTION 

General 

The  design  of  military  rigid  (plain  and  reinforced)  pavements  is  currently 
based  upon  the  classical  Westergaard  free  edge  stress  slab  theory.  It  is  this 
method  that  has  historically  been  the  basis  for  design  methods  found  in  present 
US  Army  Technical  Manuals  and  U.S.  Air  Force  Manuals.  In  addition  to  this  ap¬ 
proach,  recent  research  conducted  by  Parker  et  al.  (1)  at  the  USACE  Waterways 
Experiment  Station  has  demonstrated  that  an  equivalent,  if  not  superior,  de¬ 
sign  methodology  based  upon  multi-layer  linear  elastic  theory  is  also  another 
feasible  design  approach  for  rigid  pavements. 

At  present,  both  of  these  design  methods  are  deterministic  in  that  a 
unique  pavement  system  is  designed  for  the  specific  set  of  input  variables 
necessary  to  solve  the  problem.  In  a  deterministic  model,  all  of  these  in¬ 
put  variables  are  also  unique  in  terms  of  their  input  magnitude. 

Study  Objective 

The  objective  of  the  research  study  reported  in  this  volume  was  to  include 
the  design  parameter  variability  into  the  USACE  design  procedure  based  upon 
multi-layer  elastic  theory.  This  probabilistic  analysis  is  based  upon  the 
design  procedure  developed  by  Parker  et  al.  (1).  The  probabilistic  method¬ 
ology  for  the  Westergaard  based  design  method  is  presented  in  Volume  III  of 
this  overall  research  study.  The  mathematical  and  analytical  development 
of  both  probabilistic  approaches  used  in  this  Volume  and  Volume  III  (Wester¬ 
gaard)  are  presented  in  Volume  II  of  the  research  study. 
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Report  Organization 

The  method  of  development  used  to  establish  the  probabilistic  multi-layer 
elastic  theory  design  approach  necessitated  the  investigation  of  determining 
an  elastic  theory  based  "Composite  Modulus"  for  a  multi-layered  foundation  sys¬ 
tem  (i.e.  a  subbase  layer  over  an  existing  subgrade  soil).  Because,  this  in 
itself,  required  a  significant  effort  and  investigation,  this  volume  is  sub¬ 
divided  into  two  major  chapters.  They  are: 

Chapter  2:  Composite  Modulus 

Chapter  3:  Probabilistic  Analysis  of  PCC  Airfield  Design 

In  essence.  Chapter  2  presents  the  solution  to  the  "composite  modulus"  pro¬ 
blem,  which  is  an  integral  part  of  the  overall  probabilistic  analysis  presented 
in  Chapter  3.  The  "composite  modulus"  has  been  evaluated  from  the  BISAR  multi¬ 
layered  elastic  theory  computer  code.  As  will  be  noted,  this  parameter  is  ex¬ 
pressed  not  only  in  terms  of  the  layer  thickness  and  material  properties  but  the 
specific  loading  conditions  (aircraft)  as  well. 

Chapter  3  dealing  with  the  probabilistic  analysis  presents  the  development 
of  the  regression  equation  used  for  maximum  stress  computations  used  to  evaluate 
the  variability  of  the  pavement  design.  Two  solution  approaches  have  been  de¬ 
veloped:  (a)  the  approximate  closed  -  form  using  first  order  Taylor  series  ex¬ 
pansion  and  (b)  Monte  Carlo  simulation.  A  computer  program  has  been  developed 
to  solve  this  probabilistic  approach  and  run  example  problems,  a  User's  Guide 
and  program  listing  are  presented  in  Appendices. 
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Chapter  2 


COMPOSITE  MODULUS 

Introduction 

The  concept  of  a  composite  modulus  {of  elasticity)  was  introduced  to  take 
into  account  the  combined  effect  of  a  layered  (subbase/subgrade)  system  under¬ 
neath  the  rigid  pavement  on  its  performance.  This  parameter  is  commonly  used 
in  pavement  design  and  evaluation  and  several  methods  were  developed  for  com¬ 
puting  the  composite  modulus  as  a  function  of  layer  thicknesses  and  moduli. 

It  should  be  noted  that  the  composite  modulus,  like  the  Equivalent  Single 
Wheel  Load,  is  defined  as  an  equivalent  modulus  that  will  ultimately  lead  to 
the  same  correct  response  as  with  the  original  system  of  layers.  Therefore, 
like  the  ESWL  parameter,  the  composite  modulus  may  be  expected  to  depend  upon: 
(1)  the  pavement  response  (maximum  stress,  strain  or  deflection)  chosen  for 
equivalency  computations;  (2)  the  load  configuration  and  (3)  the  pavement  geo¬ 
metry  (layer  thickness  and  moduli).  It  is  stressed  that  all  variables  should 
be  initially  evaluated.  Deleting  anyone  variable  from  the  relation  can  be 
made  only  if  its  effect  is  found  to  be  negligible.  This  chapter  presents  the 
evaluation  of  the  composite  modulus  of  elasticity  for  the  case  of  a  layered 
system.  It  is  subdivided  into  the  following  sections: 

(1)  A  brief  description  and  discussion  of  the  method  for  deriving  the 
composite  subgrade  modulus  of  reaction  (Westergaard  model)  in  the  current  rigid 
pavement  design  system; 

(2)  A  study  of  the  effect  of  the  composite  modulus  on  rigid  pavement  per¬ 
formance  in  order  to  assess  the  degree  of  accuracy  needed  in  the  composite 


modulus  evaluation; 


(3)  A  comprehensive  study  of  the  effect  of  all  design  variables  on  the 
composite  modulus  of  elasticity.  All  variables  of  the  pavement  geometry  are 
included  in  a  regression  equation  relating  the  composite  modulus  for  one  single 
wheel  load  to  the  pavement  layer  thicknesses  and  moduli.  The  load  configuration 
effect  is  included  in  a  separate  equation  relating  the  composite  modulus  for 
the  given  gear  configuration  to  the  composite  modulus  for  one  single  wheel 
load. 

Although,  the  composite  modulus  is  not  mandatory  for  the  elastic  layered 
systems  (it  is  possible  to  compute  the  stress  for  the  original  system) ,  it  is 
of  practical  and  economic  interest  to  use  the  composite  modulus  in  order  to 
assess  the  value  of  additional  layers  and  to  simplify  the  design  procedure. 

Current  U.S.  Corp  of  Engineers  Method  for  Determining  the  Composite 

Subgrade  Modulus  of  Reaction  (kc  Value) 

In  the  current  U.S.A.C.E.  method  (2)  separate  diagrams  were  developed 
for  different  materials.  Figures  2.1  and  2.2  relate  the  composite  subgrade 
modulus  of  reaction  (denoted  kc)  for  specific  unbound  granular  materials 
(well-graded  crushed  material  and  natural  sand  and  gravel  (PI  <  8) ,  respec¬ 
tively)  to  the  base  layer  thickness  and  subgrade  modulus  of  reaction.  In 
Figure  2.3,  the  relationship  is  given  for  a  broad  range  of  subbase  materials 
which  are  characterized  through  their  modulus  of  elasticity.  As  ex¬ 
plained  in  Volume  III,  the  partially  bonded  rigid  overlay  pavement  design 
equation  for  h^c  could  also  have  been  used. 

Results  of  Figures  2.1  and  2.2  are  based  on  field  plate  tests, 
while  those  of  Figure  2.3  are  based  on  computations  using  the  elastic 
layered  theory  and  correction  procedures  (3).  In  both  diagrams,  the 
equivalency  criterion  is  deflection  or  stiffness,  i.e.,  the  composite 
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Figure  2.1  COMPOSITE  SUBGRADE  MODULUS  OF  REACTION 
FOR  WELL-GRADED  CRUSHED  MATERIAL 
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Figure  2.2  COMPOSITE  SUBGRADE  MODULUS  OF  REACTION 
FOR  NATURAL  SAND  AND  GRAVEL 
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COMPOSITE  MOOULUS  OP  SOIL  OF  FOUNDATION  BENEATH 
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2.3  COMPOSITE  SUBGRADE  MODULUS  OF  REACTION 


FOR  DIFFERENT  BASE- SUBBASE  MATERIALS 


modulus  of  subgrade  corresponds  to  equal  stiffnesses  of  the  equivalent 
subgrade  and  that  of  the  original  base- subgrade  system.  Therefore,  the 
composite  modulus  computation  scheme  does  not  include  the  effect  of  slab 
thickness  and  modulus,  neither  that  of  the  load  configuration. 

It  should  be  remembered  that  the  conventional  military  rigid  pavement  de¬ 
sign  method  is  based  on  the  Westergaard  theory  where  the  subgrade  is  represented 
by  either  a  dense  liquid  or  a  spring  (i.e.  two  layer  slab  system).  There  is  no 
possibility  of  taking  care  of  the  multilayer  system  underneath  the  pavement  slab. 
Furthermore,  the  design  method  is  based  on  limiting  the  maximum  tensile  stress 
at  the  bottom  of  the  slab  edge.  It  is  not  obvious  whether  the  composite  subgrade 
modulus  evaluated  using  the  above  procedure  will  also  result  in  equal  maximum 
tensile  stresses.  In  order  to  determine  the  composite  subgrade  modulus  for 
equal  maximum  stresses,  it  would  be  necessary  to  conduct  field  tests  on  slabs, 
with  stress  or  deflection  bowl  measurements.  Such  an  approach  would  not  be 
practical.  However,  with  elastic  multi-layered  theory  (unlike  the  Westergaard 
theory),  it  is  possible  to  directly  determine  an  equivalent  or  composite  mod¬ 
ulus,  with  any  equivalency  criterion  parameter.  In  the  ensuing  analysis,  the 
maximum  tensile  stress  has  been  used  as  the  equivalency  criterion  for  obvious 
reasons . 

Study  of  the  effect  of  Composite  Modulus  on  Rigid  Pavement  Performance 

An  alternative  design  procedure,  using  multi-layered  elastic  systems  and 
center  loading,  to  the  existing  one  using  Westergaard  theory  and  edge  loading 
has  been  developed  by  the  U.S.A.C.E.  (1_).  The  salient  features  of  this  design 
analysis  include: 
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(1)  Computation  of  the  maximum  tensile  stress  (o)  at  the  bottom  of  the 
concrete  slab,  using  elastic  layered  system  theory.  In  these  computations, 
it  was  assumed  that  (a)  the  interface  between  the  concrete  slab  and  the  base 
or  subgrade  layer  is  frictionless  and  (b)  the  subgrade  depth  is  finite  and 
equals  20  ft;  (c)  a  rigid  (stiff)  semi-infinite  layer  lies  below  the  subgrade. 

(2)  Computation  of  the  design  factor  (DF)  is  defined  as: 

DF  =  MR/o  (2.1) 

where 

MR  =  modulus  of  rupture  of  concrete  determined  at  age  of 
90  days. 

(3)  Computation  of  the  allowable  number  of  coverages  (N)  is  obtained 
from  the  following  relationships: 

DF  =  0.58901  +  0.35486  log1Q  (N)  (2.2) 

or 

fn  N  =  6.4887  MR/o  -  3.82192  (2.3) 

The  data  and  the  correlation  of  DF  and  N  are  shown  in  Figure  2.4. 

The  above  design  method  will  be  used  to  study  the  effect  of  overesti¬ 
mating  or  underestimating  the  composite  modulus  on  pavement  performance. 

This  can  be  expressed  in  terms  of  variation  of  N  (the  number  of  coverages) 
as  a  function  of  the  variation  of  Ecofflp  (the  composite  modulus) .  Substitut¬ 
ing  the  function  of  o  (a  function  of  the  slab  thickness  and  modulus,  the  com¬ 
posite  or  subgrade  modulus  and  the  loading  conditions)  into  Equation  2.3 
leads  to  the  relationship  between  N  and  Ecomp* 

Expressing  o  in  its  general  form  as: 

o  -  o( hj  .  Ej/E^^,  LOAD)  (2.4) 
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Figure  2.4  PERFORMANCE  CRITERIA  BASED  ON  DESIGN  FACTOR 


where 


hj  =  slab  thickness 

Ej  ■  concrete  modulus 

Ecomp  =  composite  modulus  subgrade  modulus 

LOAD  =  loading  parameters 

and  combining  equations  2.3  and  2.4  lead  to: 

4nN  •  6.4887  - — -  .  3.82192  (2.5) 

0(hl-  VEcof  WAD> 


It  is  seen  that  the  effect  of  varying  E  on  N  is  not  a  simple  function. 

comp  r 

It  will  vary  with  all  pavement  and  load  variables.  Equation  2.5  is  used  to 


compute  the  deviation  of  N  due  to  a  deviation  of  Ecojnp  in  the  following  way: 


(a)  For  a  given  aircraft,  slab  thickness,  modulus  and  Ecomp>  compute 


V 


(b)  Compute  Nj  and  N2  at  (0.8  Ecomp)  and  (1.2  Ecomp)  corresponding  to 

20%  deviation  is  E  _ ; 

comp  N 

(c)  Compute  the  relative  deviation  100  (N2  -  N^)  +  (-^)  ; 

(d)  Repeat  computations  (b)  and  (c)  for  different  percentage  15,  10  and 

5%  of  E  _  deviations; 

comp 

(e)  Repeat  computations  for  different  sets  of  pavement  variables. 

(f)  Repeat  computations  for  different  aircrafts. 


Results  of  this  analysis  are  given  in  Tables  2.1. a  to  2.1.e  for  5  dif¬ 
ferent  aircrafts  AG-2,  AG-9,  AG- 10,  AG-12  and  AG-13  where  Ej  and  MR  are  as¬ 
sumed  equal  to  4,000,000  and  700  psi  respectively.  It  can  be  seen  that  the 
relative  variation  of  N  for  N  values  between  100  and  100000  is  1.25  to  5.0 
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times  that  of  E 


.  The  magnifying  factor  increases  with  increasing 


comp 

N  (corresponding  also  to  increasing  slab  thickness  and/or  to  increasing 

Ecomp) *  *s  *arSer  ^or  the  heavier  aircrafts  (AG-9  to  13)  than  for 

the  lighter  aircraft  (AG-2).  Therefore,  it  is  concluded  that  the 

design  procedure  is  very  sensitive  to  E^  and  more  attention  should 

be  given  to  the  determination  of  E  _ .  Care  should  be  taken  in  the 

comp 

use  of  E  since  there  is  no  maximum  value  of  E  which  may  be  used 
comp  comp 

as  there  is  with  the  composite  modulus  of  subgrade  reaction  (k  =  500  pci). 

max 

It  should  be  noted,  however,  that  h  is  still  the  major  governing  factor. 

For  example,  if  the  computation  of  N  is  kept  within  50%  accuracy, 

determination  of  E  should  be  limited  within  less  than  40%  accuracy 
comp 

for  the  light  aircraft  and  within  less  than  10%  accuracy  for  the  heavier 
aircrafts.  The  above  conclusion  must  be  taken  into  account  in  the  design 
procedure  and  any  applications  of  the  procedure  as  discussed  in  the 
following  paragraphs. 
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TABLE  2.1a: 


Effect  of  Composite  Modulus 
Deviation  on  Number  of  Cover¬ 
ages  for  AG-2 


Slab 

Thickness 

in. 

Composite  or 
Subgrade 
Modulus , 
psi 

Number  of 
Coverages -N 

Relative  Deviation  of  N  for  Dif¬ 
ferent  Relative  Derivations  of 

E  % 

comp 

20  15  10  5 

6 

50,000 

25 

31 

23 

16 

8 

8 

50,000 

980 

42 

31 

21 

10 

10 

50,000 

82,620 

55 

41 

27 

14 

S 

20,000 

190 

31 

23 

15 

8 

10 

20,000 

9500 

41 

31 

20 

10 

12 

20,000 

1,070,000 

54 

40 

27 

13 

8 

10,000 

73 

25 

18 

12 

6 

10 

10,000 

2,650 

33 

25 

16 

8 

12 

10,000 

203,000 

43 

32 

21 

11 

8 

5,000 

34 

19 

14 

10 

5 

10 

5,000 

967 

25 

19 

13 

6 

12 

5,000 

54,950 

32 

24 

16 

8 

TABLE  2.1b:  Effect  of  Composite  Modulus  Deviation 
on  Number  of  Coverages  for  AG- 9 


Slab 

Thickness, 

in. 

Compute  or 
Subgrade 
Modulus , 
psi 

Number  of 
Coverage-N 

Relative  Deivation  of  N  for  Dif¬ 
ferent  Relative  Deviations  of 

E  % 

comp 

20  15  10  5 

8 

50,000 

200 

64 

48 

32 

16 

10 

2,100 

84 

61 

40 

20 

12 

50,000 

2,450 

105 

76 

49 

24 

14 

50,000 

322,000 

130 

91 

58 

28 

8 

20,000 

16 

47 

35 

23 

11 

10 

20,000 

95 

57 

42 

28 

14 

12 

20,000 

630 

67 

49 

33 

16 

14 

20,000 

4,700 

78 

57 

37 

19 

16 

20,000 

41,000 

90 

65 

42 

21 

18 

20,000 

427,000 

102 

73 

47 

23 

10 

10,000 

18 

41 

30 

20 

10 

12 

10,000 

92 

48 

36 

24 

12 

14 

10,000 

520 

56 

41 

27 

14 

16 

10,000 

3,500 

63 

46 

31 

15 

18 

10,000 

28,200 

69 

51 

33 

17 

20 

10,000 

292,000 

74 

54 

36 

18 

12 

5,000 

22 

37 

28 

18 

9 

14 

100 

42 

32 

21 

11 

16 

560 

47 

35 

23 

12 

18 

3930 

50 

37 

24 

12 

20 

38,400 

49 

36 

24 

12 

14 


TABLE  2.1c:  Effect  of  Composite  Modulus  Deviation 
on  Number  of  Coverages  for  AG-10 


Slab 

Thickness, 

in 

Composite  or 
Subgrade 
Modulus, 
psi 

Number  of 
Coverages 

N 

Relative  Deviation  of  N  for  Dif¬ 
ferent  Relative  Deviations  of 

E  ,  % 

comp’ 

20  15  10  5 

8 

50,000 

5,000 

89 

65 

43 

21 

10 

50,000 

102,500 

119 

85 

55 

27 

8 

20,000 

170 

64 

48 

31 

16 

10 

20,000 

1,580 

80 

59 

39 

19 

12 

20,000 

16,400 

99 

71 

46 

23 

14 

20,000 

188,000 

120 

85 

55 

27 

8 

10,000 

26 

48 

36 

24 

12 

10 

10,000 

157 

60 

45 

29 

15 

12 

10,000 

1,040 

73 

54 

35 

18 

14 

10,000 

7,470 

87 

63 

41 

20 

16 

10,000 

59,200 

101 

73 

47 

23 

18 

10,000 

520,000 

116 

82 

53 

26 

10 

5,000 

26 

46 

35 

25 

11 

12 

5,000 

121 

55 

41 

27 

13 

14 

5,000 

630 

63 

46 

31 

15 

16 

5,000 

3,800 

68 

50 

33 

16 

18 

5,000 

26,600 

70 

51 

54 

17 

20 

5,000 

236,000 

66 

48 

31 

15 

TABLE  2. Id:  Effect  of  Composite  Modulus  Deviation 
on  Number  of  Coverages  for  AG-12 


Slab 

Thickness, 

in. 

Composite  or 
Subgrade 
Modulus, 
psi 

Number  of 
Coverages 

N 

Relative  Deviation  of  N  for  Dif¬ 
ferent  Relative  Deviations  of 

E  ,  % 

comp 

20  15  10  5 

1  $ 

1,350 

71 

53 

35 

17 

16,800 

92 

67 

44 

22 

14 

50,000 

243,000 

117 

84 

54 

26 

10 

20,000 

82 

53 

39 

26 

13 

12 

20,000 

550 

65 

48 

32 

16 

14 

20,000 

4,060 

79 

58 

38 

19 

16 

20,000 

33,000 

94 

68 

45 

22 

18 

20,000 

294,000 

112 

80 

52 

25 

10 

10,000 

17 

40 

30 

20 

10 

12 

10,000 

80 

50 

37 

24 

12 

14 

10,000 

410 

60 

44 

29 

15 

16 

10,000 

2,300 

71 

52 

34 

17 

18 

10,000 

13,600 

83 

61 

40 

20 

20 

10,000 

87,000 

96 

70 

45 

22 

12 

5,000 

18 

39 

29 

20 

10 

14 

5,000 

68 

47 

35 

23 

12 

16 

5,000 

276 

55 

41 

27 

14 

18 

5,000 

1,200 

64 

47 

31 

16 

20 

5,000 

5,700 

72 

53 

35 

17 

22 

5,000 

29,200 

79 

58 

38 

19 

24 

5,000 

166,000 

85 

62 

40 

20 

TABLE  2.1e:  Effect  of  Composite  Modulus  Deviation 
on  Number  of  Coverages  for  AG- 13 


Slab 

Thickness 

in. 

Composite  or 
Subgrade 
Modulus , 
psi 

Number  of 
Coverages 

N 

Relative  Deviation  of  N  for  Dif¬ 
ferent  Relative  Deviation  of 

E  ,  % 

comp’ 

20  15  10  5 

Hi 

50,000 

52 

52 

39 

26 

13 

50,000 

277 

62 

46 

30 

15 

50,000 

1,550 

72 

53 

35 

17 

17 

50,000 

9,400 

82 

60 

39 

20 

19 

50,000 

63,300 

94 

68 

44 

22 

13 

20,000 

28 

41 

31 

21 

10 

15 

20,000 

112 

48 

35 

23 

12 

17 

20,000 

488 

54 

40 

26 

13 

19 

20,000 

2,330 

60 

45 

30 

15 

21 

20,000 

12,500 

68 

50 

33 

16 

23 

20,000 

77,700 

76 

55 

36 

18 

15 

10,000 

27 

35 

26 

18 

9 

17 

10,000 

100 

40 

30 

20 

10 

19 

10,000 

400 

45 

33 

22 

11 

21 

10,000 

1,780 

50 

37 

24 

12 

23 

10,000 

9,080 

55 

40 

27 

13 

25 

10,000 

55,300 

60 

44 

29 

15 

17 

5,000 

30 

31 

23 

16 

8 

19 

5,000 

103 

35 

26 

17 

9 

21 

5,000 

402 

38 

28 

19 

9 

23 

5,000 

1,800 

41 

31 

20 

10 

25 

5,000 

9,750 

43 

32 

22 

11 
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Determination  of  the  Composite  Modulus 


The  composite  modulus  is  determined  on  the  basis  of  equal  maximum  ten¬ 
sile  stress,  i.e.  the  same  maximum  tensile  stress  will  occur  when  analyzing 
either  the  original  system  composed  of  a  concrete  layer  on  a  base-subbase- 
subgrade  system  or  an  equivalent  system  composed  of  the  same  concrete  layer 
on  the  top  of  the  "equivalent"  composite  subgrade.  Like  the  ESWL,  all  pave¬ 
ment  variables  (h^,  E^)  and  the  loading  (gear)  variables  may  be  expected  to 
affect  the  composite  modulus.  In  the  following.it  was  assumed  that  the  Econip 
function  was  of  the  following  mathematical  form. 


E  =  E_  f, 
comp  3  1 


(2.6) 


where 


fj  =  a  function  of  the  pavement  geometry  only,  under  the 
application  of  one  wheel  load  only 

f2  *  a  function  of  the  loading  (gear)  conditions  only 


Composite  Modulus  for  One  Wheel  Load 

The  following  range  of  variables  were  included  in  the  computation  of 
the  maximum  tensile  stress  and  computation  of  the  Ecomp  for  one  wheel  load. 

(a)  Wheel  load  of  AG-2,  AG-4,  AG-6,  AG-9,  AG-10,  AG-12  and  AG-13; 

(b)  hj  -  thickness  of  concrete  layer  varying  from  6  in.  to  25  in.; 
for  each  aircraft  group,  four  thicknesses  covering  the  whole  range  of  ex¬ 
pected  design  thicknesses  (see  Table  2.2); 

(c)  h2  -  thickness  of  base-subbase  layer  of  6  and  12  in; 

(d)  E2  -  modulus  of  base-subbase  layer  of  50,000,  200,000  and  800,000 
psi  to  include  stabilized  material,  as  well  as  unbound  materials; 

(e)  Ej  -  modulus  of  subgrade  of  4000,  12,000  and  35,000  psi. 
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The  form  of  the  function  f  j  was  chosen  to  include  the  limiting  case  of 
1*2=0  and  E2=E3  where  no  base-subbase  layer  exists.  Several  functions  of  the 
pavement  geometry  variables  were  assumed.  The  following  equation  gave  the  most 
satisfactory  correlation  coefficient  and  standard  error  of  estimate: 


£n(E  /E,)  =  a.  Zn==-  * 

v  comp  3-j  1  Ej 


E  E 

Ina  (a2£nh1  +  a^-dn^i  +  a^nfl^+l)  ♦  a5-Cn^- 
3  3 


E  E 

+  agfinhp2  +  a^dnl^&i^i  ♦  ajjCfn— )  )  ♦ 


2  2 
(£na)  (a^dnh^  +  a^lnO^  ♦  1)  ♦  ag£ng=- } 


(2.7) 


where 


“  -  V  V  01 2  +  15 

(E  /E_)  »  the  ratio  between  E„  _  and  E,  for  one 

comp  3  j  comp  3 

wheel  load 

al  ***  all  ”  'egression  coefficients 

ax  *  0.434783,  a2  *  -0.134317,  a3  «  0.0816432, 

a.  =-0.0807450,  a.  =  -0.190090,  a,  =-0.164512, 

4  5  o 

a?  «  0.119533,  ag  =  0.0550536,  ag  =  0.0943027, 

a1Q=  0.0102537,  an«  -0.0118395 

7 

The  R  of  the  multiple  correlation  was  0.994  and  the  standard  error  of 
estimate  0.0776  corresponded  to  a  standard  deviation  of  8.0%  in  the  Ecomp  ~ 
evaluation.  It  should  be  noted  that  the  equation  has  a  zero  intercept  and 


.9 


fno,  &»=- 
fc3 


equal  0  when  h2  =  0  and  E2  *  E^. 


A  simpler  equation,  with  5  terms  only  was  also  derived: 


-tnffj/Ej)  =  (£na)  {aj  £nhj  ♦  a2  £n(h2+l)  ♦  a^  in—  }  ♦  (2.8) 

E  2 

£na  {  a^C&ih^2  ♦  a^(£n^)  } 

where 

a|  =  -0.161317,  a2  =  0.0977117,  a^  =  0.0358696 
a'4  =  0.0802967,  a^  =-0.00196301 

2 

with  R  =  0.991  and  standard  error  of  estimate  =  0.095  corresponding  to  a  stan¬ 
dard  deviation  of  10%  in  the  E  evaluation.  It  should  be  noted  that  (1)  the 

contribution  of  the  terms  with  h^  and  E^/E^  to  the  regression  was  relatively 
substantial,  partly  due  to  the  wide  range  of  h^  and  E^/E^  included  in  the 
analysis;  (2)  the  general  trend  of  the  effect  of  pavement  geometry  on  the  E 

comp 

parameter  was  as  follows:  (a)  E  /E,  increases  as  h_  and/or  E-./E,  increases; 

comp  b  i  ib 

(b)  Ecomp/E3  decreases  as  hj  increases  and/or  E^/E^  decreases;  (3)  within  the 
range  of  the  wheel  loads  used,  it  seems  that  the  relationship  is  independent 
of  the  radius  of  the  contact  area;  (4)  for  the  378  data  points  used  in  the 
regression  the  deviation  in  evaluating  E^q  with  the  regression  equation  ex¬ 
ceeded  25%  in  some  but  few  cases.  Therefore,  depending  on  the  accuracy  level 
required  for  predicting  the  number  of  coverages,  equation  2.7  may  or  may  not 
be  adequate  for  design  purposes.  According  to  the  above  analysis  of  the  ef¬ 
fect  of  composite  modulus  deviation  on  number  of  coverages,  a  deviation  of  10 

percent  in  E  evaluation  resulted  in  about  25  and  40  percent  deviation  in 
r  comp 

the  expected  number  of  coverages  for  the  light  and  heavy  loads  respectively. 

A  procedure  for  improving  the  composite  modulus  evaluation  was  developed. 
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based  on  only  one  run  of  the  BISAR  program.  The  procedure,  presented  in 
Chapter  3  must  be  used  at  least  for  all  heavy  aircrafts;  (5)  equation  2.8 


predicts  the  dependence  of  Ecomp  upon  the  pavement  geometry  variables  quite 
well.  Because  it  is  easier  to  differentiate,  this  equation  was  used  in  the 


probabilistic  approach,  in  the  evaluation  of  the  variance. 


Composite  Modulus  for  Different  Aircrafts 


The  effect  of  loading  (gear)  conditions,  that  is,  the  number  of  wheels 


and  wheel  configuration  is  expressed  in  Equation  2.6  by  the  ^-function.  For 
one  wheel  load,  f 2  =  1  (i.e.  for  AG-1,  2  and  3).  The  loading  conditions  seem 
to  have  a  similar  effect  on  E  as  h, ,  the  concrete  layer  thickness.  From 
previous  paragraphs,  it  was  stated  that  while  all  other  variables  are  kept  con¬ 


stant,  increasing  the  h^  resulted  in  a  decrease  of  the  Ecomp  term.  This  can 
be  interpretated  as  follow:  increasing  h^  increases  the  pavement  stiffness, 
the  radius  of  relative  stiffness  (£)  and  the  size  of  the  deflection  bowl;  The 


relative  base-subbase  layer  thickness  h decreases,  entraining  a  reduction 

of  its  effect  on  E  .  The  same  relation  between  E  and  number  of  wheels 
comp  comp 

exists:  when  the  number  of  wheel  loads  increases,  the  deflection  bowl  size 


increases,  and  the  relative  base-subbase  layer  thickness  decreases.  There¬ 
fore,  the  f2 -function  may  logically  be  expected  to  decrease  as  the  number  of 
wheels  increases.  Instead  of  expressing  f^  for  each  aircraft  group,  it  was 
decided  to  simplify  the  analysis  of  f2  and  confine  it  to  three  types  of  wheel 
loading  configuration,  (namely  two,  four  and  six  gear  wheels).  Figures  2.5 
to  2.9  show  the  relationship  between  the  Ecojnp/Ej  parameter  for  one  wheel 

load  and  E  /E-  for  AG-4,  6,  9,  13  and  10.  A  simple  log-log  function  with- 
comp  o 

out  intercept  was  fitted  to  the  computation  results  for  the  two,  four  and  six 
wheel  assemblies,  (AG-4  §  6,  9  8  13  and  10  respectively): 


comp  3  two  wheels 


^comp^3^one  wheel 

Figure  2.5  THE  f2  -  FUNCTION  FOR  AG-4 
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Figure  2.6  THE  f-  -  FUNCTION  FOR  AG-6 


^comp^3^one  wheel 

Figure  2.7  THE  f2  -FUNCTION  FOR  AG-9 


Figure  2.8  THE  f,  -  FUNCTION  FOR  AG-13 


3}six  wheels 


1  2  4  6  8  10  20  30 


^comp^3^one  wheel 


Figure  2.9  THE  -  FUNCTION  FOR  AG-10 


a 

n 


£n(E  /E,)  *  a  in(E  /ET)  (2.9) 

v  comp  3'n  n  v  comp  3  j  1  J 

where 

(E  /E_)  *  the  E„  /E_  for  n-wheels 

v  comp  3 comp  3 

an  =  regression  coefficients,  for  n-wheels 

a,  «  1.0,  a,  =  0.889813,  a,  =  0.730507,  a,  =  0.695222 

12  4  O 


The  k  -coefficients  were  0.984,  0.987,  0.975  and  the  standard  errors  of 
estimate  were  0.073,  0.0493,  0.0704  for  2,  4  and  6  wheels  respectively.  It 
is  noted  that  (1)  the  composite  modulus  is  strongly  affected  by  the  loading 
conditions,  reaching  about  half  the  value  of  one  wheel  for  the  4  and  6-wheel 
loading  conditions;  (2)  the  relationship  f2  could  be  improved  if  some  pavement 
geometry  variables  were  included  in  equation  2.9:  The  an-regression  coeffici¬ 
ent  seems  to  increase  as  hj  increases.  However,  it  is  believed  that  equation 
2.9  is  quite  accurate  and  there  is  no  need  for  expanding  it. 


Evaluation  of  the  Results 


The  results  of  the  above  computations  are  compared  with  those  of  the  cur¬ 
rent  USACE  procedure  for  determining  the  composite  subgrade  modulus  of  reac¬ 
tion.  A  summary  of  the  results  is  presented  in  Table  2.3  where:  (1)  k  - 
the  subgrade  modulus  of  reaction  is  computed  according  to  Parker  et  al.  (1 ) , 


using: 


log  ksg  *  (log  Ej  -  1 .415)/! . 284 


(2.10) 


(2)  k  -  the  composite  subgrade  modulus  of  reaction  is  evaluated  from  Fig- 
'  '  comp 

ure  2.3;  (3)  k(E  )  is  evaluated  using  equation  2.10  and  substituting  E 
comp 

for  E,;  (4)  the  range  of  k(E  _ )  includes  all  aircrafts  and  concrete  layer 

3  coop 

thicknesses  given  in  Table  2.2.  It  should  be  noted  that  the  upper  limit  cor¬ 
responds  to  the  light  aircraft  (AG-2)  and  very  thin  concrete  layer  (6")  while 
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TABLE  2.2:  Summary  of  Concrete  Layer  Thickness 
Used  in  the  Analysis  of  Composite 
Modulus  for  One  Wheel 


Aircraft  Grou 

■ 

Concrete  Layer 
Thickness,  in. 

2 

6 

4 

7,10,12,15 

6 

7,10,13,17 

9 

8,12,15,19 

10 

8,12,15,20 

12 

10,15,19,24 

13 

11,16,20,25 

TABLE  2.3:  Comparison  of  Composite  Subgrade 
Moduli  of  Reaction 


h2,in 

E2,psi 

Ej.psi 

W  *pci 

k  ,pci 

comp  r 

k(E  ) ,pci 

comp' * v 

6 

50,000 

35,000 

275 

360 

278-295 

12,000 

120 

160 

124-154 

4,000 

50 

50 

53-  75 

200,000 

35,000 

275 

490 

307-389 

12,000 

120 

220 

134-207 

4,000 

50 

60 

65-113 

800,000 

35,000 

275 

(970) 

315-584 

12,000 

120 

390 

165-349 

4,000 

50 

120 

78-226 

12 

50,000 

35,000 

275 

(550) 

282-317 

12,000 

120 

134-202 

4,000 

50 

65 

63-120 

200,000 

35,000 

275 

342-577 

12,000 

120 

340 

174-397 

4,000 

50 

105 

78-285 

800,000 

35,000 

275 

(1270) 

428-1352 

12,000 

120 

(580) 

242-1125 

4,000 

50 

200 

133-984 

Note:  Numbers  in  parentheses  are  obtained  from  extrapolation 
in  Figure  2.3  .  The  value  should  be  500  pci 
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the  lower  limit  corresponds  to  heavy  aircrafts  and  thick  concrete  layer. 

The  values  of  ^comp  for  weak  subgrade  were  evaluated  by  extrapolating 
the  results  of  Figure  2.3.  For  the  case  of  E, = 4000  psi  (k  =  50  pci)  and 

J  S5 

E-  *50,000  psi,  higher  k  values  are  obtained  from  Figure  2.2  than  from 
1  comp 

extrapolating  Figure  2.3. 

Comparison  of  the  k  and  k(E  )  in  Table  2.3  shows  that:  (1)  for 
r  comp  v  compJ  v  ' 

weak  subgrades,  the  ^(^comp)  *-s  either  in  the  range  of  or  higher  than  the 

^comp  eva^uate<^  with  the  current  USACE  procedure.  However,  if  the  extreme 

cases  of  very  thin  concrete  layer  were  excluded  from  k(E  )  and  the  k 

were  derived  from  Figure  2.2,  the  discrepancy  would  be  quite  substantially 

reduced;  (2)  for  relatively  strong  subgrades,  the  k(Ecomp)  is  lower  than  the 

^comp*  exceeds  500  pci  in  few  cases  only.  This  result  seems  to  support 

the  current  USACE  procedure  where  the  subgrade  modulus  of  reaction  (or  the 

composite)  is  limited  to  a  maximum  value  of  500  psi.  It  appears  that  the 

computation  results  of  the  composite  modulus  are  in  general  agreement  with  the 

current  USACE  k  value,  ahthough  more  research  studies  are  necessary  to  fully 
IDcLX 

support  this  general  observation.  Furthermore,  the  present  derivation  takes 
into  account  all  of  the  design  variables  which  appear  to  have  different  effects 
on  the  pavement  performance. 

Summary  and  Conclusions 

A  composite  modulus  computation  scheme  has  been  presented  that  includes 
all  pavement  geometry  and  loading  condition  variables.  It  was  found  that  the 
composite  modulus  of  elasticity  is  not  only  a  function  of  the  base-subbase  and 
subgrade  layers,  but  also  a  function  of  the  concrete  layer  and  load  configura¬ 
tion. 


The  study  of  the  effect  of  E  on  allowable  number  of  coverages  (or 

comp 

pavement  performance)  clearly  showed  that  the  degree  of  accuracy  achieved  in 
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the  design  (predicted  number  of  coverages)  is  strongly  related  to  that  achieved 

in  E  _  evaluation.  The  relative  deviation  in  the  allowable  number  of  cover- 

comp 

ages  is  1.2  to  5  times  the  relative  deviation  in  E 

comp 

These  two  results-  of  the  function  Econp  being  dependent  upon  all  vari¬ 
ables  and  that  of  the  sensitivity  analysis,  emphasize  the  importance  of  the 
evaluation  of  Ecomp»  and  give  quantitative  relationships  to  determine  the  re¬ 
quired  degree  of  accuracy. 
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Chapter  3 


PROBABILISTIC  ANALYSIS  OF  PCC 
AIRFIELD  DESIGN  USING  ELASTIC 
LAYERED  THEORY 

Introduction 

The  analysis  presented  in  this  chapter  is  based  on  the  mathematical 
formulation  given  in  Volume  II,  namely  the  approximate  closed-form 
probabilistic  approach  and  the  probabilistic  simulation  approach.  Both 
methods  require  stress  computations  which  can  be  handled  with  computer 
programs  such  as  the  Shell  BISAR  code.  However,  while  the  use  of  a 
computer  program  is  justified  for  the  design  of  a  particular  pavement, 
it  does  become  unrealistic  and  uneconomical  for  the  probabilistic  approach, 
due  to  the  large  computer  time  needed.  The  analysis  was  therefore  confined 
to  specific  critical  aircraft  type  included  in  the  military  (USAF)  aircraft 
classification  scheme  (AG-1  to  13).  Solutions  for  other  aircraft  type  may 
be  developed  by  following  the  procedure  noted  below  and  developing  an 
equation  for  the  maximum  tensile  stress. 

The  probabilistic  approach  presented  in  the  following  paragraphs 
includes: 

(1)  Stress  computations  and  derivation  of  an  equation  for  the  maxi¬ 
mum  tensile  stress  for  each  of  13  aircraft  types:  USAF  classification 
AG-1  to  13; 

(2)  Derivation  of  the  relationships  between  variances  of  the  depend¬ 
ent  and  independent  variables  for  the  approximate  closed-form  approach. 

The  linear  or  first  order  Taylor  series  expansion  is  assumed  as  presented 
in  Volume  II. 

(3)  Formulation  of  a  simulation  model  which  bypasses  the  linear 
assumption  in  the  closed-form  approach. 
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(4)  Run  examples,  to  compare  between  the  closed-form  approach  and 
the  simulation  one,  at  different  variation  levels  of  the  design  parameters. 


Stress  Computations 

In  the  "normal"  linear  elastic  layered  system,  the  stress  at  any  point 
is  a  function  of  the  layer  thicknesses  (expressed  in  terms  of  the  radius 
of  contact  area)  and  of  the  modular  ratios  of  the  layers  and  the  subgrade. 
However,  because  of  the  assumptions  made  by  Parker  et  al .  (1)  concerning 
the  depth  of  the  subgrade,  the  modulus  of  the  rigid  layer  underneath  the 
subgrade  (chosen  to  be  equal  to  1, 000,000  psi),  and  the  value  of  the 
friction  coefficient  at  the  interface  between  the  concrete  and  base  on 
subgrade  layers,  it  is  not  possible  to  use  dimensionless  variables.  There¬ 
fore  stress  computation  were  conducted  for  all  13  aircrafts  of  the  USAF 
classification  scheme  (AG-1  to  13)  for  the  following  ranges  of  variables 
(see  Table  3.1) : 

(1)  Modulus  of  concrete  -  *  4,000,  000  psi  and  Poisson's  ratio 

of  0.2; 

(2)  Concrete  layer  thickness  -  hj ,  choosen  to  cover  a  wide  range  of 
design  possibilities.  Four  different  values  were  assigned  for  each  air¬ 
craft  group,  as  shown  in  Table  3.1; 

(3)  Subgrade  moduli  or  composite  subgrade  moduli  of  4,000,  12,000, 
35,000  and  100,000  psi  and  Poisson's  ratio  of  0.4. 

(4)  Subgrade  thickness  of  20  ft  =  240  in. 

(5)  Elastic  modulus  of  the  stiff  infinite  layer  beneath  the  subgrade 
of  1,000,000  psi. 

(6)  Friction  coefficient  (interface  compliance)  -at  the  interface  of 
the  concrete  and  subgrade  layer  of  1 ,000, (defined  in  the  BISAR  program). 
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Table  3.1  Summary  of  Input  Data  for 
Stress  Computations 


Radius  of  Contact 

Contact  area  Pressure 
in.  psi 


9. 3S60 
S.6419 
8. 7S86 

11.2838 


8.68S6 


98.1818 

270.0 

186.722 

87.1875 
155. 4545 
156.9684 
190.4008 


Wheel  Coordinates 


Wheel  No 


1  0.0  0. 

1  0.0  0. 

1  0.0  0. 


8.3302  183.0275 


8.1369  186.6587 


9.6738  180.3061 


8.8310  188.5714 


9.2189  233.7079 


9.5246  105.6728 


0.0  0. 

34.5  0. 

0.0  56. 

34.5  56. 

0.0  0. 

32. 5  0. 
0.0  48. 


Max. Stress  Concrete  Layer 
Under  Thicknesses 

Wheel  No  in. 


6  9  11  14 

6  9  11  14 

7  10  12  15 

7  10  12  15 

7  10  13  17 

7  10  13  17 

8  11  14  18 

8  12  15  19 


3  -53.0 

0.0 

4  -87.0 

0.0 

5  -2.50 

65.0 

6  -50.5 

65.0 

8  12  15  19 


8  12  IS  20 


10  15  19  24 


11  16  20  25 


8  12  15  20 
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The  general  form  of  the  maximum  tensile  stress  at  the  bottom  of 


the  concrete  layer  is  given  by: 

o  =  o(hj,  Ej,  Ecomp»  loading  conditions)  (3.1) 

where  Ecpm  =  the  composite  modulus  of  subgrade  is  equal  to  the  subgrade 
modulus  when  no  base-subbase  layer  exists. 

The  computation  results  for  each  aircraft  were  used  to  develop  a 
regression  equation  for  prediction  the  maximum  tensile  stress  for  any  com¬ 
bination  of  pavement  geometry  variables.  The  following  model  fitted  all 
thirteen  aircrafts: 


a  ~  a.  +  a^  +  a„E  +  a^.Eh  +  a„„fIE  h  + 


EH 


E2H 


£n8  /0  £nB  ^  ,£nB.2  2  ,,  .... 

31  IT  *  V8  a3  TST  *  a4  ‘IT 1  *  V6  (3'2) 


where 


-h  =  h. 


-E  =  E./E 

1  comp 

-B  *  h3*£ 

-hj=  concrete  layer  thickness,  inch. 

-Ej=  modulus  of  concrete,  psi. 

-Ecojup  =  composite  or  subgrade  modulus,  psi. 
-a  «  maximum  tensile  stress,  psi. 

-a^=  regression  coefficients. 


2 

The  values  of  the  regression  coefficient,  R  and  the  relative  error 
of  estimate  are  given  in  Table  3.2.  It  is  seen  that  the  regression  equa¬ 
tion  is  an  excellent  one,  giving  a  maximum  standard  error  of  0.8%. 


Comput  a 


224.6386  -4.946335  0.07221752  -0.003123284  0.  0.  -5139.840  -75798.53  7112.995  460427.8  I  .99991  0.54 


Although  equation  3.2  contains  quite  a  large  number  of  terms  in  order 
to  achieve  the  high  degree  of  accuracy  needed,  it  is  quite  easily  pro¬ 
grammable  on  micro-computers  and  calculators.  It  should  be  noted  that: 

(a)  the  maximum  tensile  stress  was  found  to  be  higher  for  a  6-wheel  gear 
than  for  the  12-wheel  gear  of  the  C-5  (AG-10);  (b)  maximum  tensile  stress 
was  found  to  correspond  to  a  4-wheel  gear  in  some  combinations  of  h  and 
E  and  to  an  8-wheel  gear  in  other  combinations  for  the  B-747  (AG-12). 
Equation  3.2  expresses  the  maximum  stress  corresponding  to  either  4  or  8- 
wheel  gear. 

Equation  3.2  covers  a  quite  wide  range  of  values  of  the  pavement 
design  parameters  and  design  analysis.  It  should  be  remembered  that 
equation  3.2  applies  for  a  two-layer  system  only.  When  the  layered  system 
underneath  the  concrete  contains  a  base-subbase  layer,  it  should  be  trans¬ 
formed  into  an  equivalent  subgrade  and  its  composite  modulus  evaluated 
according  to  Chapter  2  of  this  volume. 


The  Approximate  Closed  Form  Probabilistic  Approach 
In  the  approximate  closed-form  probabilistic  approach,  the  average  and 
the  variance  of  the  random  variables  are  given  by  the  following  equations, 
assuming  linear  or  first  order  TAylor  series  expansion  of  the  variable 
function: 

E[g(x.)J  «  g(x.)|_  (3.3a) 

Xi 

2  n/3g(x.)j  \2 

E[(g(x.)-u)  ]  =  Var[g(xi)]  -  M -57; — I  /  Var[x.]  ♦ 

i=l  i  'x. 


Covar  [xi,xk] 
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(3.3b) 


where 


g(x^)  =  random  variable  function  of  variables 
xi  =  average  of  x^ 

Higher  orders  of  expectation  of  the  random  variable  can  also  be 
computerd,  for  a  complete  description  of  the  probability  density  distri¬ 
bution.  In  the  following,  it  was  assumed  that  the  random  variable  is 
normally  distributed  and  therefore,  the  higher  order  of  expectations  are 
not  needed.  Pavement  performance,  expressed  in  terms  of  number  of  coverages, 
is  not  likely  normally  distributed  but  rather  closer  to  a  log-normal  distri¬ 
bution.  The  number  of  coverages  can  not  be  chosen  as  the  independent 
variable.  The  next  variable  in  the  design  sequence  is  the  design  factor 
whose  distribution  can  be  approximated  by  the  normal  one.  Through  a  change 
of  variable,  this  will  lead  to  log-normal  distribution  for  the  number  of 
coverages.  Therefore,  g(x^)  in  equation  3.3  represents  DF,  and  are  the 
design  variables  namely  MR,  hj ,  E^,  h^,  E3  and  LOAD.  The  design  factor 
is  expressed,  for  each  aircraft: 


DF  =  — 

o 


(3.4a) 


o  =  a  (h.,  E  /E  ) 
l  1  comp 


(3.4b) 


(%]  *  {i  h>-h4-  if  ) 


(3.4c) 


where  the  functions  a,  fj  and  f2  are  given  in  equations  3.2,  2.7  and  2.9 
respectively. 

Assuming  that  only  and  MR  are  correlated  and  substituting  equations 
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3.4  into  3.3b  lead  to  the  following  expression  of  the  coefficient  of 
variation  of  DF: 

CV2[DF]  =  CV2[MR]  ♦  Z  (  xi  3o  I  .  CV  [x.j\2  + 

i=l\a  3xJ  __  1 1 

x. 

i 

-  fl  3a_  p[E. ,MR] .CV[MR]  .  CV[E.]  (3.5) 

o  3Ej  1  1 

where  x  *  hj,  E^,  E3 


Computation  of  the  derivatives  of  DF  with  respect  to  x^  for  a  given 


aircraft  are  as 

follows : 

3DF 

MR 

3o 

(3.6a) 

ah® 

‘  2 

3h. 

1 

a 

1 

3DF 

MR 

.  3  a 

(3.6b) 

3h_  * 

'  2 

3h_ 

2 

a 

2 

3DF 

3E  * 

MR 

2 

.  3o 

3E, 

(3.6c) 

1 

a 

1 

3DF 

MR 

3o 

3E 

corap 

(3.6d) 

"SET  " 

’  2 

3E 

3E. 

2 

0 

comp 

2 

3DF 

MR 

.  3a 

3E 

.  comp 

(3.6e) 

■SET  = 

~  2 

3E 

3E_ 

3 

a 

comp 

3 

Using 

equations  3.2 

and  2.8,  the  following  final 

equations  are 

obtained : 
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The  probabilistic  analysis  of  PCC  airfield  designs,  based  upon 
multilayered  elastic  theory,  relies  upon  the  determination  of  the  probability 
density  distribution  of  N  for  a  given  aircraft  type  and  given  pavement 
geometry.  The  general  solution  sequence  is  summarized  below: 

(1)  Computer  average  stress  (c)  using  equation  3.2  and  average  values 
of  the  design  variables.  (BISAR  could  have  been  used  in  the  calculation  of  o; 
however,  equation  3.2  has  been  used  within  the  program  in  order  to  reduce  the 
additional  computer  time  which  would  have  been  required  to  run  BISAR.  BISAR 
is  used  later  in  the  program  if  a  corrected  composite  modulus  is  required.) 


(2)  Compute  average  DF  (DF)  as: 

DF  =  MR/~  (3.9) 

(3)  Compute  the  variance  of  DF  using  equation  3.5  and  the  derivations 
in  equations  3.6  to  3.8. 

(4)  Using  DF  and  Var[DF],  compute  the  cumulative  distribution  of  DF 
from  the  normal  distribution.  Subdivide  the  interval  of  DF  _+  3o  into  say 
30  DFy  values  and  compute  the  cumulative  distribution  Py  (DF  <  DFy)  corre¬ 
sponding  to  each  DF^. 

(5)  Compute  corresponding  to  DF^  using  equation  2.3.  The 

is  related  to  Py,  giving  the  cumulative  distribution  of  number  of  coverages. 

This  scheme  is  implemented  in  a  computer  program  (see  Appendices  I  and 
II),  and  makes  use  of  the  regression  equations  developed  for  this  purpose. 
However,  it  was  shown  that  the  degree  of  accuracy  achieved  in  the  prediction 
of  the  number  of  coverages  depends  upon  the  degree  of  accuracy  in  the  com¬ 
posite  modulus  and  stress  computations.  Because  most  of  the  error  is  due 

to  the  E  evaluation,  a  correction  option  for  reducing  the  error  was 
comp 

included  in  the  program.  The  correction  scheme,  which  corrects  the  value  of 

E  _  to  give  an  exact  tensile  stress  is  as  follows: 

comp 
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(1)  Input  average  design  variables  h^,  E^,  E ^  and  and  the 

design  load  using  the  format  of  the  BISAR  program. 

(2)  Compute  maximum  tensile  stress  -o'  for  the  original  layered 
system,  using  BISAR  program. 

(3)  Using  Newton-Raphson  method  for  determination  of  roots  of 

equation,  compute  E CQ  corresponding  to 


a  -  o'  =0  (5.10a) 

where  o  is  given  by  equation  3.2. 

(4)  Correct  E^  as  computed  from  equation  2.7  by  adding  a  constant 

term 


E  -  E 
comp  comp 


(3.10b) 


The  correction  scheme  corresponds  to  a  translation  of  the  E  -function. 

r  comp 

This  procedure  does  not  affect  the  derivative  of  EcQmp  with  respect  to  design 
variables.  It  is  simple  and  requires  only  one  run  of  the  BISAR  program. 

The  application  of  the  probabilistic  analysis  and  correction  schemes 
will  be  illustrated  through  the  example  runs. 


The  Simulation  Approach 

The  probabilistic  analysis  with  the  simulation  approach  includes  300 
computation  runs  where  the  independent  variables  are  randomly  generated  and 
the  dependent  variables  DF  and  N,  computed  using  equation  2.7,  2.3,  2.1 
and  3.2.  The  number  of  runs  (300)  was  choosen  in  order  to  insure  a  good 
description  of  the  probability  density  distribution  of  both  independent  and 
dependent  variables.  In  the  following  computations,  it  was  assumed  that  the 
independent  variables  are  normally  distributed. 
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It  should  be  noted  that  this  assumption  is  not  mandatory,  and  any  given 
density  distribution  can  be  easily  simulated.  The  generation  of  the  random 
variables  was  made  in  the  following  steps,  for  each  independent  variable: 

(1)  generate  n  ■  12  random  decimal  numbers  with  a  uniform  distribution 
from  0.0  to  1.9,  a  mean  of  0.5  and  a  standard  deviation  of  o  *  1//12. 

(2)  From  the  Central  Limit  Theorem,  the  random  variable  defined  as: 


k  = 


z  ri  +  (M  '  I  77TF3 


i=l 


(3.11a) 


is  a  random  observation  from  an  approximately  normal  distribution  with  mean 
y  and  standard  deviation  o.  Choosing  n  =  12  for  computation  convenience, 
and  y  =  0.0  and  o  =  1.0  for  the  standard  normal  distribution,  leads  to: 

12 

k  =  (  I  r.)-  6  (5.11b) 

i=l  1 


(3)  In  order  to  insure  closeness  to  a  normal  distribution,  any  set 
of  300-k  whose  standard  deviation,  skewness  and  kurtosis  coefficients  is  out¬ 
side  the  range  of  0.95  to  1.05,  -0.15  to  0.15  and  2.5  to  3.5  respectively 
is  disregarded,  and  a  new  number  set  is  generated. 


(4)  generate  the  independent  random  variables  x^  from: 
x.  *  x  (1.  +  k.  *  CV[x)) 


(3.12) 


where 


kj  =  given  by  equation  3.11b 


CV[x]  =  coefficient  of  variation  of  x 

The  above  generated  variables  are  substituted  into  equation  2.7,  2.3, 
2.1  and  3.2  to  give  300  random  variables  of  DF^  and  .  They  are  then 


oaaikasai^iifcitoi  -m 
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analyzed  to  derive  their  mean,  standard  deviation,  and  frequency  distribution. 
Results  of  these  analyses  will  be  presented  in  the  next  paragraph  and  compared 
with  those  of  the  approximate  closed-form  approach. 


It  should  be  noted  that  simulation  approach  uses  the  same  option  for 

correction  as  in  the  approximate  closed-form  one.  In  addition,  whenever 

I12  (the  base/subbase  layer  thickness)  equals  zero  and  ^2^3  ^e  ratio  of 

base  and  subgrade  moduli)  equals  one,  the  E  takes  the  values  of  the 

comp 

subgrade  modulus. 


Run  Examples 

Composite  Modulus  Correction 

Figures  3.1a  to  3.1e  show  the  input  data  and  output  results  of  a 
run  example,  illustrating  the  composite  modulus  correction  methodology.  The 
case  dealt  with  is:  AG-4  loading  and  a  combination  of  pavement  geometry  vari¬ 
ables  corresponding  to  20  percent  deviation  in  the  composite  modulus  evaluation. 
The  exact  computed  stress  from  the  BISAR  program  using  the  original  pavement 
is  479.62  psi.  The  composite  modulus  computed  from  equation  2.7  for  the 
given  pavement  geometry  is  34777  psi,  while  a  value  of  28980  psi  is  necessary 
in  order  to  get  the  correct  stress  of  479.62  from  equation  3.2.  Therefore 
a  correction  factor  of  aQ  *  -5796  psi  is  used  to  reduce  the  composite 
modulus  evaluated  from  equation  2.7.  When  this  correction  factor  is  applied 
in  the  approximate  closed-form  approach,  it  leads  to  exact  values  of  the 
average  stress,  design  factor  and  number  of  coverages. 

The  corrected  modulus  replaces  that  evaluated  from  equation  2.7  in  all 
subsequent  computations.  It  should  be  noted  that  the  correction  has  littlf: 
effect  on  the  derivatives  evaluated  from  the  regression  equations. 
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OUTPUT  RESULTS  FRC3M  SIMULATION  SOLUTION 


When  the  correction  factor  is  used  in  conjunction  with  the  simulation 
approach  -  i.e.,  all  300  composite  moduli  evaluated  are  shifted  by  5796  psi, 
regardless  of  their  generated  value,  the  average  stress  generated  is  483.74 
psi,  higher  than  the  one  computed  from  BISAR  program  by  0.9  percent.  It 
seems  that  the  composite  modulus  correction  scheme  is  adequate  for  reducing 
to  minimum  the  errors  associated  with  the  use  of  regression  equations. 

Comparison  of  the  Approximate  Closed-Form  and  the  Simulation  Approaches 

The  errors  involved  in  the  approximate  closed-form  approach  are  due  to 
the  linear  assumption  used  in  the  Taylor  series  expansion.  Whenever  the 
dependent  variable  is  related  to  the  independent  variables  through  a  non¬ 
linear  function,  some  error  will  result  in  the  evaluation  of  the  average  and 
of  the  variance  of  the  variable.  Therefore,  the  results  should  be  checked 
and  compared  with  other  results  obtained  from  a  different  approach  which  does 
not  assume  linearity.  In  this  paragraph,  the  results  of  the  approximate  closed- 
form  approach  are  compared  to  those  of  the  simulation  scheme  which  assumes  that 
the  independent  variables  are  normally  distributed,  as  achieved  by  the 
generation  scheme  presented  above. 

From  Figures  3. Id  and  3.1e,  it  is  seen  that  the  average  and  the 
coefficient  of  variation  values  of  the  design  factor  are  very  close:  1.46 
and  0.1229  in  the  approximate  closed-form  solution,  and  1.456  and  0.1262  in 
the  simulation  process.  These  values  were  obtained  for  realistic  values  of 
standard  deviations  of  the  independent  variables.  Figures  3.2a  and  3.2b 
show  the  input  data  and  the  relationship  between  reliability  and  number  of 
coverages  of  AG-13.  It  is  seen  that  a  slight  deviation  between  the  results 
exist  at  the  distribution  tails.  The  deviation  which  can  be  attributed  to 
the  number  generation  in  the  simulation  process  is  considered  negligible 
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Figure  3.2b  COMPARISON  OF  APPROXIMATE  CLOSED  FORM  AND  SIMULATION  SOLUTIONS 


from  an  engineering  design  point  of  view. 

Table  3.3  summarizes  study  results  of  analysis  for  AG-2  using  Taylor 

series  expansion.  The  results  are  given  in  terms  of  weighting  factors  of 

the  design  variables  in  the  coefficient  of  variation  computation,  i.e.,  when 

n 

CV2[°F]  -  ifl Cw.  CV[x.))2  (3.13) 

w^  is  the  weighting  factor  of  x^.  Note  that  w^  equals  one  for  *  MR. 

It  is  seen  that  the  weighting  factors  of  the  concrete  layer  parameters  Chj 

and  Ej)  increases  as  the  other  variables  tend  to  increase  the  relative 

stiffness  of  the  concrete  layer.  For  example,  w  increases  as  h.  increases, 

"l  1 

and  as  E 2»  h^  or  E^  decrease.  As  for  the  second  and  subgrade  layers,  their 

weighting  factors  increase  as  the  other  variables  tend  to  increase  the  stress 

in  the  layer.  For  example,  w^  increases  as  E3  increases,  and  as  hj,  h^  or 

3 

E2  decrease.  It  appears  that  the  weighting  factor  expresses  the  functional 
importance  of  the  variable  in  the  design.  It  should  be  noted  that  the  results 
were  derived  without  the  correction  option  of  the  composite  modulus  and  that 
the  weighting  factors  are  slightly  different  when  the  subgrade  layer  is  sub¬ 
divided  in  two  layers  of  equal  moduli.  This  is  attributed  to  the  form  of  the 
composite  modulus  regression  function  used.  It  is  however  stressed  that  the 
deviation  from  the  homogeneous  case  is  negligible. 

In  order  to  compare  the  results  of  Table  3.3  to  those  of  the  simu¬ 
lation  process,  only  one  variable  at  a  time  will  be  varied.  Table  3.4a 
shows  the  results  of  simulating  variation  of  hj  only.  It  is  seen  that: 

(1)  the  coefficient  of  variation  of  the  design  factor  DF  is  proportional  to 
that  of  the  concrete  layer  thickness,  i.e,,  for  the  particular  case 
studied,  the  linearity  assumption  seems  correct;  (2)  the  weighting  factors 


Table  3.3  Summary  of  Results  for  AG-2  using 
Taylor  Series  Expansion 
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Table  3.4a  Results  of  Simulation  for  AG-2  and 
Concrete  Layer  Thickness  Weighting  Factor 
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are  slightly  higher,  by  about  3  percent  in  the  simulation  process  than  in 
the  approximate  closed-form  solution.  This  is  negligible  for  engineering 
practical  purposes. 

Table  3.4b  shows  the  results  of  the  simulation  for  the  subgrade  modulus 
weighting  factor.  It  is  seen  that  (1)  increasing  the  coefficient  of 
variation  of  the  subgrade  modulus  from  0.10  to  0.25  increases  the  coefficient 
of  variation  of  the  design  factor  by  2.60,  a  slight  deviation  from  linearity; 

(2)  the  weighting  factors  are  higher  in  the  simulation  than  in  the  Taylor 
series  expansion  results  (up  to  15  percent) .  While  this  discrepancy  is 
quite  high  for  the  individual  weighting  factor,  its  effect  on  the  coefficient 
of  variation  of  the  design  factor  with  variability  of  all  design  parameters  will 
be  attenuated  in  the  summation  of  the  squared  contribution  (see  equation  3.3). 

It  appears  from  the  above  that  the  results  of  reliability  will  be 
similar  for  the  Taylor  series  expansion  and  simulation  solution.  The  assumed 
linearity  (or  first  order  expansion)  of  the  function  seems  to  be  supported 
by  the  simulation  approach. 

Example 

The  approximate  approach  was  used  to  develop  reliability-number  of 
coverages  curves  for  different  thicknesses  of  the  concrete  layer  (18,  20, 

22  and  24  inch),  see  Figure  3.3a.  Thickness  -  number  of  coverages  curves 
were  then  developed  for  different  reliability  levels  (0.9,  0.8,  0.7,  0.6 
and  0.5)  and  shown  in  Figure  3.3b.  This  run  example  will  serve  to  find 
out  the  reliability  level  used  in  the  current  design.  The  analysis  is  made 
for  two  values  of  percentile  (0.85  and  0.90)  for  computing  the  design 
parameters  in  the  current  "deterministic”  method.  The  results  are  summarized 
in  Table  3.5  where: 


Table  3.4b  Results  of  Simulation  for  AG-2  and 
Subgrade  Modulus  Weighting  Factor 
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1.0 


AND  DIFFERENT  CONCRETE  LAYER  THICKNESSES 


Table  3.5  Summary  of  Pavement  Parameters 
in  the  Current  Design  Method 


Percentile  Value 

0.85 

0.90 

Ka  in  normal  distribution 

1.03 

1.30 

Subgrade  modulus  of  elasticity,  psi 

9,600. 

8,900. 

Base  modulus  of  elasticity,  psi 

40,000. 

37,000. 

Subgrade  modulus  of  reaction,  pci 

100. 

94. 

Composite  subgrade  modulus 

of  reaction,  pci 

200. 

180. 

Modulus  of  rupture  of  concrete,  psi 

630. 

610. 

Required  Thickness,  in 

24.2 

25.1 

a-  u 
o 


,  where:  a  *  value  of  random  width 
u  ■  the  mean  value 
a  *  the  standard  deviation 
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(a)  The  subgrade  and  the  base  modulus  of  elasticity,  and  the  modulus 
of  rupture  are  computed  using: 

x.  *  7.  (1  -  Ka  CV[x.])  (3.14) 

(b)  The  average  values  of  the  variables  and  their  coefficient  of 
variation  are  given  in  Figure  3.2a. 

(c)  The  subgrade  modulus  of  reaction  is  computed  from  (Parker  et  al., 

0): 

log  k  *  (log  Mr  -  1 ,415)/1 .284  (3.15) 

(d)  The  composite  subgrade  modulus  of  reaction  is  derived  from 
Figure  2.3. 

(e)  The  required  thickness  is  calculated  using  design  curves  in  (1) . 

The  thickness  corresponds  to  traffic  area  A  and  9200  coverages. 

Entering  Figure  3.3b  with  N  *=  9200  coverages  and  h  a  24.2  and  25.1  in. 
the  corresponding  reliability  levels  of  0.88  and  0.95  are  derived.  These 
values  give  an  idea  of  the  reliability  levels  used  today  in  the  design 
procedure . 
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Chapter  4 


SUMMARY 

The  design  procedure  of  military  rigid  airfield  pavements  developed 
by  Parker  et  al  (1)  is  expressed  in  this  volume  in  probabilistic  and  re¬ 
liability  terms.  The  procedure,  based  on  the  multi-layer  elastic  theory 
necessitated  further  developments  in  order  to  make  the  analysis  feasible. 
Two  major  investigations  were  conducted:  (1)  Evaluation  of  the  Composite 
Modulus  of  Elasticity  of  the  layered  (subbase/subgrade)  system  underneath 
the  rigid  pavement  and  (2)  Evaluation  of  the  maximum  tensile  stress  at 
the  bottom  of  the  concrete  layer  for  each  of  13  aircraft  types:  USAF 
Classification  AG-1  to  13. 

Valuable  results  were  derived  from  the  investigation  of  the  composite 
modulus.  It  was  found  that: 

(1)  The  composite  modulus  depends  not  only  upon  the  design  parameters 
of  the  layered  (subbase/subgrade)  system  underneath  the  rigid  pavement, 
but  also  upon  the  pavement  thickness  (for  constant  elastic  properties 

of  the  concrete)  and  the  loading  conditions  (number  of  wheels  in  the  gear) . 

(2)  The  effect  of  deviations  in  the  composite  modulus  evaluation  on 
pavement  performance  (expressed  in  terms  of  predicted  allowable  number 
of  coverages)  is  quite  substantial.  The  error  in  the  evaluation  of 
number  of  coverages  is  about  1.2  to  5  times  that  achieved  in  the  evalua¬ 
tion  of  the  composite  modulus.  The  lower  range  corresponds  to  light  load 
aircraft  and  traffic  whereas  the  upper  range  corresponds  to  heavy  load 
aircraft  and  traffic.  This  result  emphasizes  the  need  for  an  accurate 
determination  of  the  composite  modulus. 
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(3)  The  accuracy  achieved  in  a  correlation  between  the  composite 

2 

modulus  and  the  design  parameters  is  excellent  (R  =  0.994  and  SE  =  8%). 
However,  in  the  case  of  heavy  load  aircraft  and  larger  number  of  cover¬ 
ages,  this  (accuracy)  is  insufficient  for  both  design  purposes  and 
probabilistic/reliability  analyses.  Therefore,  a  correction  methodology 
is  developed  to  reduce  the  error  to  minimum. 

A  regression  equation  of  the  maximum  tensile  stress  at  the  bottom 
of  the  concrete  layer  is  derived  for  each  of  13  aircraft  types.  The 
equation  is  highly  accurate,  allowing  its  use  without  any  correction 
for  both  design  purposes  and  probabilistic/reliability  analyses. 

The  above  derivations  of  the  composite  modulus  and  of  the  maximum 
tensile  stress  are  included  in  a  computer  program  for  probabilistic/ 
reliability  analysis  of  rigid  pavements.  Both  the  approximate  closed 
form  (Taylor  series  expansion)  and  the  simulation  solutions  are  imple¬ 
mented.  The  computer  program  can  be  used: 

(1)  In  the  analysis  in  probabilistic/reliability  terms  of  a  specific 
pavement  structure  and  loading  aircraft,  for  a  given  material  properties 
and  variabilities.  The  design  parameters  (means  and  coefficient  of 
variations)  serve  as  input  to  the  computer  program  which  produces  values 
of  number  of  coverages  and  their  corresponding  reliability  levels. 

(2)  In  the  design  of  a  rigid  pavement  at  a  specific  reliability 
level,  given  all  design  parameters.  Several  solutions  of  the  above 
under  (1)  for  different  pavement  structures  must  be  conducted,  and  the 
requested  design  is  derived. 

Several  example  runs  of  the  computer  program  are  presented,  illus¬ 
trating  (a)  the  use  of  the  correction  procedure  in  the  evaluation  of 
the  composite  modulus.  It  is  suggested  to  call  the  procedure  in  the 
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case  of  medium  and  heavy  load  and  traffic,  (b)  the  simultude  of  the 
results  obtained  with  the  approximate  closed  form  and  the  simulation 
solutions  and  (c)  the  interpretation  of  the  current  deterministic 
design  procedure  in  probabilistic/reliability  terms. 
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APPENDIX  I 


USER'S  GUIDE 


Card  1  (12,  19A4) 

1-2  IG  -  aircraft  group  number 

3-78  HED  -  text  for  title 

Card  2  (12) 

1-2  NLA  -  number  of  layers  (=3  when  base-subbase  layer  exists) 

Card  3  (4F10.0) 

I- 10  H1AV  -  average  thickness  of  concrete  layer,  inch. 

II- 20  HI STD  -  coefficient  of  variation  of  concrete  layer  thickness 

21-30  E1AV  -  average  modulus  of  elasticity  of  concrete,  psi. 

31-40  E1STD  -  coefficient  of  variation  of  concrete  modulus  of  elasticity 
Card  4  (4F10.0)  required  if  NLAs3 

I- 10  H2AV  -  average  thickness  of  base-subbase  layer,  inch. 

II- 20  H2STD  -  coefficient  of  variation  of  base-subbase  layer  thickness 

21-30  E2AV  -  modulus  of  elasticity  of  base-subbase  material,  psi. 

31-40  E2STD  -  coefficient  of  variation  of  base-subbase  modulus  of 

elasticity 

Card  5  (2F10.0) 

I- 10  ES  -  subgrade  modulus  of  elasticity,  psi. 

II- 20  ESSTD  -  coefficient  of  variation  of  subgrade  modulus  of  elasticity 

Card  6  (3F10.0) 

I- 10  MR  -  average  concrete  modulus  of  rupture,  psi. 

II- 20  MRSTD  -  c  efficient  of  variation  of  concrete  modulus  of  rupture 

21-30  ROEMR  -  regression  coefficient  between  modulus  of  elasticity  and 
modulus  of  rupture  of  concrete 
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Card  7 
1 

3 

Card  > 

Hot* 


(212) 

2  IRAN  -  flag  i.  0  for  running  Taylor  series  expansion  approach 

>  0  for  running  simulation  approach 

4  IOPT  -  flag  »  0  no  correction  for  composite  modulus  computations 

i  0  run  B1SAR  program  foT  correcting  composite  modulus 
required  if  IOPT  i  0 

Input  data  required  by  BISAR  program 

Several  problems  (Cards  1-7)  can  be  run  in  sequence. 
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APPENDIX  II 


PROGRAM  LISTING 


Appendix  II  is  a  copyrighted  program  listing.  Information  on  the  pro¬ 
gram  can  be  obtained  from  the  authors  of  this  report. 


II-l 


